High-Energy Astrophysics with Neutrino Telescopes 

T. Chiarusi and M. Spurio 
Dipartimento di Fisica dell'Universita di Bologna and INFN 
Viale Berti Pichat 6/2 - 40127 Bologna (Italy) 

Abstract 

Neutrino astrophysics offers new perspectives on the Universe investigation: high 
energy neutrinos, produced by the most energetic phenomena in our Galaxy and in 
the Universe, carry complementary (if not exclusive) information about the cosmos 
with respect to photons. While the small interaction cross section of neutrinos allows 
them to come from the core of astrophysical objects, it is also a drawback, as their 
detection requires a large target mass. This is why it is convenient put huge cosmic 
neutrino detectors in natural locations, like deep underwater or under-ice sites. In 
order to supply for such extremely hostile environmental conditions, new frontiers 
technologies are under development. The aim of this work is to review the motivations 
for high energy neutrino astrophysics, the present status of experimental results and 
the technologies used in underwater/ice Cherenkov experiments, with a special focus 
on the efforts for the construction of a km'^ scale detector in the Mediterranean Sea. 



1 



2 



T. Chiarusi and M. Spurio 



1 Introduction 

Fiat lux. It was written, and scientists never fail to observe new spectacular astrophysical 
discoveries when new experimental techniques on new photons wavelengths are available: 
from the Cosmic Microwave Background observation up to the TeV 7-ray astronomy using 
Imaging Air- Cherenkov Technique. 

Fiat neutrinos, it was never written, and Mr. Pauli itself has feared that this particle 
would never be discovered. Nevertheless, observation of the solar neutrinos and of neutri- 
nos from the supernova 1987A opened up a new observation field. High energy neutrino 
astronomy is a young discipline derived from the fundamental necessity of extending con- 
ventional astronomy beyond the usual electro-magnetic messengers. 

One of the main questions in astroparticle physics is the origin and nature of high- 
energy cosmic rays, CRs (§2). It was discovered at the beginning of the last century that 
energetic charged particles strike the Earth and produce showers of secondary particles in 
the atmosphere. While the energy spectrum of the cosmic rays can be measured up to 
very high energies, their origin remains unclear. There are many indications of the galactic 
origin of the CR bulk (protons and other nuclei up to ~ 10^^ 10^^ eV), although it is not 
possible to directly correlate the CR impinging directions on Earth to astrophysical sources 
since CRs are generally deflected by the galactic magnetic fields. 

Recent advances on ground-based 7-ray astronomy have led to the discovery of more 
than 80 sources of TeV gamma-rays, as described in ^2.2[ No definitive proof still exists 
that galactic CR originate from supernova explosions. Compelling evidences have been 
accumulated by TeV 7-ray telescopes on the possible CR acceleration by some peculiar 
galactic objects (i ]2.4l ). Assuming that at acceleration sites a fraction of the high-energy CRs 
interact with the ambient matter or photon fields, TeV 7-rays are produced by the vr decay 
while neutrinos are produced by charged pion decay. This is the so-called astrophysical 
hadronic model, ^2.3[ which describes the mechanisms which lay behind the production 
of neutrinos and high energy photons from CR interactions. In this framework, the energy 
spectrum of secondary particles follows the same power law of the progenitor CRs and it 
is possible to put constraints to the expected neutrino flux from sources where 7-rays are 
detected. However, almost all observed objects emitting in the TeV 7-ray band are also 
sources of non-thermal X-rays, presumably of synchrotron origin, radiated by multi-TeV 
electrons. Since the same electrons can also radiate TeV 7-rays through inverse Compton 
scattering, this leptonic model represents a competitive process for TeV radiation. Only the 
coincident measurement of neutrinos from the source would give a uncontroversial proof 
of the discovery of the galactic CR acceleration sites. 

The highest energy CRs are probably originated from extragalactic sources, as indicated 
by recent measurements ( i]3.1b . Protons with E > 10^^ eV interact with the cosmic mi- 
crowave background. This effect, known as the Greisen-Zatsepin-Kuzmin cutoff, restricts 
the origin of high energy protons seen on Earth to a small fraction of the Universe, of the 
order of 100 Mpc. The prediction of high energy neutrino sources of extra-galactic origin 
is a direct consequence of the UHE CR observations, ^3.21 This connection between CRs, 
neutrinos and 7-rays can also be used (i ]3.3l ) to put upper bounds on the expected neutrino 
flux from extragalactic sources, since the neutrino energy generation rate will never exceed 
the generation rate of high energy protons. These predictions impose that the scale of the 
neutrino detectors will be of the order of 1 km'^. 

A high-energy neutrino detector behaves as a telescope when the neutrino direction is 
reconstructed with an angular precision of \° or better. This is the case for high energy 
charged current interactions. The accurate measurement of the direction (which can 
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reach 0.2° in water for energies larger than few tens of TeV, © allows the association 
with (known) sources. It allows also the neutrino telescopes to face some of the most funda- 
mental questions on HE physics beyond the standard model, ^ the nature of Dark Matter 
through the indirect search for WIMPs; the study of sub-dominant effects on neutrino os- 
cillations, as those possibly induced by the violation of the Lorentz invariance; the study 
of relic particles (magnetic monopoles, nuclearites) in the cosmic radiation; the coincident 
neutrino emission with gravitational waves. 

The small interaction cross section of neutrinos allows them to come from far away, but 
it is also a drawback, as their detection requires a large target mass. The idea of a neutrino 
telescope based on the detection of the secondary particles produced in neutrino interactions 
was first formulated in the 1960s by Markov 1 1 ). He proposed to install detectors deep in 
a lake or in the sea and to determine the direction of the charged particles with the help of 
Cherenkov radiation. As we will show in Jj6l starting from the Markov idea and from the 
present knowledge of TeV 7-rays sources, the challenge to detect galactic neutrinos is open 
for a kilometer-scale apparatus. We will use the fact that high-energy muons retain informa- 
tion on the direction of the incident neutrino and can pass through several kilometers of ice 
or water, ^5.31 Along their trajectory, the muons emit Cherenkov light. From the measured 
arrival time of the Cherenkov light ^5.4[ the direction of the muon can be determined. This 
process is referred to as muon track reconstruction. We will also derive in a simple way that 
the number of optical sensors required to reconstruct muon tracks is of the order of 5000. 

Neutrino production in astrophysical sites through -k or K decay leads to a flavor ratio at 
sources of fg : t'^ : z^r = 1 : 2 : 0, which is changed by the neutrino oscillation mechanism 
to 1 : 1 : 1 on Earth. As discussed in ^5. H and ^5.2[ the measurements of showers induced 
by very and ultra high energy and v^- is another very important challenge for large volume 
neutrino detectors, although the neutrino direction measurement is poorer for these flavors 
with respect to the t'^ channel. The extragalactic CRs-neutrinos connection [,2J sets also the 
scale of the detectors to 1 km^. 

The properties of water and ice connected to the possibility of detecting high energy 
neutrinos are discussed in f|2l The pioneering project for the construction of an under- 
water neutrino telescope was due to the DUMAND collaboration 111, which attempted to 
deploy a detector off the coast of Hawaii in the 1980s. At the time technology was not ad- 
vanced enough to overcome these challenges and the project was cancelled. In parallel, the 
BAIKAL collaboration [4] started to work in order to realize a workable detector systems 
under the surface of the Baikal lake (©. 

Regarding deep ice, a major step towards the construction of a large neutrino detector 
(see Sj9ll is due to the AMANDA collaboration Q. AMANDA deployed and operated the 
optical sensors in the ice layer of the Antarctic starting from 1993. After the completion 
of the detector in 2000, the AMANDA collaboration proceeded with the construction of a 
much larger apparatus, IceCube. 59 of the 80 scheduled strings (April 2009) are already 
buried in the ice. Completion of this detector is expected to be in 201 1. 

In water, the pioneering DUMAND experience is being continued in the Mediterranean 
Sea by the ANTARES [6|, NEMO [7] and NESTOR [8] collaborations, which demonstrated 
the detection technique (see ^flOl ). The ANTARES collaboration has completed (May 2008) 
the construction of the largest neutrino telescope {r^ 0.1 km^) in the Northern hemisphere. 
The ANTARES detector currently take data. These projects have lead to a common design 
study towards the construction of a km^-scale detector in the Mediterranean Sea (^T\\. 
KM3NeT |9| is an European deep-sea research infrastructure, which will host a neutrino 
telescope with a volume of at least one cubic kilometer at the bottom of the Mediterranean 
Sea that will open a new window on the Universe. 
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2 The connection among primary Cosmic Rays, 7-rays and neu- 
trino. Our Galaxy. 



2.1 Primary Cosmic Rays 

Cosmic Rays (CRs) are mainly high energy protons (Fig. [B and heavier nuclei which are 
constantly hitting the upper shells of the Earth's atmosphere. The energy spectrum spans 
from ~ 10^ eV to more than lO^'' eV, is of non-thermal origin and follows a broken power- 
law of the form: 
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Figure 1: Cosmic Ray spectrum from 10^ to 10^^ eV as measured on Earth, from l\10\l . Note 
that the vertical scale has been multiplied by E"^. On the low-energy domain, when the mea- 
surements are available, the contribution of protons, electrons, positrons and antiprotons it 
is also reported. See HlOV for the reference to the experiments. 
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Direct or indirect techniques are used to measure the CR spectrum. The measured 
power law spectrum of CRs (eq. [B is characterized by an index a = 2.7 up to energies of 
roughly 3 x 10^^ eV. Beyond 3 x 10^^ eV, the index becomes a = 3.1. This feature in the 
energy spectrum is known as the knee. There is no consensus on a preferred accelerator 
model for energies above the knee up to 10^^ eV, where there is a flattening in the spectrum, 
denoted as the ankle. 

The highest CRs exceed even 10^'^ eV. After the ankle, it is generally assumed that CR 
sources are of extragalactic origin. The experimental search for sources of these ultra high 
energy CRs is recently entered a hot phase. Detailed reviews of the theory and measurement 
of the primary CR spectrum can be found in lfTTllT2l[T3ll . 

Up to energies of ~ 10^'^ eV, the CR spectrum is directly measured above the atmo- 
sphere. Stratospheric balloons or satellites have provided the most relevant information 
about the composition of CRs in the Galaxy and have contributed to establish the standard 
model of galactic CRs. Measurements show that ~ 90% are protons, 9% are Helium 
nuclei and ~ 1% are heavier nuclei. 

In this energy range, the mechanism responsible for the acceleration of particles is the 
Fermi mechanism [14, 15 1. This mechanism explains the particle acceleration by iterative 
scattering processes of charged particles in a shock-wave. These shock-waves are originated 
in environments of exceptional disruptive events, like stellar gravitational collapses. In 
each scattering process, a particle with energy E gets an energy gain of AE ^ (3E, where 
P ~ 10~^. Due to the magnetic fields confinement, the scattered particles are trapped inside 
the acceleration region and they have a small probability to escape. This iterative process 
of acceleration is a very appealing scheme for the origin of CRs, since it naturally explains 
the power law tendency in the spectrum. 

Supernova remnants (SNR) in the Galaxy are the most accredited site of acceleration 
of CRs up to the knee |16|, although this theory is not free from some difficulties |17|. 
The Fermi mechanism in the SNR 1 18|, predicts a power law differential energy spectrum 
~ i?^^ and fits correctly to the energy power involved in the galactic cosmic rays of ~ 
5 X 10^° erg/s. 

The measured spectral index (a ~ 2.7) is steeper than the expected spectrum near the 
sources, because of the energy dependence of the CR diffusion out of the Galaxy, as ex- 
plained by the so called leaky box |[T9l . In the leaky box model, particles are confined by 
galactic magnetic fields (B ~ 3/iG) and have a small probability to escape. The gyro- 
magnetic radius for a particle with charge Z, energy E, in a magnetic field B is i? ~ 
During propagation, high energy particles (at a fixed value of Ze) have larger probabil- 
ity to escape from the Galaxy due to their larger gyromagnetic radii. As a consequence, 
an energy-dependent diffusion probability P can be defined. P is experimentally estimated 
through the measurement of the ratio between light isotopes produced by spallation of heav- 
ier nuclei. It was found that P{E) ~ i?"^, with the diffusion exponent ao ~ 0.6 [11]. 
The differential CR flux at the sources can be estimated as the convolution of the measured 
spectrum ([T]) and the CR escape probability P: 
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with acB. = a — ~ 2, as predicted by the Fermi model. 

The knee of the CR spectrum is still an open question and different models have been 
proposed to explain this feature |[20l . Some models invoke astrophysical reasons: due to the 
iterative scattering processes involved in the acceleration sites, a maximum energy for the 
CRs is expected. This maximum energy depends on the nucleus charge Ze, and this leads 
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to the prediction of a different energy cutoff for every nucleus type. As a consequence, CRs 
composition is expected to be proton-rich before the knee, and iron-rich after. Other more 
exotic models try to explain the steepening in the CR flux, for instance the hypothesis of 
new particle processes in the atmosphere 121]. 

Above ^ 10^^ eV, CR measurements are only accessible from ground detection infras- 
tructures. The showers of seconddary particles created by interaction of primary CRs in 
the atmosphere are distributed in a large area, enough to be detected by detector arrays. 
The energy region around the knee has been explored by different experiments, as for in- 
stance KASCADE [22|. Although the experimental techniques are very difficult and have 
poor resolution, observations of this region of the energy spectrum seem to indicate that the 
average mass of CRs increases when passing the knee. 

The SNR models cannot explain the CRs flux above 10^^ eV, but there is no con- 
sensus on a preferred accelerator model up to 10^^ eV. CRs can be accelerated beyond the 
knee if, for instance, the central core of the supernova hosts a rotating neutron star. Already 
accelerated particles can also suffer additional acceleration due to the neutron star strong 
variable magnetic field. The maximum energy cannot exceed ~ 10^^ eV. 

2.2 High energy 7-rays 

Some galactic accelerators must exist to explain the presence of CRs with energies up to 
the ankle. These sources can be potentially interesting for a neutrino telescope. Apart from 
details, it is expected that galactic accelerators are related to the final stage of the evolution 
of massive, bright and relatively short-lived stellar progenitors. 

Due to the influence of galactic magnetic fields, charged particles do not point to the 
sources. Neutral particles (gamma-rays and neutrinos) do not suffer the effect of magnetic 
fields: they represent the decay products of accelerated charged particles but cannot be 
directly accelerated. 

Photons in the MeV-GeV energy range were detected by the Energetic Gamma-ray 
Experiment Telescope (EGRET) f23\ on board of the CGRO satellite in the 1990s. The last 
EGRET catalogue contains 271 detections with high significance, from which 170 are not 
identified yet. 

Following its launch in June 2008 1*241, the Fermi Gamma-ray Space Telescope (Fermi) 
began a sky survey in August. Its Large Area Telescope (LAT) has produced, in 3 months, a 
deeper and better-resolved map of the 7-ray sky than any previous space mission. The initial 
result for energies above 100 MeV 125 ] regards the 205 most significant 7-ray sources, 
which are the best-characterized and best-localized ones. Most of them are in the galactic 
plane, and were associated with known pulsars. 

Gamma-rays above 100 GeV are detected on ground, using the Imaging Air-Cherenkov 
Technique (lACT). High-energy 7-rays are absorbed when reaching the Earth atmosphere, 
and the absorption process proceeds by creation of a cascade (shower) of high-energy rela- 
tivistic secondary particles. These emit Cherenkov radiation, at a characteristic angle in the 
visible and UV range, which passes through the atmosphere. As a result of Cherenkov light 
collection by a suitable mirror in a camera, the showers can be observed on the surface of 
the Earth. 

The pioneering ground based 7-ray experiment was built by the Whipple collaboration 
||26l. During the last decade, several ground-based 7-ray detectors were developed, both in 
the North |27 | and South [28, 29 1 Earth hemisphere. At present, the new generation appa- 
ratus are the H.E.S.S. 1301 and VERITAS |31| telescope arrays and the MAGIC telescopes 
||32il . A full and detailed review of VHE astrophysics with the ground-based 7-ray detectors 
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can be found in |[33l[34l . These lACT telescopes have produced a catalogue of TeV 7-ray 
sources. A (partial) sky map can be seen in Fig. [2l Of particular interest (mainly for a 
neutrino detector placed in the North hemisphere) is the great population of new TeV 7-ray 
sources in the galactic centre region discovered by the H.E.S.S. telescope. A list of more 
than 70 galactic and extra-galactic sources is in ||331 . 

Both electrons and protons can be accelerated by astrophysical objects. We refer re- 
spectively to a leptonic model ll33l [34l when electrons are accelerated, and to a hadronic 
model [35 1 when protons or other nuclei are accelerated. The most important process which 
produces high energy 7-rays in the leptonic model is the Inverse Compton (IC) scattering. 
IC 7-rays are produced in the interactions of energetic electrons with ambient background 
photon fields: the CMB, and the diffuse galactic radiation of star light. This process is 
very efficient for producing 7-rays, since low energy photons are found in all astrophys- 
ical objects. Multi-TeV electrons producing 7-rays of TeV energies via IC, produce also 
synchrotron radiation in the X-ray band as well fj6^|. Therefore, measurements of the syn- 
chrotron X-ray flux from a source is a signal that the accompanying 7-rays are likely pro- 
duced by leptonic processes. 

Both models, the leptonic model and the hadronic model fST] could provide an adequate 
description of the present experimental situation. If high energy photons are produced in 
the hadronic models, high energy neutrinos will be produced as well. Most of observed TeV 
7-ray galactic sources have a power law energy spectrum E^"-'' , where a-y ~ 2.0 2.5. 
The values of the spectral index are very close to the expected spectral index of CR sources, 
acR- This lead to the conclusion that sources of TeV 7-rays can also be the sources of 
galactic CRs. 
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Some of the most promising candidate neutrino sources in our Galaxy are extremely in- 
teresting, due to the recent results from 7-ray detectors. A neutrino telescope in the North- 
ern hemisphere (as a detector in the Mediterranean sea) is looking at the same Southern 
field-of-view as the H.E.S.S. and CANGAROO Imaging Air Cherenkov telescopes, while 
the neutrino telescope in the South Pole is looking at the Northern sky. 

2.3 TeV 7-rays and neutrinos from hadronic processes 

The astrophysical production of high energy neutrinos is mainly supposed via the decay of 
charged pions in the beam dump of energetic protons in dense matter or photons field. 

Accelerated protons will interact in the surroundings of the CRs emitter with photons 
predominantly via the A+ resonance: 

p + 7^A+^7r+ + n (3) 

Protons will interact also with ambient matter (protons, neutrons and nuclei), giving 
rise to the production of charged and neutral mesons. The relationship between sources of 
VHE 7-ray (Ey > 100 MeV) and neutrinos is the meson-decay channel. Neutral mesons 
decay in photons (observed at Earth as 7-rays): 

TT° 77 (4) 

while charged mesons decay in neutrinos: 



^ /u+ ^ 17^ + z^e + e+ 

^ fi~ ^ i'f^+Ve + e~ (5) 

Therefore, in the framework of the hadronic model and in the case of transparent 
sources, the energy escaping from the source is distributed between CRs, 7-rays and neutri- 
nos. A transparent source is defined as a source of a much larger size that the proton mean 
free path, but smaller than the meson decay length. For these sources, protons have large 
probability of interacting once, and most secondary mesons can decay. 

Because the mechanisms that produce cosmic rays can produce also neutrinos and high- 
energy photons (from eqs. |4]and|5]l, candidates for neutrino sources are in general also 7-ray 
sources. In the hadronic model there is a strong relationship between the spectral index of 
the CR energy spectrum E~"^^ , and the one of 7-rays and neutrinos. It is expected |[38l 
that near the sources, the spectral index of secondary 7 and u should be almost identical 
to that of parent primary CRs: acR ~ a^/ ~ a^. Hence 7-ray measurements give crucial 
information about primary CRs, and they constrain (see ^2 41) the expected neutrino flux. 

2.4 Prediction of HE neutrino flux from astrophysical sources 

Here, we present some proposed mechanisms for the production of cosmic high energy 
neutrinos. Some of them seem to be guaranteed, since complementary observations of TeV 
7-rays can hardly be explained by leptonic models alone. The expected neutrino fluxes at 
Earth, however, are uncertain and predictions differ in some cases up to orders of magnitude. 
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2.4.1 Shell-type supernova remnants 

Particles can be accelerated in the supernova remnants (SNRs) via the Fermi mechanism. 
If the final product of the SN is a neutron star, already accelerated particles can gain addi- 
tional energy, due to the neutron star strong variable magnetic field. Shell-type SNRs are 
considered to be the most likely sites of galactic CR acceleration, hypothesis supported by 
recent observations from the TeV 7-ray lACT. 

Of particular interest is the supernova remnant in the Vela Jr. (RX J0852.0-4622). This 
SNR is one of the brightest objects in the southern TeV sky. Recent observations of 7- 
rays exceeding 10 TeV in the spectrum of this SNR by H.E.S.S. |39| have strengthened the 
hypothesis that the hadronic acceleration is the process that is needed to explain the hard and 
intense TeV 7-ray spectrum. H.E.S.S. has observed that the 7-ray TeV emission originates 
from several separated parts of a region of apparent size of ~ 2°. The angular resolution of 
neutrino telescopes for the channel is much better than 2° (©. The expected neutrino- 
induced muon rate leads, in some calculations BOl . to encouraging results for a 1 km'^ 
detector in the Mediterranean sea. 
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Figure 3 : Measured gamma-ray flux from RX Jl 71 3. 7-3946 and estimated neutrino flux 
with their error bands. The atmospheric neutrino flux ( which represents the background) is 
integrated over the search window and averaged over one day. 



A second important source is the SNR RX J17 13.7-3946, which has been the sub- 
ject of large debates about the nature of the process (leptonic or hadronic) that originates 
its gamma-ray spectrum Ell. RX J1713.7-3946 was first observed by the CANGAROO 
experiment which firstly claimed a leptonic origin [42]. Successive observations with 
CANGAROO-II H3l disfavour purely electromagnetic processes as the only source of the 
observed 7-ray spectrum. Neutrino flux calculations based on this result have predicted 
large event rates, also in neutrino telescopes with size smaller than 1 km'^ li44l l45l . This 
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source has successively been observed with higher statistics by the H.E.S.S. telescope B6l . 
supporting the hadronic origin. The measured spectrum deviates from a pure power law 
spectrum. It can be reasonably well described by a power law with an exponential cutoff, 

d^^/dE^ = Io{E^/l TeV)"'-' exp{-E^/Ec), where the cutoff parameter is = 12 ± 2 
TeV |33|. 

The neutrino flux calculation (shown in Fig. [3]) based on the H.E.S.S. result, with the 
exponential cutoff in spectrum and other assumptions, lead to the prediction that the source 
should be marginally detectable in a kilometer-scale Mediterranean detector. This result 
strongly depends on the assumed cutoff value. Without cutoff, the event rate increases by a 
significant factor, making these sources easily accessible to neutrino telescopes. 

2.4.2 Pulsar wind nebulae (PWNe) 

PWNe are also called Crab-like remnants, since they resemble the Crab Nebula ( ^6.11 ). 
which is the youngest and most energetic known object of this type. PWNe differ from the 
shell-type SNRs because there is a pulsar in the center which blows out equatorial winds 
and, in some cases, jets of very fast-moving material into the nebula. The radio, optical 
and X-ray observations suggest a synchrotron origin for these emissions. H.E.S.S. has also 
detected TeV 7-ray emission from the Vela PWN, named Vela X. This emission is likely to 
be produced by the inverse Compton mechanism. The possibility of a hadronic origin for 
the observed 7-ray spectrum, with the consequent flux of neutrinos, was also considered 

izi. 

The neutrino flux calculated for a few PWNe in the framework on a hadronic production 
of the observed TeV 7-rays (such as the Crab, the Vela X, the PWN around PSR1706-44 
and the nebula surrounding PSR1509-58) agrees with the conclusions that all these PWNe 
could be detected by a kilometer-scale neutrino telescope. For instance, from 6 to 12 events 
are predicted in 1 y, with 1 background event due to the atmospheric neutrinos P8l . Others 
||38l . which assume an exponential cutoff in the energy spectrum, give more pessimistic 
results (10 events/5 y from the Vela X source, with 4.6 background events). 

2.4.3 The galactic Centre (GC) 

The galactic centre is probably the most interesting region of our Galaxy, also regarding 
the emission of neutrinos. It is specially appealing for a Mediterranean neutrino telescope 
since it is within the sky view of a telescope located at such latitude. The interest in it has 
increased after the recent discoveries of H.E.S.S.. 

Early H.E.S.S. observations of the GC region showed a point-like source at the gravita- 
tional centre of the Galaxy (HESS J 1745 -290 \A9l ) coincident with the supermassive black 
hole Sagittarius A* and the SNR Sgr A East. In 2004, a more sensitive campaign revealed 
a second source, the PWN G 0.9-1-0.1 (50]- 

Thanks to the good sensitivity of the H.E.S.S. telescope, it is possible to subtract the 
GC sources and search for the diffuse 7- ray emission which spans the galactic coordinates 
|/| < 0.8°, |6| < 0.3". This diffuse emission of 7-ray with energies greater than 100 GeV 
is correlated with a complex of giant molecular clouds in the central 200 pc of the Milky 
Way |[5T1 . The measured 7-ray spectrum in the GC region is well described by a power 
law with index of ~ 2.3. The photon index of the 7-rays, which closely traces back the 
spectral index of the CR, indicates in the galactic centre a local CR spectrum that is much 
harder and denser than that measured on Earth, as shown in Fig. |4] Thus it is likely that 
an additional component of the CR population is present in the galactic centre, above the 



High-Energy Astrophysics with Neutrino Telescopes 



11 



diffuse CR concentration which fills the whole Galaxy. The fact that cosmic accelerators 
are very close to the GC, and therefore the possibility of neglecting the CR diffusion loss 
due to propagation (see eq. |2ll, gives a natural explanation for the harder observed spectrum, 
which is closer to the intrinsic value of the CR spectral index. In 151] it is suggested that the 
central source HESS J1745-290 is likely to be the source of these CR protons (and thus of 
neutrinos), with two candidates for CR accelerations in its proximity: the SNR Sgr A East 
(estimated age around 10^ yrs), and the black hole Sgr A*. 



GC region 




Energy (TeV) 



Figure 4: H.E.S.S. [51 / measurement of j-ray flux per unit solid angle in the Sgr B region, 
near the galactic center (open circle data points). The measured spectrum in the galactic 
region \l\ < 0.8", \b\ < 0.3° is shown using fij.ll circles. These data can be described by 
a power law with spectral index ~ 2.3. In comparison, the expected ^ -ray flux assuming 
a CR spectrum as measured in the solar neighborhood is shown as a shaded region. The 
measured ^ -ray flux (>1 TeV) implies a high-energy cosmic ray density which is 4 to 10 
times higher than our solar neighborhood value. The spectrum of the source HESS Jl 745- 
290 is also shown for comparison. 



2.4.4 Microquasars 

Microquasars are galactic X-ray binary systems, which exhibit relativistic radio jets, ob- 
served in the radio band [52 1. The name is due to the fact that they result morphologically 
similar to the AGN, since the presence of jets makes them similar to small quasars. This 
resemblance could be more than morphological: the physical processes that govern the for- 
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mation of the accretion disk and the plasma ejection in microquasars are probably the same 
ones as in large AGN. 

Microquasars have been proposed as galactic acceleration sites of charged particles up 
to £' ~ 10^^ eV. The hypothesis was strengthened by the discovery of the presence of 
relativistic nuclei in microquasars jets like those of SS 433. This was inferred from the 
observation of iron X-ray lines [53 ] . 

Two microquasars, LS I +61 303 and LS 5039, have been detected as 7-ray sources 
above 100 MeV and listed in the third EGRET Catalogue. They are also detected in the 
TeV energy range l !54ll55l . 

There is yet uncertainty as to what kind of compact object lies in LS I +61 303 (observed 
by the MAGIC telescope). Recently, a multi wavelength campaign including the MAGIC 
telescope, XMM -Newton, and Swift was conducted during 60% of an orbit in 2007. A 
simultaneous outburst at X-ray and TeV 7-ray bands, with the peak at phase 0.62 and a 
similar shape at both wavelengths, gives conclusive indication of variability also in the 7- 
ray emission. The X-ray over TeV 7-ray flux ratio favors leptonic models [56]. Because 
the source is located in the Northern sky, it is specially appealing for a neutrino telescope 
located in the Southern hemisphere as IceCube, which will be able to detect (or rule out) 
neutrinos coming from this source ll57l . 

Microquasar LS 5039 (detected by H.E.S.S. in the Southern sky) has features similar to 
LS I +61 303, and the observed flux still does not allow an unequivocal conclusion about the 
variability of the source. Different astrophysical scenarios have been proposed to explain 
the TeV 7-ray emission, which involve leptonic and/or hadronic interactions. In particular, 
the leptonic model is strongly disfavored in 1.58,1 . and it is expected that LS 5039 could 
produce between 0.1 0.3 events/year in a detector like ANTARES (see ^10.11 ). The event 
rate depends on the assumed neutrino spectrum (power law with index ranging from 1.5 to 
2.0), and two energy cutoff (Ec = 10 TeV and 100 TeV [58]). The expected rate is ~ 25 
times higher for a 1 km^ detector in the Mediterranean sea. 

Other microquasars were considered in |59|. The best candidates as neutrino sources 
are the steady microquasars SS433 and GX339-4. Assuming reasonable scenarios for TeV 
neutrino production, a 1 km'^-scale neutrino telescope in the Mediterranean sea could iden- 
tify microquasars in a few years of data taking, with the possibility of a 5a level detection. 
In case of no-observation, the result would strongly constrain the neutrino production mod- 
els and the source parameters. 

2.4.5 Neutrinos from the galactic plane 

In addition to stars, the Galaxy contains interstellar thermal gas, magnetic fields and CRs 
which have roughly the same energy density. The inhomogeneous magnetic fields confine 
the CRs within the Galaxy. Hadronic interactions of CRs with the interstellar material 
produce a diffuse flux of 7-rays and neutrinos (expected to be equal, within a factor of ~ 
2). The fluency at Earth is expected to be correlated to the gas column density in the Galaxy: 
the largest emission is expected from directions along the line of sight which intersects most 
matter. 

Recently, the MILAGRO collaboration has reported the detection of extended multi- 
TeV gamma emission from the Cygnus region 1601 . which is well correlated to the gas 
density and strongly supports the hadronic origin of the radiation. The MILAGRO obser- 
vations are inconsistent with an extrapolation of the EGRET flux measured at energies of 
tens of GeV This supports the hypothesis that in some areas of the galactic disk the CR 
spectrum might be significantly harder that the local one. 
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With the assumption that the observed 7-ray emission comes from hadronic processes, 
it is possible to obtain an upper limit on the diffuse flux of neutrinos from the galactic 
plane. The KM3NeT consortium [61] has made an estimate of the neutrino flux from the 
inner Galaxy, assuming that the emission is equal to that observed from the direction of 
the Cygnus region. The expected signal rate for a km'^ neutrino telescope located in the 
Mediterranean Sea is between 4 and 9 events/year for the soft (index a = 2.55) and hard 
(a = 2.10) spectrum respectively, with an atmospheric neutrino background of about 12 
events per year. 

2.4.6 Unknowns 

In addition to SNR, PWNe and microquasars, there are other theoretical environments in 
which hadronic acceleration processes could take place with production of a neutrino flux. 
For instance, neutron stars in binary systems and magnetars [|62l might be sources of an 
observable neutrino flux. 

New improvements in the GeV- TeV scale 7-ray astronomy are expected in the next 
years. In particular practically all the lACT telescopes are improving their apparatus. News 
are also expected from the ARGO |[63l and MILAGRO [64] large field of view observato- 
ries. Finally, it is also worth remarking that a non-negligible number of VHE 7-ray sources 
detected by H.E.S.S. do not have a known counterpart in other wavelengths. The origin of 
such sources is a theoretical challenge in which neutrino astronomy may yield some insight. 

Although not certainly inspired by neutrino astronomy, it is interesting to quote this 
sentence from the former US Secretary of Defense, Donald Rumsfeld. The sentence is the 
exact words as taken from the official transcripts on the Defense Department Web site [ 65 1 : 
The Unknown. As we know, / There are known knowns. / There are things we know we 
know. / We also know / There are known unknowns. / That is to say / We know there are 
some things/ We do not know. /But there are also unknown unknowns, / The ones we don't 
know / We don 't know. 

3 The connection among extragalactic sources of primary Cos- 
mic Rays, 7 rays and neutrino. 

3.1 IMeasurements of the UHECR 

The CR flux above 10^^ eV, still dominated by protons and nuclei |[66l . is one particle per 
kilometer square per year per stereoradian. It has long been assumed (6T\ that ultra high 
energy cosmic rays (UHECR) are extragalactic in origin [68], and can be detected only by 
very large ground-based installations. Therefore, the structure in the CR spectrum above 
~ 10^^ eV (the ankle) is usually associated with the appearance of this flatter contribution 
of extra-galactic CRs. In fact, above the ankle the gyroradius of a proton in the galactic 
magnetic field exceeds the size of the Galaxy disk (300 pc). 

Fig. [5] [69] shows a diagram first produced by Hillas (1984). Hillas derived the maxi- 
mum energy which a particle of a given charge can reach, independently of the acceleration 
mechanism. It was obtained from the simple argument that the Larmor radius of the par- 
ticle should be smaller than the size Rkpc (in kpc) of the acceleration region. This energy 
-E(EeV) (in units of 10^^ eV) is given by: 



EiEeV) ~ pZB^aRkpc 



(6) 
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Figure 5: The Hillas diagram (drawn by Murat Boratav). Acceleration of cosmic rays up to 
a given energy requires magnetic fields and sizes above the respective line. Some sources 
candidates are still controversial (1 EeV=l0^^eV, 1 ZeV=10^^eV). 



where /? is the velocity of the shock wave in the Fermi model or in any other acceleration 
mechanism. Fig. |5] gives the relation between the dimensions of the astrophysical objects 
and the magnetic fields needed to contain the accelerating particle, in order that protons can 
reach up to 10^'' eV (dashed line) or 10^^ eV (upper full line), and iron nuclei up to lO^'^ 
eV (lower full line). As can be seen from the Hillas plot, plausible acceleration sites may 
be the radio lobes or hot spots of powerful active galaxies. 

The search for UHECR sources must take into account another effect, the Greisen- 
Zatsepin-Kuzmin cutoff (GZK) fTOl ItTI . which imposes a theoretical upper limit on the 
energy of cosmic rays from distant sources. Above a threshold of few 10^^ eV, protons 
interact with the 2.7° K cosmic microwave background radiation (CMB) and lose energy 
through the resonant pion production of eq. [3] Due to the GZK cutoff, protons above that 
threshold cannot travel distances further than few tens of Mpc. 

From the astrophysical point of view, this cutoff is very important because it limits the 
existence of standard astrophysical UHECR emitters inside our local super-cluster of galax- 
ies. The GZK cutoff has stimulated important debate, since there were two contradictory 
measurements in the region between lO" ^ lO^o eV, made by the AGASA and by the 
High Resolution Fly's Eye (HiRes) |73| experiments. 

Nowadays, the largest experiment is the Auger Observatory 1741, which combines the 
measurement of extensive air showers and light fluorescence detection. Auger has pub- 
lished |75 1 the result of the first data set, rejecting the hypothesis that the cosmic ray spec- 
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trum continues in the form of a power law above 10^^'® eV with 6 sigma significance. In 
addition, Auger has reported the first hints of association of CRs with E > 6x 10^^ eV and 
nearby (less than 100 Mpc) concentration of matter and AGN [76]. Although its statistical 
significance is still limited, the results suggest that regions of matter with AGN can be the 
source candidates for UHECR acceleration. 

3.2 Extra- galactic CR and u sources 

The prediction of high energy neutrino sources of extra-galactic origin is a direct conse- 
quence of the CR observations. As for the origin of UHE Cosmic Rays, Active Galactic 
Nuclei (AGN) are the principal candidates as neutrino sources. Other potentially promising 
particle accelerators are 7-ray bursts (GRBs). Here we consider these two astrophysical 
classes of objects, with a particular attention to the possible neutrino production mecha- 
nism. Finally, radio observation of starburst galaxies has motivated the idea of the existence 
of hidden sources of CR. These sources can represent pure high energy neutrino injectors, 
and some predictions are presented. 

Extra-galactic sources are very far and the possibility of a individual discovery in a 
km^ scale neutrino telescope is expected only in particular theoretical models, or using the 
source stacking methods: it is a combined analysis for different classes of objects which 
enhance the neutrinos detection probability [77j|. 

An alternative way to prove the existence of extragalactic neutrino sources is through 
the measurement of the cumulative flux in the whole sky. The only way to detect this diffuse 
flux of high energy neutrinos is looking for an excess oi^ high energy events in the measured 
energy spectrum induced by atmospheric neutrinos. 

Theoretical models constrain the neutrino diffuse flux, ^3.3[ These upper bounds are 
derived from the observation of the diffuse fluxes of 7-rays and UHECR. One of them 
(the Waxman-Bahcall, shortened as W&B) is used as the reference limit to the predicted 
neutrino flux coming from different extra-galactic sources. 

In addition to neutrinos generated by high energy cosmic accelerators, there are high 
energy neutrinos induced by the propagation of CRs in the local Universe [781. The sub- 
sequent pions decay will produce a neutrino flux (called GZK or cosmological neutrinos) 
similar to the W&B bound above 5 x 10^^ eV Ii79.l . since neutrinos carry approximately 5% 
of the proton energy. 

3.2.1 Active Galactic Nuclei (AGN) 

Active Galactic Nuclei (or AGN) are galaxies with a very bright core of emission embed- 
ded in their centre, where a supermassive black hole (10° ^ 10^ solar masses) is probably 
present. As outlined in ^3.1[ the Auger observatory has reported the first hints of correla- 
tion between CR directions and nearby concentrations of matter in which AGN are present. 
This measurement (although still controversial) suggests that AGN are the most promising 
candidates for UHECR emission. 

The supermassive black hole in the centre of AGN would attract material onto it, re- 
leasing a large amount of gravitational energy. According to some models [80], the energy 
rate generated with this mechanism by the brightest AGNs can be L > 10^^ erg s~^. Early 
models fSTl [82l Wi \. postulating the hadronic acceleration in the AGN cores, predicted a 
production of secondary neutrinos well above the W&B upper limit, and the prediction 
from some of these models has been experimentally disproved by AMANDA 111121 . More 
recent models [ 84 1 predict fluxes close to the W&B bound. For instance, a prediction has 
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recently been carried out for the Centaurus A Galaxy, which is only 3 Mpc away. In ll85l the 
estimate neutrino flux from hadronic process is E'^d^^/dE < 5 x 10^ TeV^^ cm~^s~^. 
By varying the power-law indices between -2.0 and -3.0, they obtaine d be tween 0.8 and 
0.02 events/year for a generic neutrino detector of effective muon area ( ^6.31) of ~ 1 km^. 

A particular class of AGN (called hlazars) has their jet axis aligned close to the line of 
sight of the observer. Blazars present the best chance of detecting AGNs as individual point 
sources of neutrinos because of a significant flux enhancement in the jet through Doppler 
broadening. Blazars exhibit non-thermal continuum emission from radio to VHE frequen- 
cies and are highly variable, with fluxes varying by factors of around 10 over timescales 
from less than 1 hour to months. The third EGRET catalog |[86]l contains a list of 66 
blazars, plus 27 additional candidates, and 119 are in the recent Fermi LAT bright gamma- 
rays source list [25|; an increasing population of TeV blazars at higher redshifts is being 
detected by the latest generation of 7-ray lACT; so far 1 8 blazars have been discovered over 
a range of red-shifts from 0.03 to >0.3 [34|. 

In hadronic blazar models, the TeV radiation is produced by highly relativistic baryons 
in jets interacting with radiation fields and the ambient matter. Owing to the low mat- 
ter density in relativistic jets, photo-production of pions is commonly believed to be the 
most important energy loss channel, followed by proton synchrotron radiation. The 7 rays 
from neutral pion decay induce electromagnetic cascades, disrupting the strict neutrino-to- 
gamma-ray ratio of pion decay kinematics for the emerging radiation. Another important ef- 
fect to take into account is that the observed TeV 7-ray spectrum from extragalactic sources 
is steepened due to absorption by the Extragalactic Background Light (EBL). In the case of 
a distant blazar, such as lESl 101 at z=0.186, the observed spectral index of 2.9 is estimated 
to correspond to a spectral index as hard as 1.5 near the source ll87l . Neutrinos, however, are 
unaffected by the EBL. As a consequence of the effective hardening of the spectrum, some 
TeV-7 bright blazars, in some models, are expected to produce fluxes exceeding the at- 
mospheric neutrino background in a cubic kilometer neutrino telescope. H.E.S.S. recently 
reported also highly variable emission from the blazar PKS2155-304 [88l. A two order 
of magnitude flux increase, reaching 10 Crab Units (C.U., defined in ^6.11) was observed 
during a one hour period. Such flaring episodes are interesting targets of opportunity for 
neutrino telescopes. Assuming that half of the 7-rays are accompanied by the production 
of neutrinos, a flare of 10 C.U. lasting around 2.5 days would result in a neutrino detection 
at the significance level of 3 sigma 116 ill . 



3.2.2 Gamma ray bursts (GRBs) 

GRBs are short flashes of 7-rays, lasting typically from milliseconds to tens of seconds, 
and carrying most of their energy in photons of MeV scale. The likely origin of the GRBs 
with duration of tens of seconds is the collapse of massive stars to black holes. Observa- 
tions suggest that the formation of the central compact object is associated with Ib/c type 
supernovae 1 89, 90. 9 1 1 . 

GRBs also produce X-ray, optical and radio emission subsequent to the initial burst 
(the so called afterglow of the GRB). The detection of the afterglow is performed with 
sensitive instruments that detect photons at wavelengths smaller than MeV 7-rays. In 1997 
the Beppo-Sax 192] satellite obtained for the first time high-resolution X-ray images of the 
GRB970228 afterglow, followed by successive observations in optical and longer wave- 
lengths with an angular resolution of arcminute. This accurate angular resolution allowed 
the redshift measurement and the identification of the host galaxy. It was the first step to 
demonstrate the cosmological origin of GRBs. 
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Leading models assume that a fireball, produced in the collapse, expands with an 
highly relativistic velocity (Lorentz factor F ~ 10^'^) powered by radiation pressure. Pro- 
tons accelerated in the fireball internal shocks lose energy through photo-meson interaction 
with ambient photons (the same process of eq. In the observer frame, the condition 
required to the resonant production of the A+ is E^Ep = 0.2 GeV^F^. For the production 
of gamma-rays with ~ 1 MeV the characteristic proton energy required is Ep = 10^^ 
eV, if F ~ 10^^. The interaction rate between photons and protons is high due to the 
high density of ambient photons and yields a significant production of pions. The charged 
ones decay in neutrinos, typically carrying 5% of the proton energy. Hence, neutrinos with 
E,^ ~ 10^^ eV are expected [93 1. Other neutrinos with lower energies can also be produced 
in different regions or stages where GRB 7-rays are originated. Depending on models, a 
different contribution of neutrinos is expected at every time stage of the GRB. For instance, 
the neutrino emission from early afterglows of GRBs, due to dissipation made by the exter- 
nal shock with the surrounding medium or by the shock internal dissipation, was discussed 
in ||94] . Here, the implications of recent Swift [95 1 satellite observations concerning the 
possible neutrino signals in a neutrino telescope were also considered in detail. 

Some calculations of the neutrino flux |W| from GRB show that a kilometer-scale neu- 
trino telescope can be sufficient to allow detection. The average energy of these neutrinos 
(100 TeV) corresponds to a value for which neutrino telescopes are highly efficient. Nev- 
ertheless, being transient sources, GRBs detection has the advantage of being practically 
background free, since neutrino events are correlated both in time and direction with 7- 
rays. As for the case of the ANTARES detector |[97ll98l . unfiltered data can also be stored 
in the occurrence of a GRB alert from a satellite or a groundbased telescope (see Fig. [6] 
||99|). The analysis of collected data around a GRB alert can be carried out some time later, 
with the advantage of using very precise astronomical data, improved by later observations 
of the afterglow with optical telescopes. 



3.2.3 Starburst or neutrino factories 

Radio observations have motivated the idea of the existence of regions with an abnormally 
high rate of star formation, in the so-called starburst galaxies, which are common through- 
out the Universe. These regions of massive bursts of star-formation can dramatically alter 
the structure of the galaxy and can input large amounts of energy and mass into the inter- 
galactic medium. Supernovae explosions are expected to enrich the dense star forming re- 
gion with relativistic protons and electrons [100, 101 1 . These relativistic charged particles, 
injected into the starburst interstellar medium, would lose energy through pion production. 
Part of the proton energy would be converted into neutrinos by charged meson decays and 
part into 7-rays by neutral meson decays. Very recently, the (relatively) nearby NGC 253 
galaxy in the southern hemisphere and the M82 galaxy in the northern hemisphere were 
identified as starburst galaxies by, respectively, the H.E.S.S. [102| and VERITAS Iil031 
telescopes. The 7-ray flux above 220 GeV measured from NGC 253 imply a cosmic -ray 
density about three orders of magnitude larger than that in the center of our Galaxy. Such 
hidden accelerators of CRs are thus intense neutrino sources, since mainly neutrinos would 
be able to escape from these dense regions. A cumulative flux of GeV neutrinos from star- 
burst galaxies was calculated in |104| as E^<^>^ ~ lO"'^ GeV cm ^s ^sr ^ , a level which 
can be detected by a km'^-scale neutrino detector. 
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Figure 6: Sketch of an external alert from a GRB satellite to the ANTARES neutrino tele- 
scope. When an alert from the GRB Coordinates Network (GCN, which include active 
satellites) is received by ANTARES, all raw data in memory covering a few minutes are 
saved to disk. Any possible neutrino signal from the GRB (before, during, and after the 
photon detection by the satellite) is stored on disk Il99\l . 

3.3 The upper limits for transparent sources 

The observation of diffuse flux of gamma-rays and of UHE CRs can be used to set theoret- 
ical upper bounds on the total flux of neutrino from extragalactic sources (diffuse neutrino 
flux). High energy 7-rays can be produced in astrophysical acceleration sites by decay of 
the neutral pion (eq. |4l). Neutrinos will be produced in parallel from decay of the charged 
pions and they will escape from the source without further interactions, due to their low 
cross section. High-energy photons from ir" decay, on the contrary, will develop electro- 
magnetic cascades when interacting with the intergalactic radiation field. Most of the 7-ray 
energy will be released in the 1 MeV-100 GeV range. Therefore, the observable neutrino 
flux (within a factor of two due to the branching ratios and kinematics at production of 
charged and neutral pions) is limited by the bolometric observation of the gamma-ray flux 
in this energy band. 

The diffuse gamma-ray background spectrum above 30 MeV was measured by the 
EGRET experiment as [105]: 

E^I^{E) = (1.37 ± 0.06) X 10"^ (GeV cm^^sr^^s"^) (7) 

If nucleons escape from a cosmic source, a similar bound can be derived from the 
measured flux of CR from extragalactic origin. Fermi acceleration mechanism can take 
place when protons are magnetically confined near the source. Neutrons produced by photo- 
production interactions of protons with radiation fields (eq. [3]l can escape from transparent 
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sources and decay into cosmic protons outside the region of the magnetic field of the host 
accelerator. 

Some additional factors have to be considered before establishing a relationship be- 
tween CR and neutrino fluxes. These factors take into account the production kinematics, 
the opacity of the source to neutrons and the effect of propagation. This last factor is the 
subject to the larger uncertainties, because it has a strong dependence on galactic evolution 
and on the poorly-known magnetic fields in the Universe. There is some controversy about 
how to use relationships to constrain the neutrino flux limit. There are however two relevant 
predictions: 

• The Waxman-Bahcall upper bound. Following Mannheim fET), the upper bound 
proposed by Waxman-Bahcall 111061 (W&B) takes the cosmic-ray observations at 
EcR ~ 10^^ eV to constrain neutrino flux. With a simple inspection of Fig. [T] 
we can see that E'^dN/dE ~ 10"® GeV cm^^s^^sr^^ at 10^^ eV. This flux is two 
orders of magnitude lower than the limit provided by the extragalactic MeV-GeV 
7-ray background (eq. |7]l. 

In the computation of the upper bound, several hypothesis are made: it is assumed that 
neutrinos are produced by interaction of protons with ambient radiation or matter; that 
the sources are transparent to high energy neutrons; that the 10^^ eV CRs produced 
by neutron decay are not deflected by magnetic fields; finally (and most important) 
that the spectral shape of CRs up to the GZK cutoff is dN/dE oc E~'^, as typically 
expected from the Fermi mechanism. The upper limit that they obtain is: 



Although this limit may be surpassed by hidden or optically thick sources for protons 
to p7 or pp(n) interactions, it represents the "reference" threshold to be reached by 
large volume neutrino detectors (see Fig. |7]l. 

• Mannheim-Protheroe-Rachen (MPR) upper bound. The W&B limit was criti- 
cized as not completely model-independent. In particular, the main observation was 
about the choice of the spectral index a = 2. In 1 107 ] a new upper bound was de- 
rived using as a constraint not only the CRs observed on Earth, but also the observed 
gamma-ray diffuse flux. The two cases of sources opaque or transparent to neu- 
trons are considered; the intermediate case of source partially transparent to neutrons 
give intermediate limits. 

The limit for sources opaque to neutrons is: 



This is two orders of magnitude higher than the W&B limit and similar to the EGRET 
limit on diffuse gamma rays (eq. |7]), because a source opaque to neutrons produces 
very few CRs (neutrons cannot escape and cannot decay outside the source), but it is 
transparent to neutrinos and 7-rays. This limit was already excluded in a wide energy 
range by the AMANDA-II experiment, as shown in Fig. |7] 

The limit for sources transparent to neutrons decreases from the value of eq. |9] at 
Ey ~ 10*^ GeV to the value of eq. ^at E^, ~ 10^ GeV. Above this energy, the hmit 
increases again due to poor observational information. 



Eld^/dEy < 4.5 X 10"® (GeV cm-^sr-^s"^) 



(8) 



Eld^/dEy < 2 X 10"^ ( GeVcm^^sr^^s^^) 



(9) 



20 



T. Chiarusi and M. Spurio 



The W&B and the MPR Umits for neutrino of one flavor are reported in Fig. |7] The orig- 
inal values are divided by two, to take into account the neutrino oscillations from the source 
to the Earth (see ^4.3I ). Experimental upper limits are indicated as solid lines, ANTARES 
11108 II and IceCube 90% C.L. sensitivities for 3 years with dashed lines. Frejus |1091, 
MACRO fTTO], Amanda-II 2000-03 [1121 limits refer to muon neutrinos. Baikal [111] 
and Amanda-II UHE 2000-02 II113I refer to neutrinos of all-flavors. In this case, the orig- 
inal upper limits are divided by three. In fact, due to neutrino oscillations, we expect at 
Earth a flux of cosmic neutrinos of all flavors in the same proportion. The red line inside 
the shadowed band represents the Bartol [157] atmospheric neutrino flux. The lowest limit 
of the band represents the flux from the vertical direction, with a negligible contribution 
from prompt neutrinos. The upper limit of the band represents the flux from the horizontal 
direction, with one of the prompt model which gives the maximum contribution [ 1 14]. 

Extragalactic neutrino sources at the MPR-limit above 100 TeV cannot be experimen- 
tally excluded yet owing to the Earth occultation. If the main neutrino production mecha- 
nism is the photo-meson production due to the interaction of protons with radiation fields, 
the neutrino output becomes maximal at energies above 100 TeV. Due to the increasing 
neutrino interaction cross-section, the Earth becomes optically thick (see ^6.21 ) to neutri- 
nos above this energy. Ultra High Energy neutrino events must thus be looked for near 
the horizontal direction. As we will discuss in ^6.31 the experimental neutrino detection 
probability is expressed by the effective neutrino area. This quantity strongly depends from 
the background suppression capability of the detector for horizontal or downward-going 
events. 



4 Particle and fundamental physics with neutrino telescopes. 

Neutrino detectors can contribute to the multi-messenger astronomy, to solve some of the 
outstanding problems of high energy astrophysics described in the previous sections. In 
addition, these experiment will address some of the fundamental questions of high energy 
physics beyond the standard model: search for relic particles in the cosmic radiation; what 
is the nature of the Dark Matter; neutrino oscillations through the "standard" mass-flavors 
mechanism and with possible subdominant effects, as those induced by the violation of the 
Lorentz invariance or of the equivalence principle. 

4.1 Relic particles in the cosmic radiation 
4.1.1 Magnetic Monopoles 

Most of the Grand Unified Theories (GUTs) predict the creation of magnetic monopoles 
(MM) in the early Universe II115II116I . MM are topologically stable and carry a magnetic 
charge defined as a multiple integer of the Dirac charge §£, = hc/2e, where e is the elemen- 
tary electric charge, c the speed of light in vacuum and fi the Planck constant. Depending 
on the GUT group, the masses inferred for magnetic monopoles can take range over many 
orders of magnitude, from 10*^ to 10^^ GeV. 

MM are stable particles and they would have survived until now, diluted in the Universe, 
as predicted by theoretical studies which set limits on their fluxes, like the Parker flux limit 
|[Tl7[. Stringent limit s ove r a very wide MM velocity range were set by the underground 
MACRO experiment ifTTSl . 

Neutrino detectors can be used to search for fast MM. Since fast MM have a large 
interaction with matter, they can lose large amounts of energy in the terrestrial environment. 
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Figure 7: Sensitivities and upper limits for a E^'^ diffuse flux of high energy neutrinos 
of one flavor, see text. Experimental upper limits are indicated as solid lines, ANTARES 
and IceCube 90% C.E sensitivities with dashed lines. Upper limits obtained from all- 
flavor analyses (Baikal and Amanda-II UHE 2000-02) are not directly comparable to the 
Ufj^ upper limits. However, for the assumed astrophysical neutrino production models and 
for a wide range of neutrino oscillation parameters (' ^4.3^ , the flavor flux ratio at Earth can 
be assumed to be Vf, '■ i^fi '■ i^r = ^ ■ 1 : 1. In that case, either a single flavor limit 
can be multiplied by three (and compared to an all-flavor result) or an all-flavor limit can 
be divided by three and compared to a single-flavor result, as we did in the flgure. For 
reference, the W&B and MPR98 limits for transparent sources are also shown. Both upper 
bounds are divided by two, to take into account the neutrino oscillation effects. 



The total energy loss of a relativistic MM with one Dirac charge is of the order of 10 
GeV [11201 after having crossed the full Earth diameter. Because MM are expected to be 
accel erated in the galactic coherent magnetic field domain to energies of about 10^^ GeV 
II121I . some could be able to cross the Earth and produce an upgoing signals in a neutrino 
detector. 



The monopole magnetic charge g = mgjj can be expressed as an equivalent electric 
charge g = 68.5 me, where m is an integer. In a medium with a refractive index n, a 
MM with j3 = V / c> 1/n emits a factor (gn/ ze)"^ more Cherenkov light than an electrical 
charge ze with the same velocity. Thus relativistic MM with /3 > 0.74 carrying one Dirac 
charge will emit a large amount of direct Cherenkov light when traveling through the neu- 
trino detector; in water (n = 1.35) it gives rise to ~ 8500 more intense li ght th an a muon. 
Fig. [8] show the limits on MM from the Baikal |122| and AMANDA jllll Ch erenkov 
neutrino telescope, together with the best limits from other experiments II1181I124I . 



22 



T. Chiarusi and M. Spurio 



10 



B 

a 



10 



Parker 



Baikal (2003) 



Ohya (1991) 



MACRO (2002) 



AMANDA (2001) 



0.2 0.4 0.6 0.8 1 

Figure 8: Experimental upper limits (90% C.L.) on the MM flux as a function of the MM 
velocity f3c (solid lines). The Parker fTTT^ bound is also shown (dashed line). From IJ191 



4.1.2 Nuclearites 

Nuclearites are hypothetical nuggets of strange quark matter that could be present in cosmic 
radiation. Their origin is related to energetic astrophysical phenomena. Down-going nucle- 
arites could reach a neutrino telescope with velocities ~ 300 km/s, emitting blackbody ra- 
diation at visible wavelengths while traversing water/ice. Nuclearites with M > 10^^ GeV 
would be electrically neutral; the small positive electric charge of the quark core would be 
neutralized by electrons forming an electronic cloud. 

The relevant energy loss mechanism is represented by the elastic collisions with the 
atoms of the traversed media fl251. Nuclearites moving into the water/ice could be de- 
tected because of the black-body radiation emitted by the expanding cylindrical thermal 
shock wave. The luminous efficiency (defined as the fraction of dissipated energy appear- 
ing as light) was estimated, in the case of water, to be ~ 3 x 10^^ 1 125 1. The best limits on 
the search for nuclearites in the cosmic radiation are from the MACRO experiment 1 126 1. 
Because nuclearites are slowly moving particles, light coming from the passage of a nucle- 
arite on a km-scale neutrino detector may span into intervals from tens of /is up to 1 ms. 
Preliminary results were presented by the ANTARES collaboration II127II . 



4.2 Indirect dark matter searches 

The existence of nonbaryonic Dark Matter (DM) in our Universe is supported by strong 
cosmological observational evidences |128]. The presence of non- visible mass in galaxies 
is also motivated by the fact that dark matter halos seem to help stabilize spiral disk struc- 
ture. Nonbaryonic DM may consist of Weakly Interacting Massive Particles (WIMPs). A 
long list of nonbaryonic cold (i.e. non-relativistic) dark matter candidates has been sug- 
gested, among which the supersymmetric (SUSY) neutralino and the axion seem to be the 
most promising. SUSY postulates a symmetry between bosons and fermions predicting 
SUSY partners |J29il . In theories where R parity is conserved there exists a stable lightest 
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supersymmetric particle (LSP). If the neutralino is the LSP, it is a natural WIMP candidate: 
it is a weakly interacting particle with a mass between roughly a GeV and a TeV and would 
be expected to have a significant relic density. 

From the experimental point of view, direct and indirect methods H I 301 exist for de- 
tecting the WIMPs in the galactic halo. The methods can probe complementary regions of 
the supersymmetric parameter space, even when more extensive LHC results will become 
available. Direct methods detect Weakly Interacting Particles via the elastic scattering of the 
WIMP with a nucleus: the energy deposited in a low-background detector can be measured. 
Indirect methods look for by-products of WIMP decay or annihilation such as neutrinos re- 
sulting from the annihilation of WIMPs. Neutrino telescopes can perform indirect WIMP 
searches looking for high energy neutrinos from WIMP annihilation in the core of the Earth 
or the Sun. 

Dark matter WIMPs existing in the galactic halo can be captured in a celestial body 
by losing energy through elastic collisions and becoming gravitationally trapped. As the 
WIMP density increases in the core of the body, the WIMP annihilation rate increases 
until equilibrium is achieved between capture and annihilation. High energy neutrinos are 
produced via the hadronization and decay of the annihilation products (mostly fermion- 
antifermion pairs, weak and Higgs bosons) and may be detected as upward-going muons in 
neutrino telescopes. 

The capture rate for an astrophysical body depends, apart from the mass of the celestial 
body and from the escape velocity, on several poorly known factors: the WIMP mean halo 
velocity, the WIMP local density and the WIMP scattering cross section. The WIMP may 
scatter from nuclei with spin (hydrogen in the Sun) via an axial- vector spin-dependent inter- 
action. In this case, the WIMP couples to the spin of the nucleus or via a scalar interaction 
in which the WIMP couples to the nuclear mass. Elastic scattering is most efficient when 
the mass of the WIMP is similar to the mass of the scattered nucleus. Hence, the heavy 
nuclei in the Earth core make it very efficient in capturing WIMPs with mx < 100 GeV. 
Nuclei in the Sun, in contrast, have a smaller average mass; the Sun is nevertheless efficient 
in capturing WIMPs due to the larger value of the escape velocity. 

As studied by underground experiment like Super- Kamiokande fl311, MACRO II132I 
and others, WIMP annihilation signal would appear in a neutrino telescope as a statistically 
significant excess of upward-going muon events from the direction of the Sun or of the 
Earth among the background of atmospheric neutrino-induced upward-going muons. The 
precise direction measurement allows a restriction of the search for WIMP annihilation 
neutrinos to a narrow cone pointing from the Earth center or from the Sun, greatly reducing 
the background. This detection method achieves an increasingly better signal to noise ratio 
for high WIMP masses, because of the increase in neutrino cross section with energy and 
longer range of high energy muons. 

At present AMANDA 1 133ll an d Baikal pi 341 have published results on indirect search 
of WIMPs, while ANTARES 111351 has presented preliminary results. Instead of the general 
supersymmetry scenario, in ANTARES the more constrained approach of minimal super- 
gravity (mSUGRA) was used. mSUGRA models are characterized by four free parameters 
and a sign: mi/2, rrio, Aq, tan(/3) and sgn(^). They investigated mSUGRA models that 
possess a relic neutralino density that is compatible with the cold dark matter density Qx 
as measured by the WMAP experiment. No excess is found, and the sensitivity has been 
sufficient to put constraints on parts of the mSUGRA parameter space. 

Both the IceCube and a cubic kilometer experiment in the Mediterranean sea would be 
sensitive to a wide part of the mSUGRA parameter phase-space. 
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4.3 Neutrino oscillations 

In recent years, neutrino oscillation became a well known phenomenon, which plays also 
an important role on determining the flavor on Earth of neutrinos of cosmic origin. Neu- 
trino oscillations were observed in atmospheric neutrinos, in solar neutrino experiments 
and on Earth based accelerator and reactor experiments. A complete review about neutrino 
oscillations can be found in 1 136|. 

As already mentioned, high energy neutrinos are produced in astrophysical sources 
mainly through the decay of charged pions, in p^, pp, pn interactions (eq. [5]l. Therefore, 
neutrino fluxes of different flavors are expected to be at the source in the ratio: 

Ue:v^,:ur = l : 2 : Q (10) 

Neutrino oscillations will induce flavor changes while neutrinos propagate through the 
Universe. One has to consider mass eigenstates Um = vi, U2, fs in the propagation, 
instead of weak flavor eigenstates ui = Ue, z/^, Vt- The weak flavor eigenstates ui are 
linear combinations of the mass eigenstates Um through the elements of the mixing matrix 

Ulm- 

3 

'^1= ^ Ulm (11) 
m=l 

Because mixing angles are large, the flavor eigenstates are well separated from those of 
mass. The oscillation probability in the simple case of only two flavors, for instance 
(f^, Vt) and one mixing angle 023> is: 

f f ^ = z^2Cos 6*23 + i^3sin 6I23 ,^2) 
I = -1/2 sin 023 + cos 023 

The survival probability for a pure a beam: 

o „ / 1.27Am2 • l\ 
P{yi, ^v^,) = l- sin' 2023 sm' I ^ J (13) 

where Am' = m\ — m' (eV'), L (km) is the distance travelled by the neutrino from 
production to detection and Ey (GeV) the neutrino energy. 023 and Am^ may be experi- 
mentally determined from the variation of P(t'^ ^ 1/^) as a function of the zenith angle or 
from the variation mL/ Ey. 

With three neutrino flavors, three mass differences can be defined (two linearly inde- 
pendent). The mass difference measured with atmospheric neutrinos is Am^^^ ~ Am23 
= ±2.5 X 10-3 eV'. For the mixing angle, 023 ~ 45° (that correspond to maximal mix- 
ing), while 013 is small. The values of Amfg and of the other mixing angle 0i2 are de- 
termined by the solar neutrino experiments and KamLand. The most recent data favor 
II137II very clearly the solution with a best fit: Am^^^ ~ Am'2 = 7.6 x lQ~^eV^ and 
sin' 012 = 0.32 (012 = 34.4°). 

According to these neutrino oscillation parameters, the ratio of fluxes of neutrinos from 
astrophysical origin (i.e. very large baseline L) in eq. [TO] changes to a flux ratio at Earth 
IIT38I as: 

Ve:Va:yT = l ■■ 1 : 1 (14) 
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Theoretical predictions which does not take into account neutrino oscillations must be 
corrected to include this effect (as we did in Fig. |7] ). In particular, muon neutrinos are 
reduced at Earth by a factor of two. 

Neutrino telescopes can detect thousand of atmospheric per year. For instance, the 
ANTARES detector is efficient for > few tens of GeV, were atmospheric upward-going 
are still suppressed by flavor oscillations. Neutrino telescopes can test also non-standard 
oscillations. If the standard mass-induced oscillation is assumed as the leading process 
for flavor changes, other mechanisms can be tested for flavor transitions as a subdominant 
effect. As an example, we discuss in the following the case of the Lorentz invariance 1, 139 J . 

4.4 Violation of the Lorentz invariance 

Quantum gravity theories assume that the spacetime take a foamy nature 111401 . Interactions 
with this spacetime foam may lead to the breaking of CPT symmetry, leading to the viola- 
tion of Lorentz invariance Ill41i In addition, some theories of quantum gravity predict that 
there is a minimum length scale, of order the Planck length (10^'^^ m): the existence of a 
fundamental length scale may also induce the violation of Lorentz invariance (VLI). 

Lorentz invariance violation may manifest itself in many different ways and can be 
tested with many different experimental systems, in particular with a modified neutrino 
oscillation lengtlQ. 

The VLI subdominant effect can be studied using atmospheric neutrino data collected 
by neutrino telescopes, after correction for the known mass-induced neutrino oscillation. 
In this scenario [142], neutrinos can be described in terms of three distinct bases: flavor 
eigenstates, mass eigenstates and velocity eigenstates, the latter being characterized by dif- 
ferent maximum attainable velocities (MAVs) in the limit of infinite momentum. Here, both 
mass-induced oscillations and VLI transitions are treated in the two-family approximation. 
It is also assumed that mass and velocity mixings occur inside the same families (e.g., 1/2 
and u^). In the case of mass-flavor oscillation, the survival probability of muon neutrinos at 
a distance L from production is given by eq. [13] In the VLI case, the z/^ survival probability 
is: 

P{j^t, ^i^f,) = l- sin^ 2(9^ sin^ (2.54 x lO^^AvLE^^ (15) 

where Av = {v^.^ — Vy^) is the neutrino MAV difference in units of c and 9^ is the mixing 
angle. Notice that neutrino flavor oscillations induced by VLI are characterized by an LE^ 
dependence of the oscillation probability (eq. [T5l ). to be compared with the L/Ey behavior 
of mass-induced oscillations (eq. [T3l ). 

When both mass-induced transitions and VLI induced transitions are considered simul- 
taneously, the muon neutrino survival probability can be expressed as 

P{v^, ^v^) = l- sin^ 26 sin^ (16) 



where 20 = atan{ai/a2) and VL = ^0^ + 02, with oi and 02 which depend from 

Am^, ^23, Aw, 9^, L and E^ lfT42l. 

The same formalism also applies to violation of the equivalence principle, after substi- 
tuting Av/2 with the adimensional product | cj) \ A7 ; A7 is the difference of the coupling 
constants for neutrinos of different types to the gravitational potential 1 144]. 

' In the literature, neutrino oscillations induced by the simplest models of VLI and violation of the equiva- 
lence principle (VEP) are described within the same formalism. In the following we will mention only this VLI 
formalism for simplicity. We refers to 1 140] for the more general case. 
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The most conse rvative bounds from underground experiments were obtained by Super- 
Kamiokande [139] and MACRO II 1421 . In particular, the 90% confidence level Umits ob- 
tained by MACRO are \Av\ < 6 x IQ-^'^ at sin2(9„ = and \Av\ < 2.5 - 5 x lO"^'^ at 
sm29y = ±1. The typical energy scale of upward thoroughgoing neutrino-induced muons 
measured from these experiments is of the order of E'iy ~ 10 ^ 100 GeV. Neutrino tele- 
scopes are sensitive to much higher neutrino energies. The larger number of atmospheric 
neutrino events and the greater average energy allows neutrino telescopes to be much more 
sensible. The AMANDA-II expected sensitivity (90% C.L.) for maximal mixing and six 
years of simulated data is \Av\ < 2.1 x 10^^^ ||143|| . A null observation would be able to 
place very stringent bounds on quantum decoherence effects, and on VLI parameters which 
modify the dispersion relation for massive neutrinos. 

4.5 HE neutrinos in coincidence with gravitational waves 

High-energy neutrinos and gravitational waves (GW), contrarily to high-energy photons 
and charged cosmic rays may escape from dense astrophysical regions and travel over large 
distances without being absorbed, pointing back to their emitter. 

It is expected that many astrophysical sources produce both gravitational waves and HE 
neutrinos II 145I . A possible coincident detection will then provide important information 
on the processes at work in the astrophysical accelerators. Furthermore, if a mechanism (as 
for instance the gamma-ray bursts) allows a precise measurement of the time delay between 
neutrinos and GW signals, some quantum-gravity effects can be tested with the possibility 
to constrain some dark energy models 1 146 1. 

The search for coincident detection is motivated by the advent, in association of neutrino 
telescopes, of a new generation of GW detectors VIRGO 111471 and LIGO II148II (which are 
now part of the same experimental collaboration). The neutrino/GW association requires 
the measurement of the event time, arrival direction and associated angular uncertainties. 
Each candidate is obtained by the combination of reconstruction algorithms specific to each 
experiment, and quality cuts used to optimize the signal-to-background ratio. A joint neu- 
trino/GW analysis program is planned for the ANTARES neutrino telescope and VIRGO-i- 
LIGO [ 145 ]. The coincident (non-)observation shall play a critical role in our understanding 
of the most energetic sources of cosmic radiation and in constraining existing models. They 
could also reveal new, hidden sources unobserved so far by conventional photon astronomy. 



5 Neutrino detection principle 

The basic idea for a neutrino telescope is to build a matrix of light detectors inside a trans- 
parent medium. This medium, such as deep ice or water: 

• offers large volume of free target for neutrino interactions; 

• provide shielding against secondary particles produced by CRs; 

• allows transmission of Cherenkov photons emitted by relativistic particles produced 
by the neutrino interaction. 

Other possibilities, such as detecting acoustic or radio signals generated by EeV (10^^ 
eV) neutrinos in a huge volume of water or ice are not considered in this review |149'|. 
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High energy neutrino interact with a nucleon N of the nucleus, via either charged cur- 
rent (CC) weak interactions 

i^i + N^l + X (17) 
or neutral current (NC) weak interactions 

ui + N^ui + X. (18) 

At energies of interest for neutrino astronomy, the leading order differential cross sec- 
tion for the uiN ^ IX CC interactions is given by II150II 



dxdy TT (Q2 + M^Y 

X [xq{x,Q^)+x{l-yfq{x,Q'')] (19) 

where x = Q'^ jlmj^iEy — Ei) and y = {By — Ei)/Ey are the so-called scale variables 
or Fenyman-Bjorken variables, is the square of the momentum transferred between the 
neutrino and the lepton, mj\[ is the nucleon mass, Mw is the mass of the W boson, and 
Gp is the Fermi coupling constant. The functions q{x,Q'^) and q{x,Q'^) are the parton 
distributions for quarks and antiquarks. Fig. |9] shows the and cross sections as a 
function of the neutrino energy. As can be seen, at low energies the neutrino cross section 
rises linearly with Ei, up to ~ 10^ GeV. For higher energies, the invariant mass = 
2mNEyxy could be larger than the W-boson rest mass, reducing the increase of the total 
cross section. Since there is not data which constrain the structure functions at very small x, 
outside the range measured with high precision at the HERA collider, some uncertainties are 
estimated on the total cross section at large energies [ 151 1. Computer libraries [ 152] provide 
a collection of parton distribution function (PDF) to model the neutrino cross section also 
at very high energies. 

Cosmic neutrino detectors are not background free. Showers induced by interactions 
of cosmic rays with the Earth's atmosphere produce the so-called atmospheric muons and 
atmospheric neutrinos. Atmospheric muons can penetrate the atmosphere and up to several 
kilometers of ice/water. Neutrino detectors must be located deeply under a large amount of 
shielding in order to reduce the background. The flux of down-going atmospheric muons 
exceeds the flux induced by atmospheric neutrino interactions by many orders of magnitude, 
decreasing with increasing detector depth, as is shown in Fig. (TO] The previous generation 
of experiments which had looked also for astrophysical neutrinos (MACRO [153|, Super- 
Kamiokande 1154|) was located under mountains, and has reached almost the maximum 
possible size for underground detectors. 

Charged particles travel through the medium until they either decay or interact. The 
mean length of the distance travelled is called the path length of the particle and it depends 
on its energy loss in the medium. If the path length exceeds the spatial resolution of the 
detector, so that the trajectory of the particle can be resolved, one have a track. In a high 
energy neutrino detector, one can distinguish between two main event classes: events with 
a track, and events without a track (showers). 

Relativistic charged particles emit Cherenkov radiation in the transparent medium. A 
detector must measure with high precision the number and arrival time of these photons 
on a three-dimensional array of Photo Multiplier Tubes (PMTs), from which some of the 
properties of the neutrino (flavor, direction, energy) can be inferred. 
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Figure 9: Cross section for andV^ as a function of the (anti)-neutrino energy according 
to CTEQ6-DIS USI]! parton distributions. 

In order to behave as a neutrino telescope, a neutrino detector must be able to point at 
a specific celestial region if a signal excess over the background is found. Neutrino tele- 
scopes must have the same peculiarities of GeV-TeV 7-ray experiments (satellites, imaging 
Cherenkov) to associate some of the signal excesses to objects known in other electromag- 
netic bands. In order to achieve an angular resolution of a fraction of degree, only the CC 
f^j interaction can be used. The angular resolution for other flavors and for NC is so poor 
that there is no possibility to perform associations. For the same reason, the particle physics 
and general physics open questions which can be afforded with a neutrino telescope largely 
rely on the i/^j channel. 

On the other hand, a high energy neutrino detector is motivated by discovery and must 
be designed to detect neutrinos of all flavors over a wide energy range and with the best 
energy resolution. This is of particular interest for the case of the neutrino diffuse flux from 
extragalactic sources. In addition, the neutrino oscillation changes the source admixture 
from z^e : : = 1 : 2 : to 1 : 1 : 1. While above hundreds of TeV muon and electron 
neutrinos become absorbed by the Earth, the tau neutrino is regenerated \ 158 J : high energy 
Ur will produce a secondary Ur of lower energy, lowering its energy down to 10^^ eV, where 
the Earth is transparent. 

Schematic views of a z^g , i^fi and Ur CC events and of a NC event are shown in Fig. [TT] 
Neutrino and anti-neutrino reactions are not distinguishable; thus, no separation between 
particles and anti-particles can be made. Showers occur in all event categories shown in 
Figure. However, for CC f^, often only the muon track is detected, as the path length of a 
muon in water exceeds that of a shower by more than 3 orders of magnitude for energies 
above 2 TeV. Therefore, such an event might very well be detected even if the interaction has 
taken place several km outside the instrumented volume, provided that the muon traverses 
the detector. 

Neutrino telescopes, at the contrary of usual optical telescopes, are 'looking downward'. 
Up-going muons can only be produced by interactions of (up-going) neutrinos. From the 
bottom hemisphere, the neutrino signal is almost background-free. Only atmospheric neu- 
trinos that have traversed the Earth, represent the irreducible background for the study of 
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Figure 10: Different contributions (as a function of the cosine of the zenith angle) of the: i) 
atmospheric muons (computed according to l[156\l } for two different depths; ii) atmospheric 
neutrino induced muons (from I d 57] }. for two different muon energy thresholds. 

cosmic neutrinos. The rejection of this background depends upon the pointing capability of 
the telescope and its possibility to estimate the parent neutrino energy. As we will discuss 
in f|7J either water or ice is used as media. A deep sea-water telescope has some advantages 
over ice and lake-water experiments due to the better optical properties of the medium. 
However, serious technological challenges must be overcome to deploy and operate a de- 
tector in deep sea, as we will discuss in ^ 

5.1 Electron neutrino detection 

A high energy electron resulting from a charged current interaction has a high proba- 
bility to radiate a photon via bremsstrahlung after few tens of cm of water/ice (the water 
radiation lenght is ~ 36 cm). The following process of e+e^ pair productions, and sub- 
sequent bremsstrahlung, rapidly produce an electromagnetic (EM) shower until the energy 
of the constituents falls below the critical energy Ec and the shower production stops; the 
remaining energy is then dissipated by ionization and excitation. 

The Cherenkov light from EM showers is emitted isotropically in azimuth with re- 
spect to the shower axis. The lateral extension of a EM shower is of the order of 10 cm 
P1371 and therefore negligible compared to the longitudinal one. Thus, the EM shower 
is described by the longitudinal shower profile (a parametric formula is given in II137II ). 
The longitudinal profiles are used to parameterize the total shower length L as a function 
of the initial shower energy. The shower length is defined as the distance within which 



95% of the total shower energy has been deposited. For salt water II159II it is found that 
L[m] = 3.04 + 1.091ogio(-E^[GeV]). For a 10 TeV electron L is 7.4 m 
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Figure 1 1 : Some event signature topologies for different neutrino flavors and interactions: 
a) CC interaction of a produces a muon and a hadronic shower; b) CC interaction 
of a u-T produces a r that decays into a Vr, tracing the double bang event signature, c) 
CC interaction of Uf. produces both an EM and a hadronic shower; d) a NC interaction 
produces a hadronic shower Particles and anti-particles cannot be distinguish in neutrino 
telescopes. From U59\l 

A showers size of order of 10 m is small compared to the spacing of the PMTs in any 
existing or proposed detector. EM showers represent, to a good approximation, a point 
source of Cherenkov photons. Pointing accuracy for showers is inferior to that can be 
achieved for the channel. Reconstruction of Monte Carlo simulated events performed in 
the framework of the IceCube and ANTARES collaborations shows a precision of the order 
of ~ 10°, with the possibility to reduce it to few degrees for a small subsample of events. 

Finally, we should mention (only for completeness) two effects which has some conse- 
quences on EM shower: 

- the Glashow resonance 1116011 . which affects the Fg through the resonant process Ve+e~ — > 
W~ — > q + q, {ui + i). The resonance peak is for a I7g energy of 6.3 PeV. This resonant 
channel constitutes only a small portion to the overall cross section in the energy range 
between 100 GeV and 100 PeV and it must be taken into account in Mont Carlo simulations; 

- the Landau- Pomeranchuk - Migdal (LPM) effect fl 611ll62ll . which has an influence on 
showers of ultra high energies {E > 10^^ eV). The LPM effect suppresses the radiative 
energy losses of the particles in the shower; in this case, the longitudinal development of 
EM or hadronic cascades can be largely enhanced. 

5,1,1 Neutral currents interactions 

The NC channel gives the same signature for all neutrino flavors. In this channel, a part 
of the interaction energy is always carried away unobserved by the outgoing neutrino, and 
therefore the error on the reconstructed energy of the primary neutrino increases accord- 
ingly. Even though EM and hadronic showers are different from each other in principle, 
the Ve CC and the NC channels are not distinguishable in reality, because any proposed 
detector is too sparsely instrumented. 

Hadronic cascades suffer event-to-event fluctuations which are much more important 
with respect to the EM ones |159I . The dominant secondary particles in a hadronic shower 
are pions; kaons, protons or neutrons occur in variable fractions. Muons (from pions de- 
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Figure 12: Path length of particles produced by neutrino interactions in water: muons, taus, 
electromagnetic and hadronic showers, versus their respective energy. The shower lengths 
are calculated using a shower profile parameterization as described in H159\l . 



cay) can be present as well: as these usually leave the shower producing long tracks, they 
contribute significantly to the fluctuations. 

Monte Carlo simulations (for instance, from the ANTARES collaboration ||163|) shows 
that above 1 TeV of shower energy, the largest part of the Cherenkov light is generated 
by EM sub-showers. Referring to Fig. [T2j the EM shower has a shorter length than the 
hadronic shower below 100 TeV. The electromagnetic component (produced by the vr*^ de- 
cay) in the hadronic shower increases with increasing shower energy. The longitudinal 
profile of hadronic showers can be parameterized in the same way as the EM one. The 
shower length Lh is again defined as the distance within which 95% of the total shower en- 
ergy is deposited. It can be described as: LhIpo] = 5.28 + 0.701ogxo(^[GeV]). For what 
concerns the measurement of the incoming neutrino direction, the angular difference be- 
tween the shower and the neutrino falls below 2° for E,^ above ~ 1 TeV. It is thus negligible 
with respect to the precision of the shower direction measurement. 

5.2 Tau neutrino detection 

For Ur CC interactions, the produced r-lepton travels some distance (depending on its en- 
ergy) before it decays and produces a second shower. The Cherenkov light emitted by the 
charged particles in the showers can be detected if both the Ur interaction and the r decay 
occur inside the instrumented volume of the detector. Below 1 PeV, also the Ur CC channels 
(except for the case where the r produces a muon) belong to the class of events without a 
track, because the r track cannot be resolved. 

The r lepton has a short lifetime, and in the energy range of interest it travels (depending 
on the Lorentz factor Er/mrC^) from a few meters up to a few kilometers before it decays 
(see Fig. [T2l ). Radiation losses, on the other hand, play a much smaller role than for the 
muon, because of its 17 times larger mass. Most of the possible r decay modes include the 
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generation of a hadronic or an electromagnetic cascade. Thus, if the track of the r is long 
enough to distinguish between the primary interaction of the Ur and the decay of the tau 
(typically for r energies above 1 PeV, see Fig. [121) . the expected signatures for the Ur CC 
events are that of a shower, plus a track, plus another shower. This signature is called double 
bang event, if the Cherenkov light emitted by the charged particles in the first shower can 
be detected and separated from the light emitted by the particles produced in the r decay. 

Alternatively, if the r starts or ends outside the instrumented volume, a track plus a 
shower can be detected. This signature is called lollipop event. In small size neutrino 
detector (like ANTARES) the expected lollipop event rate above 1 PeV is far below 1 event 
per year. For larger detector (IceCube, the 1 km'^ Mediterranean sea telescope) the optimal 
Ur energy value for double bang events is around 10^^ eV, because the tau path length 
rapidly exceeds the dimensions of the detectors for increasing energies. If the r decays into 
a muon, the event is presumably not distinguishable from an original CC interaction. 



5.3 Muon neutrino detection 




(GeV) 

Figure 13: Effective muon range as a function of the initial energy E^^. Curves correspond 
(from top to bottom) to different threshold energies -E^^^ of the muon arriving at the detector 
(E>^f^^. = 1, 102, 10^, 10"^, 10^ lO'^ GeV). From / f765l/ . 

Muon neutrinos are especially interesting in a search for cosmic point sources of neutri- 
nos with energies larger than ~ 1 TeV. In this energy range, interaction can occur outside 
the detector volume, while in most cases muons are energetic enough to completely traverse 
the detector. This gives a clean experimental signal which allows accurate reconstruction 
of muon direction, closely correlated with the neutrino direction. 

The relation between neutrino and muon directions is essential for the concept of a 
neutrino telescope. Since neutrinos are not deflected by (extra-) galactic magnetic fields, it 
is possible to trace the muon back to the neutrino source. This is equivalent to traditional 
astronomy where photons point back to their source. The average angle 9^^ between the 
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incident neutrino and the outgoing muon can be approximated by: 

0.6° 



'^^ < -^=== (20) 



where is the neutrino energy (see Fig. | 

Muon energy losses are due to several processes such as ionization, pair production, 
bremsstrahlung and photonuclear interactions 111641 . The total energy loss per unit length 
can be written in a parameterized formula as: 

dE^/dx = a{E^) + P{E^) ■ E^ (21) 

where a{E^) is an almost constant term that accounts for ionization, and (3{E^) takes into 
account the radiative losses. 

Fig. [13] shows the effective muon range Reff in water. Reff represents the range after 
which a muon of initial energy E^ has still a residual energy E^^^^ at the detector. As an 
example, a muon with initial energy E^ = 10 TeV travels more than 4 km in water and 
arrives with more than 1 TeV of residual energy. The event will be detected even if the 
neutrino interaction vertex is outside the instrumented detector volume of the telescope. 

A muon travelling through rock, ice or water is subject to multiple scattering. The 
devia tion of the muon direction due to this process after travelhng a distance x is given by 

imi: 

^ 13.6{MeV) 



x/Xq[1 + 0.0038 ln(x/Xo)] (22) 



where Xq is the radiation length of the medium. At energies and distances that concern us, 
6ms is smaller than 9u^ and the effect can be neglected. Muon direction can be measured 
with an intrinsic resolution which depends on many factors, and in particular on the propa- 
gation medium. From Monte Carlo simulations, the precision is of the order of less than 1° 
in ice, and ~ 0.2° in water (see 



5.4 Cherenkov radiation 

Any operating or proposed neutrino telescope in the TeV-PeV range is working by collect- 
ing the optical photons produced by the Cherenkov effect of relativistic particles. The light 
is collected by a three-dimensional anay of photomultiplier tubes (PMTs). The information 
provided by the number of photons detected and their arrival times are used to infer the 
neutrino flavor, direction and energy. 

Cherenkov radiation is emitted by charged particles crossing an insulator medium with 
speed exceeding that of light in the medium II166I . The charged particle polarizes the 
molecules along the particle trajectory, but only when the particle moves faster than the 
speed of light in the medium, an overall dipole moment is present. Light is emitted when the 
electrons of the insulator restore themselves to equilibrium after the disruption has passed, 
creating a coherent radiation emitted in a cone with a characteristic angle Oq given by 

cosOc = ^ = 7^ (23) 

where n is the refracting index of the medium and (3 is the particle speed in units of c. For 
relativistic particles (/3 ~ 1) in seawater (n ~ 1.364) the Cherenkov angle is 9c ~ 43°. 
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Figure 14: (a) Sketch of an ANTARES optical module ( OM). Those used by IceCube, NEMO 
and NESTOR experiments are similar A large hemispherical (10 inches in diameter) pho- 
tomultiplier (PMT) is protected by a pressure-resistant glass sphere. The outer diameter of 
the sphere is 43.2 cm. A mu-metal cage protects the PMT from the Earth magnetic field. An 
internal LED is used for the calibration, (b) The quantum efficiency for PMTs commonly 
used in ice or water (from Hamamatsu ). 



The number of Cherenkov photons, N^, emitted per unit wavelength interval, dX and 
unit distance travelled, dx, by a charged particle of charge e is given by 

^ = ^(1 - ^) (24) 
dxdX 137A2^ n^(3^' ^ ' 

where A is the wavelength of the radiation. From this formula it can be seen that shorter 
wavelengths contribute more significantly to the Cherenkov radiation. The light absorption 
by water/ice will strongly suppress some wavelengths, see ^ 

Fig. [l4]shows one example of an optical module used in ice and water experiments (see 
f|9]and ^ flOl ). The PMT quantum efficiency is large in the wavelength range between 300- 
600 nm, matching very well the region in which ice/water are transparent to light. Typically, 
in the wavelength range between 300-600 nm, the number of Cherenkov photons emitted 
per meter is about 3.5 x 10^. 



6 Why a km^ telescope 

In this section we develop a toy model for a neutrino telescope, dedicate to non-expert 
readers. Our aim is to derive, using a simple calculation, why a cubic kilometer scale 
detector is needed, and what is the number of optical sensors required to detect neutrinos 
in the instrumented volume. In ^2.41 and ^3.21 we have presented some theoretical models 
of, respectively, galactic and extragalactic neutrino sources. Each model is characterized 
by the differential neutrino energy spectrum (TeV~^ cm^^s^^). The rate of detectable 
events is simply derived using the neutrino effective area. This quantity (i ]6.3l) is computed 
by each experiment, and can be used to compare the sensitivity of the different experiments. 
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6.1 The reference flux 

As an ingredient in our toy model, we need a reference flux. TeV 7-rays astronomers found 
convenient to use the Crab flux as reference flux or standard candle for their measurements. 
For practical reasons, this is useful also for neutrino astronomy. 

The Crab nebulifl was discovered at TeV energies in 1989 1 169| and it is conventionally 
used as a reference source of TeV 7-rays, due to its relative stability and high intensity. A 
reference flux equal to 1 C.U. (Crab Unit) is defined as the flux similar to that of the Crab, 
assuming no energy cutoff and a spectral index F = 2 : 

E^^ = 10"^^ cm-^s-^TeV = 1 C.U. (25) 

The recent Crab TeV 7-rays spectrum measured by different experiments shows a steepen- 
ing; a better fit of the data is provided by a power law (energy E in TeV) with an exponential 

cutoff: ^ = IoE-'^e^-^/^'=\ This formula, with T = 2.39±0.03stat and a cutoff energy 
£"2 = (14.3 ± 2.1 Stat) TeV, gives a differential flux normalization at 1 TeV of II170II : 
lo = (3.76 ± 0.07stat) X 10"" cm-2s-iTeV~i. 

Regarding neutrino telescopes, the Crab is in a particular sky position (see Fig. 
because it can be seen by both telescopes located in the Northern and in the Southern Earth 
hemisphere. In order to perform a back-of-the-envelope calculation, we assume a reference 

neutrino flux equal to that of the gamma flux from the Crab, E^^^ ~ E^^^^. This is 

only true if hadronic processes are involved in the Crab TeV 7-ray emission (and for this 
particular source it seems to be disfavored with respect to the leptonic model). A detailed 
prediction of the event rate in a neutrino telescope as big as the IceCube experiment is 
N{E^ > 1 TeV) 2.8 yr^^ 1,40,1 . Our simple calculation in the next section can be 
compared with this result. 



6.2 Event rate for the reference flux 

The CC interactions give the possibility of a measurement of the neutrino direction 
within 1° or better (it can reach up to ~ 0.2° in water). The channel allows also an 
enhancement of the effective volume of the detector. The effective volume corresponds 
to the product of the detector effective area and the muon range R^, it is larger than the 
instrumented volume and it increases with increasing neutrino energy. 
The event rate in a neutrino telescope can be expressed in terms of: 

= ldE..^^{E.).A.P.,iE.,Ei^,^.) 

. ^-a{E,.)pNAZ{8) /26) 



^The Crab pulsar, with a rotational period of 33 ms and a spin-down luminosity L — 5x 10^* erg/second, is 
a particular important source for high energy astrophysics. In fact, the pulsar powers a surrounding synchrotron 
nebula which has been detected in radio and X-ray wavelengths fT67 1. It is believed that the rotational energy of 
the pulsar is mostly carried away by a relativistic wind of electrons and positrons. The interaction of this wind 
with the surrounding medium creates a relativistic shock wave, where the leptons are thought to be accelerated 
to high energies [1681 . 
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where A and T are the detector area and observation time, respectively. P^^ is the proba- 
bility that the neutrino gives an observable muon (see below). The effective neutrino area 
is the quantity: 

Alff{E,) = P,^{E,, E'i^^) ■ A ■ e-(^^)P^A^(^?) (27) 

The term e^'^^^'''>P^^^'^^\ where a{E^) is the total neutrino cross section, Na the Avogadro 
number, (pNa) the target nucleon density and 9 the neutrino direction with respect to the 
nadir, takes into account the absorption of neutrinos along the Earth path Z{6). From the 
nadir Z{0) = 6.4 x 10^ cm, the absorption becomes not negligible for a > 10~^^ cm^ (or 
equivalendy E,, > 50 TeV, see Fig. when e'^f^^^ ~ e"° °^ ~ 0.95. 

PufiiEu, E^f^^) represents the probability that a neutrino with energy E^ produces a 
muon of energy E^ which survives with energy > E"^^ after the propagation from the 
interaction point to the detector. It can be expressed in terms of: 

X Reff{E^,E';^^)dE^ (28) 

where day/dE^j, is the differential neutrino cross section to produce a muon of energy E"^, 
and R(.ff{E^, ^t/ir) effective muon range. One can tabulate P^^ for a given muon en- 
ergy threshold. Fig. [15] shows P^^ for two values of muon threshold energy: 1 GeV (which 
was the characteristic muon threshold for large area underground detector, like MACRO) 
and 1 TeV (which is the characteristic value for neutrino telescopes). 




1 2 3 4 5 6 7 

Logio E„, GeV 



Figure 15: P^^ HI 7 IV for two values of the muon threshold energy: 1 GeV and 1 TeV. The 
solid line are neutrinos and the dashed line antineutrinos. The dotted lines show a power 
law approximation. 
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Note the dependence of Pj^^ for > 1 GeV on neutrino energy lfT2l : it can be 
approximated with Py^ oc El for Ey < I TeV, and Pj,^ oc Ey for 1 < E",^ < 10^ TeV. 
The two energy regimes reflect the energy dependence oi^ the neutrino cross section (which 
depends almost linearly on energy) and effective muon range (which depends linearly on 
muon energy up to ~ 1 TeV, when muon radiative losses become dominant). 

To solve analytically eq. [26] in our simplified model, we use the approximation for 
neutrino energies larger than 1 TeV: Py^{Ey,E^^^) ~ Po-E^'® = 10"^^°® {E in TeV) 
II171I . With this approximation, the event rate for a neutrino source equivalent to 1 C.U. 
(eq. |25] ) can be analytically computed in the range between 1 10'^ TeV, neglecting the 
Earth's absorption: 



rlO^ TeV 

/ dEy ■ (KE-^) ■ A ■ (PoE°-«^ 

Jl TeV 

5 X 10"^^ • A cmr'^s'^ (29) 



where K = 10"^^ TeV"^ cm s . Here, the area A is the surface surrounding the 
instrumented volume where muons of energy E^^^ can be detected. To give a rough energy- 
independent estimation, a 1 TeV muon has a range (in water) of RixeV ~ ^-^^ 1 164]. 
The cross sectional area of a sphere of diameter rp^y is A ~ 5 km^. Assuming this value 
in eq. |29l the number of expected events is ~ 1/year, in a rough agreement with the more 
detailed computation of [40 1. 



6.3 The neutrino effective area 

The energy-dependent effective area Alff{Ey) must be computed using Monte Carlo (MC) 
simulations. Because it depends from many factors, including the detector characteristics, 
Al,f f is strongly detector-dependent. The effective area allows a more detailed computation 
of the rate of expected events for a given neutrino flux and the sensitivity of each experiment 
can be expressed in terms of the effective area |172|. In the MC calculation, a flux of 
neutrino is generated, which interacts in a huge volume Vgen of material surrounding the 
instrumented volume. The effective neutrino area is then computed as: 

Al^f = -pL- X F,en X (pNa) X a{Ey) x e-(^'')''^-^(^) (30) 

-^^ gen 

where Nx is the number of detected (triggered, reconstructed) events and Ngen is the num- 
ber of generated events. The other quantities have been already defined. Fig. [16] shows the 
effective neutrino area for IceCube, ANTARES and for two different configurations of a 
cubic kilometer detector in the Mediterranean sea, presented in the KM3Net design report 
||6T| . The effective area represents a useful tool to compare different experiments, because 
for a given source it is straightforward to compute the number of expected events as: 

In eq. [30] depends on neutrino energy and direction. This has two consequences: 
i) sources with a similar fluency (number of neutrino per unit of area and unit of time) but 
different spectral index produce a different response in neutrino telescopes (harder is the 
spectral index, better the source is seen); 
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log,o(E[GeV]) 

Figure 16: Neutrino effective area as a function of the true simulated neutrino energy. 
Two different versions of a km^ f^ZD detector in the Mediterranean sea are compared to 
the IceCube ( ^9.3i detector in the Antarctic South Pole. The ANTARES ( mO.li neutrino 
effective area is shown for comparison. 



if) due to the Earth motion, the position in the detector frame of a given source in the 
sky changes with daytime. The effective area must be computed for each declination, by 
averaging over the local coordinates (zenith and azimuth angle). 



6.4 Number of optical sensors in a neutrino telescope 

The instrumented volume for a neutrino telescope must be at least of the order of a km^. 
How many optical sensor (PMTs) are needed? This is one of the major impact factor on the 
cost of an experiment. 

In this computation, we assume that the telescope use a PMT with a 10" diameter, 
detection area Apmt ~ 0.05 and quantum efficienc>0 epmt — 0.25, see Fig. [14] Similar 
PMTs, which have the advantage to fit inside commercial pressure-resistant glass spheres 
(optical module, OM), have been chosen by the IceCube, ANTARES, NEMO and NESTOR 
collaborations. The overall efficiency of the OM is somewhat reduced with respect to that 
of the PMT, due to the presence of the glass, glue between glass and the PMT, and mu-metal 
cage for magnetic shield: eom — O.Sepmt = 0.2. 

As we will discuss in ^7.2[ ice or seawater absorption length Xabs is larger than 50 m 



'New 10" PMTs with higher quantum efficiency, ep^t — 0.35, are also considered in the KM3NeT consor- 
tium 
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for Ught in the 400-500 nm range. A photon falling inside the effective PMT volume Vpmt= 
■Apmt X ^ahs — 2.5 Can produce a photoelectron (p.e.) with a probability eom — 0.2. 

Let us call Npmt the number of optical sensors inside the instrumented volume (it is the 
number to be determined). The rate R between the effective PMT volume of Npmt and the 
instrumented volume is: 

R = ^""^^^ = X 10-' Npmt (32) 

The total number of Cherenkov photons emitted by a 1 km length muon track in the 
wavelength range of PMTs sensitivity ( ^5.41 ) is Nj ~ 3.5 x 10''. The fraction converted 
into photoelectrons which gives a signal is: 

Np.e. = X RX Com 

~ (3.5 X 10^) . (2.5 X IQ-^Npmt) ■ eom 

= 1.8xlO-^Npmt (33) 

The number of PMTs needed to reconstruct a muon track is of the order of a few tens. 
We must take into account that in most cases many photons arrive on the same PMT during 
the integration window of the electronics (which is of the order of 20-50 ns). For this reason, 
on average around Np^e. ~ 100 p.e. are necessary to affect few tens of OM. The number of 
optical sensors in a neutrino detector follow straightforward from eq. [33] 

Npmt ^ 100/1.8 X 10^2 ^ 5Q00 (34) 

As we will shown in ^9.3[ the IceCube collaboration is in an advanced status to bury 
under the ice 4800 OMs; the neutrino telescope in the Mediterranean sea, ^TT\ plan to 
deploy between 5000 and 10000 OMs, depending on the financial budget. 



7 Water and Ice properties 

The effects of the medium (water or ice) on light propagation are absorption and scattering 
of photons. These affect the reconstruction capabilities of the telescope. In fact, absorption 
reduces the amplitude of the Cherenkov wavefront, i.e. the total amount of light on PMTs. 
Scattering changes the direction of propagation of the Cherenkov photons and the distribu- 
tion of their arrival time on the PMTs; this degrades the measurement of the direction of the 
incoming neutrino. We define direct photons those which arrive on a PMT in the Cherenkov 
wavefront, without be scattered; otherwise, we define them indirect photons. 

The propagation of light in a transparent medium is quantified for a given wavelength A, 
by the medium inherent optical properties: absorption a(A), scattering 6(A) and attenuation 
c(A) = a(A) + 6(A) coefficients, or, alternately, absorption La{X) = a(A)^^, scattering 
Lfe(A) = 6(A)^^ and attenuation Lc{X) = c(A)^^ lengths. Each of these lengths represents 
the path after which a beam on initial intensity /q at wavelength A is reduced in intensity by 
a factor of 1/e through absorption and scattering, according to the following relation: 

/i(x,A) = /o(A)e-^/^'(^); i = a,b,c (35) 
where x (in meters) is the optical path traversed by the light. 
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Figure 17: Depth dependence of scattering coefficient b(.{\) = l/Ll^J^ {X) as measured by 
the IceCube collaboration I J73] for 4 different wavelengths. 

A complete description of light scattering would require, in addition to the geomet- 
ric scattering length ^^(A), the knowledge of the scattering angular distribution. Gustav 
Mie developed (1908) an analytical solution of the Maxwell equations for scattering of 
electromagnetic radiation by spherical particles, which is appropriate for modeling light 
scattering in transparent media. In particular for ice the predominant scattering centers are 
sub-millimeter sized air bubbles and micron sized dust particles. 

Generally, light can be scattered multiple times before it reaches an optical sensor. The 
average cosine of the light field of photons that have undergone multiple (= n times) scat- 
tering obeys a simple relationship: 

(cos0)„ = (cos^)" (36) 

On average, per step, a photon advances at an angle of (cos 9) a distance of L},{X) between 
each scatter. Hence after n scatters, a photon has moved in the incident direction: 

Llff{X) = ±{cos9y 3^1^ (37) 

Experimental measurements are generally expressed in terms of the effective light scatter- 
ing length L^J^ {\), instead of the (strongly correlated) values of average scattering angle 
(cos 6) and geometric scattering length L(,(A). 

In the following discussion, we will point out that seawater has a smaller value of La {X) 
with respect to ice (which is more transparent). By referring to the discussion in %.3[ the 
same instrumented volume of ice corresponds to a larger effective volume with respect to 

seawater. On the other hand, the effective scattering length L^-^-^ for ice is smaller than 
water. This is a cause of a larger degradation of the angular resolution of the detected 
neutrino-induced muons in ice with respect to the water (see ^7.41) . 
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Another difference between ice and water is that optical modules in seawater suffer 
some background from the natural radioactivity of elements (mainly ^^K) and from lumi- 
nescence produced by organisms living in the deep sea ( m.3\i . Ice is (almost) background- 
free. 



7.1 Ice properties 



Ice model extrapolation with EDML ice core data 
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Figure 18: Depth dependence of scattering coefficient be{X = 400 nm) and absorption 
coefficient adust{^ = 400 nm) as used in Monte Carlo simulations in the Antarctica exper- 
iments hi 74\l . 

The ice in which AMANDA and IceCube (see ^9.31 ) are embedded has optical properties 
that vary significantly with depth and that need to be accurately modeled. Impurities trapped 
in the ice depend on the quality of the air present in the snows: Antarctica ice is laid down 
through a process of snowfall, hence trapping bubbles of air as it compacts itself. This 
happened over roughly the last 10^ years. Because of variations in the long-term dust level 
in the atmosphere during this period, as well as occasional volcanic eruptions, impurity 
concentrations are depth dependent. IceCube and AMANDA detectors have both pulsed 
and steady light sources located at different depths under the ice. These sources are used to 
measure both the attenuation length and scattering length. This is done by measuring the 
arrival time distributions of photons at different distances from a light source. 

The scattering centers for light propagation in IceCube are air bubbles and dust particles 
of various types 11731 . 

Air bubbles play a major role in light scattering down to depths of approximately 1250 
m below the surface. Then, the pressure of ice layers above compact these air bubbles into 
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Figure 19: The absorption length measured in the underwater sites of Toulon (blue) and 
Capo Passero (red), of the ANTARES ( m0.1\l and NEMO (see ^70.21 ) experiments, respec- 
tively. Left: La{\) with A = 440 nm as function of the depth. Right: La{X) as function of 
the wavelength, compared to the behaviour of pure seawater (solid line). 



air hydrate crystals, which have an index of refraction nearly identical to that of ice. 

Four main components of the dust were individuated: insoluble mineral grains, sea salt 
crystals, liquid acid drops and soot. Sea salt crystals and liquid acid drops contribute neg- 
ligibly to absorption, sea salt being the strongest scattering component. Insoluble mineral 
grains are the most common component, and contribute to both absorption and scattering, 
while soot contributes mainly to absorption. The relati ve abundance of each of these com- 
ponents was derived from ice core measurements II175II . 

Fig. [17] shows the effective scattering coefficient be (A) measured in Antarctica. The 
strong drop off from depths of around 1250 m is due to the transition from the region where 
air bubbles are dominant to the region where the four main dust peaks are present in the ice. 

The effects of ice properties on photon propagation and arrival times on PMTs are eval- 
uated by the IceCube collaboration through Monte Carlo (MC) simulations. MC parameters 
are adjusted until an agreement is met between simulation and real data for photon timing 
distributions. Some models were developed inside the AMANDA/ IceCube collaborations 
II174II . The depths below AMANDA, which are now included in the active volume of Ice- 
Cube, were described with a model corrected using the measurements in ice cores taken at 
other sites in Antarctica, in particular at East Dronning Maud Land (EDML). Absorption 
(due to dust) and scattering coefficients obtained with different simulations and compared 
with the measurements are shown in Fig. [18] 
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7.2 Optical properties of water 

Optical properties of water depend on many factors. Environmental parameters such as 
water temperature and salinity are indicators of the aggregation state of H2O molecules, 
which biases the diffusion of light. Water absorption and scattering depend also on the 
density and the size of the floating particulate, which affects the telescope response also in 
terms of detector aging: due to bio-fouling and sediments sticking on the optical modules, 
efficiency of the photon detection can be compromised. 

For these reasons, together with the strong necessity of reducing the atmospheric muon 
background, a site for a neutrino telescope must be located at great depth. In this condition, 
high pressure and extremely slow water currents make the site characteristics stable. 
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Figure 20: Average attenuation length measured with the AC9 in the Capo Passero site 
nW . at depth 2850-3250 mfrom year 1999 to 2003. Only statistical errors are plotted. 
The solid line indicates the values of Lc(X) for optically pure seawater reported by Smith 
and Baker H177\l . 



Apart from the BAIKAL experiment ( ^8.21) . situated in the Siberian lake Baikal at a 
depth of approximately 1 km, submarine sites have been preferred in order to reach deeper 
locations. The preference for undersea sites is not free from drawbacks: because of the salts 
into the water, they present an irreducible optical background due to the radioactive decay 
of and to the bioluminescence, which strongly depends on environmental factors. 

In order to minimize the bias induced by external agents, the telescope sites must be 
in addition far enough from continental shelf breaks and river estuaries, which can induce 
turbulent currents and spoil the purity of water. At the same time, the neutrino telescope 
should be close to scientific and logistic infrastructures on shore. With such requirements, 
the Mediterranean Sea offers optimal condition on a worldwide scale. In water the absorp- 
tion and attenuation coefficients a(A) and c(A) are directly measured by means of dedicated 
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instruments, like the AC9 manufactured by WETLabs II178II . Water optical properties are 
strongly dependent on the wavelength: light transmission is extremely favored in the range 
350-550 nm [ 179 ], where the photomultipliers used to detect Cherenkov radiation reach the 
highest quantum efficiency (see Fig. [T4l) . 

In natural seawater, optical properties are also function of water temperature, salinity 
and dissolved particulate [,18nil82]| . Measurements of the profiles of temperature, salinity, 
attenuation and absorption lengths performed by the NEMO collaboration ( ^10.21 ) in the 
site named Capo Passero, 100 km off shore from the coast of Sicily, during various sea 
campaigns from year 1999 to the end of year 20 03 sh ow that such quantities are stable 
and constant at depths greater than 1500H-2000 m II176II (see Fig. [3l]for the geographical 
location of sites). 
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Figure 21: Median rates (in kHz.) measured with the 10" PMTs of the ANTARES experi- 
ment, on optical modules at two different depths (2037 m and 2386 m). Data from March 
2006 up to May 2008 l\180t The contribution of the "^^K decay is evaluated to be almost 
constant to ~ 30 — 40 kHz. 



The nature of particulate, either organic or inorganic, its dimension and concentration, 
affect light propagation. All these environmental parameters may vary significantly, for 
each marine site, as a function of depth and time. Moreover, it is known that seasonal effects 
like the increase of surface biological activity (typically during spring) or the precipitation 
of sediments transported by flooding rivers, enlarge the amount of dissolved and suspended 
particulate, worsening the water transparency. 

Fig. [19] (left) shows the absorption length La{X = 440 nm) as a function of depth 
measured in the site of the ANTARES experiment (near Toulon, France) and in the Capo 
Passero site. The right side of the figure shows the measured absorption length (in the range 
of 350 nm < A <700 nm, circle and square dots) in the two sites, compared to the model 
of pure seawater reported by Smith and Baker U77il (solid black line). 

Fig. |20] shows the mean attenuation length Lc(A) measured in different seasons in the 
Capo Passero site; such averages are obtained with measurements from 2850 m to 3250 m 
depth. 
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7.3 Optical background in water 

The background in seawater has two main natural contributions: the decay of radioactive 
elements in water, and the luminescence produced by organisms, the so called biolumines- 
cence. 

The is by far the dominant of all radioactive isotopes present in natural seawater. 
decay channels are: 

40^ ^ 40 + e" + (BR = 80.3%) 
^^K + e- ^ ^^Ar + i/^ + 7 (BR = 10.7%) 

and both contribute to the production of optical noise. A large fraction of electrons produced 
in the first reaction is above the threshold for Cherenkov light production. The photon orig- 
inating in the second reaction has an energy of 1.46 MeV and can therefore lead (through 
Compton scattering) to electrons above the Cherenkov threshold. 

The intensity of Cherenkov light from ^'^K radioactive decays depends mostly on the 
^'^K concentration in sea water. Since salinity in the Mediterranean Sea has small geograph- 
ical variation, this Cherenkov light intensity is largely site independent. 

Bioluminescence is ubiquitous in oceans and there are two sources in deep sea: steady 
glow of bacteria and flashes produced by animals. These can give rise to an optical back- 
ground up to several orders of magnitude more intense than the one due to ^*^K (see Fig. 

The two components of optical background described above are clearly visible. Bursts 
observed in the counting rates are probably due to the passage of light emitting organisms 
close to the detector. 

The phenomena of bioluminescence are not yet fully understood. The typical spectrum 
of bioluminescence light is centered around 470-480 nm ||183[ 11841 . the wavelength of 
maximal transparency of water. The distribution of luminescent organisms in deep-sea 
varies with location, depth, and time but there is a general pattern of decrease in abundance 
with depth. Fig. [22] shows the amount of luminescent cultivable bacteria as a function of 
depth, measured in the Capo Passero site. Such measurements show a bioluminescence that 
is significantly lower with respect to what discovered at similar depth in the Atlantic Ocean 
lfT85l[l86l . 

Deep sea currents were monitored at the ANTARES, NEMO and Nestor sites for long 
time periods. ANTARES discovered that the baseline component is neither correlated with 
sea cuiTcnt, nor with burst frequency; however, long-term variations of the baseline were 
observed. Periods of high burst activity are not correlated with variations of the baseline 
component, suggesting that each of the two contributions is caused by a different popula- 
tion. Moreover, a strong correlation is observed between bioluminescence phenomena and 
the sea current velocity, as shown in Fig. |23l 

7.4 Track reconstruction in water and ice 

In a 1 km^ scale detector, most of the high energy muons produce tracks which are visi- 
ble over more than 1 km. This long lever arm allows for good directional reconstruction, 
depending on the medium (water or ice), number and orientation of the optical sensors. A 
rough estimate of the muon energies it is also possible, either by the length of their tracks, or 
by measuring the specific energy loss; at energies above 1 TeV, muon energy loss (dE/dx) 
is proportional to muon energy. 

Muon reconstruction is done by maximum likelihood methods. The fitter finds the 
likelihood for different track positions and directions, and, optionally, energy. To do this. 
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Figure 22: The amount of luminescent bacteria per unit of seawater volume, sampled in the 
Capo Passero site at different depths. The bacteria were cultived at atmospheric pressure. 



it uses functions which model the hght propagation, giving the Probabihty Distribution 
Function for a photon, radiated from a track with a given orientation, to reach a PMT at 
a given distance and orientation as a function of time. Usually, these functions are pre- 
calculated using a simulation that tracks photons through the medium 1 187, 189, 190 1. 

Fig. [24] shows the angular resolution in ice and water resulted from MC calculations. 
For neutrino-induced muons (up-going) with > 1 TeV, the directional resolution is of 
the order of few tenths of degree in water and around 1 degree in ice. and low-energy 
Ut (below 1 PeV) interactions, or neutral current interactions, occurring inside the detector 
instrumented volume produce a shower, a compact deposition of energy, ^5.11 ^5.2[ 

7.4.1 Direct and indirect photons 

The effect on the tracking algorithms due to the arrival of direct and indirect Cherenkov 
photons is presented in Fig. |25](made for a simulation in water using 1 TeV muons, II190II ). 
Direct photons are emitted by a muon almost exactly at the Cherenkov angle and they arrive 
at the PMT without be scattered by the propagation medium. They caiTy the most precise 
timing information: the arrival time is only smeared by the dispersion and the transit time 
spread of the PMTs. A muon reconstruction algorithm which can ideally use only direct 
photons, is able to obtain the particle direction with the highest precision. This can be 
proved with Monte Carlo simulations, where direct and indirect photons can be tagged. 
The difference between the expected and the measured arrival time (the time residual r) for 
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Figure 23: Correlation between the burst fraction and the seawater current velocity as 
measured by the ANTARES detector. The burst fraction is the fraction of time with count 
rates on OMs exceeding 120% of the baseline rate. 



these photons has a Gaussian shape distribution (full line in Fig. [25] ). 

Photons that originate from secondary electrons or that have scattered, are often delayed 
with respect to this time. However, also for these photons the distribution of time residuals 
peaks at zero, which means that they can still be used in the reconstruction process. This is 
the reason why the measurement of the direction of electromagnetic and hadronic showers 
in Ve CC and in Vx NC interactions is much less precise with respect to the z/^ CC interaction 
channel. 



7.4.2 The background of atmospheric muons 

The bulk (see Fig. [TOl ) of reconstructed events in any neutrino telescope are downward 
going muons produced in cosmic-ray air showers flSS]. Atmospheric muons can be used 
for a real-time monitoring of the detector status, of the time variation of the PMTs efficien- 
cies, and for detector calibration. Atmospheric muons can also be used for the study of the 
telescope pointing capability through the measurement of the moon shadow. On the other 
side, muons (especially in bundles) are a major background source: downward going parti- 
cles wrongly reconstructed as upward going and simultaneous muons produced by different 
cosmic ray primaries could mimic high energy neutrino interactions. 

Muons in bundle are very common: they represent the majority of the events which 
trigger a neutrino telescope. All the muons in a bundle are almost parallel, and propagate 
in a plane perpendicular to the shower axis at the same time. Because they are produced 
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Figure 24: (Left) Angular resolution (evaluated with MC) for 9 strings and for the full 
IceCube array as a function of event energy l{188\l . It is shown here the differences be- 
tween true and reconstructed muon track. (Right) The same, for the underwater ANTARES 
detector In this case, it is also shown the difference with respect to neutrino direction. 



in the decay of secondary mesons, they follow the transverse and longitudinal momentum 
distribution of the parent mesons. As a consequence, the most energetic muons are expected 
to be closer to the axis shower. Most of the muons in the bundle have a radial distance with 
respect to the shower axis smaller than 10 meters 111561 [1911. i.e. much smaller than the 
grid size of the detectors. For this reason, usually muon reconstruction algorithms find a 
track which is in reality due to the sum of the signals induced by the muons in the shower. 
The muon multiplicity in the bundle is not measured (although works are in progress inside 
collaborations to obtain a rough estimate of the bundle muon multiplicity). 

Because of the high rate of downward going muons, to distinguish z/-induced events it 
is not enough to select events with the most likely reconstruction as upward going. Fairly 
stringent cuts must be applied to eliminate tracks with reasonable likelihoods for being 
downward going. This can be done by cutting on estimated errors from the likelihood fit, 
or using other quality estimators (see for instance Fig. |32ll. The exact cuts depend on 
the medium (water or ice); cuts are also analysis-dependent, since different analyses are 
interested in signals from different energy ranges and zenith angles. 

Another important possible background source is due to random coincident muon 
events. This happens when two muon bundles from independent cosmic-ray air showers 
traverse the detector in the same time-window (few ^s) when the detector is read out. The 
effect depends strongly from the depth and size of the detectors. For instance, the upper 
level of IceCube ( ^9.31 ) is at the depth of ^ 1450 m and the frequency of these coincident 
events is relevant. For th e sm aller-sized AMANDA-II detector, the estimated trigger rate 
of these events is 2-3 Hz [112], compared with the total trigger rate of 80 Hz. Running (or 
proposed) underwater experiments are deeper and this effect is largely reduced |192|. In 
IceCube, specific algorithms have been developed to find and reject coincident events, by 
separating hits from the two tracks based on their separation in space and/or time. 
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Figure 25: Distribution of the difference r between the expected and measured photon ar- 
rival times for 1 TeV muon in water A measurement of the light intensity and arrival time 
in one PMT is denoted as "hit". Contributions are shown for scattered and direct photons 
originating from the muon itself and from secondary electrons and positrons. All hits occur- 
ring within a distance of 100 mfrom the track are included in the figure. The distribution 
for background photons corresponds to a flat background rate of 60 kHz.. Adapted from 



8 The pioneers: DUMAND and Lake Baikal experiments 
8.1 The prototype: DUMAND 

The project to realize the Markov idea started in 1973 during the Cosmic Ray Conferen ce in 
Den ver. The Deep Underwater Muon And Neutrino Detection (DUMAND) project 111931 
11941 born in 1976 and existed through 1995. The goal was the construction of the first 
deep ocean neutrino detector, to be placed at a 4800 m depth in the Pacific Ocean off 
Keahole Point on the Big Island of Hawaii. Many preliminary studies were carried out, from 
technology to ocean water properties. A prototype vertical string of instruments suspended 
from a special ship was employed to demonstrate the technology, and to measure the cosmic 
ray muon flux at various depths (2000-4000 m, in steps of 500 m) in deep ocean 131 • 

A major operation took place in December 1993, when one string of photo-detectors, 
a string of environmental instruments and a junction box were placed on the ocean bottom 
and cabled to shore. While the cable laying was successful, short circuits soon developed in 
the instruments and it was no longer possible to communicate with the installed apparatus. 
In 1995 the US DOE cancelled further efforts on DUMAND. 

All subsequent designs for underwater experiments have taken advantage of this expe- 
rience. Some reasons for the long DUMAND development time were: i) huge depth of the 
chosen site; ii) lack of advanced fibre-optics technology for data transmission; in) lack of 
reliable pressure-resistant underwater connectors; iv) lack of Remotely Operated Vehicle 
(ROV) for underwater connections; v) limited funding. 

The Baikal group has been working in Lake Baikal in Siberia for about the same time 
as the DUMAND group in Hawaii. Initially, in the mid-1970's, the groups worked together. 
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However, political problems developed after the Soviet invasion of Afghanistan, and the 
US DUMAND group was told by its government (Reagan administration) that no funding 
would be available to collaborate with the Soviet groups. Hence the teams reluctantly took 
separate paths. 

8.2 The experiment in Lake Baikal 

The possibility to build a neutrino telescope in the Russian Lake Baikal was born with 
the basic idea of using the winter ice cover as a platform for assembly and deployment of 
instruments, instead of using a ship |195|. After initial small size tests, in 1984-90 single- 
string arrays equipped with 12 - 36 PMTs were deployed and operated via a shore cable. 
During this period, underwater and ice technologies were developed, optical properties of 
the Baikal water as well as the long-term variations of the water luminescence were in- 
vestigated in great detail. Deep Baikal water is characterized by an absorption length of 

La(480nm) = 20 24 m, an effective scattering length of V^J^ = 30 H- 70 m and a strongly 
anisotropic scattering function with a mean cosine of scattering angle (cos 6) = 0.85 0.9. 

The Baikal Neutrino Telescope NT- 200 was a second generation detector, deployed in 
Lake Baikal 3.6 km from shore at a depth of 1.1 km. It consists of 192 optical modules 
(OMs). In April 1993, the first part of NT-200, the detector NT-36 with 36 OMs at 3 short 
strings, was set into operation. A 72-OMs array (NT-72) ran in 1995-96. In 1996 it was 
replaced by the four-string array (NT-96). Since April 1997 a six-string array with 144 
OMs, take data in Lake Baikal (NT- 144). NT-200 array was completed in April, 1998 (Fig. 
|26l ). NT200 plus the new external strings form NT200-I-. An umbrella like frame carries 

8 strings 72 m long, each with 24 pair of optical modules (OMs). The OMs contain 37 
cm diameter PMT QUASAR370 and are grouped in pairs along the strings. The PMs of 
a pair are switc hed in coincidence in order to suppress background from bioluminescence 
and PMT noise |[T96l . 

The search strategy for high energy neutrinos relies on the detection of the Cherenkov 
light emitted by cascades, produced by neutrino interactions in a large volume below NT- 
200. The results of a search for high energy f g + f ^ + Vj- Ulllll gives an upper limit of 
E-'^d^/dE < 8.1 X 10"'^ GeV cm"^ ^-i sr^^ in the energy range of 2x10^ to 5 x 10^ 
GeV, see Fig. El 

9 Detectors in the South Pole ice 

An experiment at the South Pole, at the Amundsen-Scott station where the ice is about 2800 
m deep, was pioneered by the AMANDA collaboration. They drilled holes in the ice using 
a hot water drill, and lowered strings of optical sensors before the water in the hole refreeze. 
The first AMANDA string was deployed in 1993, at a depth of 800-1000 m. It was quickly 
found that at that depth the ice had a very short scattering length, less than 50 cm (Fig. [TT] ). 
In 1995-96 AMANDA deployed 4 strings between 1500 and 2000 m deep. These detectors 
worked as expected, and AMANDA detected its first neutrinos [5]. This success led to 
AMANDA-II, which consisted of 19 strings holding 677 optical sensors. 

AMANDA was limited by its small size and some technological problems II197II . Its op- 
tical sensors consisted of photomultipliers (PMTs) with resistive bases in a pressure vessel. 
Not all of the optical modules survived the high pressures present when the water in the drill 
holes froze and AMANDA consumed considerable electrical power and required manpo- 
wer-intensive calibrations yearly. High voltage was generated on the surface, and analogue 
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To Shore ^ 




Figure 26: Sketch of the NT-200 Baikal experiment. The expansion left-hand shows 2 pairs 
of optical modules (svjaska) with the svjaska electronics module, which houses part of the 
readout and control electronics. 



signals were returned to the surface. Several transmission media were tried: coaxial cables, 
twisted pairs, and optical fibers. The 2.5 km long coaxial cables and twisted pairs dispersed 
the PMT pulses, and also the optical fibers (in roughly half of the OMs) had a very limited 
dynamic range. 

IceCube was designed to avoid these problems and to be much simpler to deploy, op- 
erate and calibrate. When completed in 2011, it will consist on a deep detector (Inlce) and 
a surface detector (IceTop), see Fig. [271 The design of the main Inlce part of the detector 
||T98 1 consists of 80 strings, buried 1450 to 2450 meters under the surface of the ice, each 
bearing 60 Digital Optical Modules (DOMs), with 17 m spacing. The strings are placed on 
a 125 m hexagonal grid, providing a 1 km^ instrumented volume. The surface electronics 
are in a counting house located in the center of the array. 

The IceTop surface air-shower is an array of 80 stations II 19911 for the study of extensive 
air showers. Each IceTop station, located above an IceCube string, consists of two tanks 
filled with ice. Each of those tanks contains two DOMs of same design as the one used for 
the Inlce part of the detector. The surface array can be operated looking for anti-coincidence 
with the Inlce events to reject downgoing muons. It can also be used in coincidence, to 
provide a useful tool for cosmic ray composition studies. The array covers an area of about 
1 km^ as shown in Fig. [27] The completed detector will be operated for 20 years. 

Because of the Antarctic weather, high altitude and remote location of the South Pole, 
logistics is a key issue. The construction season lasts from November through mid-February 
(during the Austral summer). Everything needed must be flown to the Pole on ski-equipped 
LC-130 transports planes. The main task in IceCube construction is drilling holes for the 
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Deep Core 



Figure 27: The IceCube detector side view. Currently there are 59 buried Inlce strings. The 
AMANDA detector appears in the right part of Inlce. The IceTop surface array and the 
DeepCore (^9.4) are also shown. 



strings, see Fig. |28] This is done with a 5 MW hot-water drill, which generates a stream 
of ~ 800 liters/minute of 88°C water. This water is propelled through a 1.8 cm diameter 
nozzle, melting a hole through the ice. Drilling a 2500 m deep, 60 cm diameter hole takes 
about 40 hours. Deploying a string of DOMs takes about another 12 hours. 

Data acquisition with the partially finished IceCube detector is running smoothly and 
the detector is operating as expected. The detector began taking data in 2006 with a nine- 
string configuration (IC-9) and with a 22-string configuration in 2007 (IC-22) II200I . Dur- 
ing the Austral summer 2007/08 18 more strings were buried, and 19 during the summer 
2008/09. Currendy (May 2009) 59 Inlce strings with 3540 DOMS are deployed. 

9.1 The IceCube Data Acquisition System 

Each DOM used by IceCube comprises a 10" PMT (Hamamatsu R7081-02) housed in a 
glass pressure vessel and data acquisition (DAQ) electronics which reads out, digitizes, 
processes and buffers the signals from the PMT. When individual trigger conditions are met 
at the DOM, it transfers fully digitized waveforms to a software-based trigger and event 
builder on the surface. The electronics acquire in parallel, on Analog Transient Waveform 
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Figure 28: The hot water hose and support cables disappear down one of one of boreholes 
drilled into the Antarctic ice to construct the IceCube Neutrino Observatory. Photo: Jim 
Haugen l[20J]l 



Digitizers (ATWDs) at 300 MHz, sampling over a 425 ns window. In addition the electron- 
ics also record the signal with a coarser 40 MHz sampling over a 6.4 fis window to record 
the late part of the signals. Two parallel sets of ATWDs on each DOM operate to reduce 
the dead-time: as one is active and ready to acquire, the other is read out. 

Data is transmitted to the surface via a single twisted copper cable pair, which also pro- 
vides power. Each DOM consumes about 3.5 W. The cable also includes local coincidence 
circuitry, whereby DOMs communicate with their nearest neighbors. A more robust con- 
nector is used than in AMANDA, and a higher fraction of IceCube OMs survive the freezing 
of the ice. The main requirement for the IceCube hardware is high reliability without main- 
tenance. Once deployed, it is impossible to repair a DOM. About 98% of the DOMs survive 
deployment and freeze-in completely; another 1% have lost their local coincidence connec- 
tions, but they are usable. On the surface, the cables are connected to a custom PCI card in 
a PC; the remainder of the system is off-the-shelf. 

Each DOM also contains a 'flasher' board, which has 12 blue (405 nm) LEDs mounted 
around its edges. These LEDs aie used for calibrations, to measure light transmission and 
timing between different DOMs, to check the DOM-to-DOM relative timing and study the 
optical properties of the ice. The time calibration yields a timing resolution with a RMS 
narrower than 2 ns for the signal sent by the DOM to the surface II202I . The noise rate 
due to random hits observed for Inlce DOMs is of the order of 300 Hz, which gives the 
possibility to monitor the DOM hit rates. This very low value makes the detector able to 
have a sensitivity to low energy (MeV) neutrinos from supernova core collapse in the Milky 
Way and in the Large Magellanic Cloud |203 |. 

Data from the DOMs are time-sorted, combined into a single stream, and then mon- 
itored by a software trigger. The main trigger is based on multiplicity: it requires eight 
DOMs (with local coincidences) fired within 5 /iS. This collects most of the neutrino events. 
Starting from 2008, a string trigger which improved sensitivity for low energy, requiring five 
out of seven adjacent DOMs fired within 1.5 /is, was added. When any trigger occurs, all 
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data within the zblO^us trigger window is saved, becoming an event. If multiple trigger 
windows overlap, then all of the data from the ORed time intervals are saved as a single 
event. 

The total trigger rate for 40 strings was about 1.4 kHz. The majority of the triggers 
(about 1 kHz) are due to atmospheric muons. A fast on-line filter system reduces the trig- 
gered events (6% survives, for a data rate of ~30 Gbytes/day), and selected events are 
transmitted via satellite to the Northern hemisphere. The rest of the data is stored on tapes 
at the South Pole station, and tapes are carried North during the Austral summer. 

9.2 Summary of the AMANDA-II results 




Reconstructed Zenith Angle 6 (deg) Reconstructed Azimuth Angle ii> (deg) 

Figure 29: Reconstructed zenith (on the left) and azimuth (on the right) angle distributions 
for the final sample ofIC-9 events. A zenith angle of 90" (180°) corresponds to horizontal 
(straight up-going) event. The shadowed area indicates expectations with systematic errors. 
The error bars are statistical only. The configuration of the IceCube strings seen from the 
top is also plotted (box on top of the right). The preferred axis of this configuration explains 
the features observed in the azimuth angle distribution. 

IceCube has also integrated its predecessor, the AMANDA-II detector (the final con- 
figuration of the AMANDA detector as an independent entity). AMANDA, with 677 OMs 
distributed on 19 strings, is now surrounded by IceCube (see Fig. [27] ). For relatively low 
energy events, the dense configuration of AMANDA gives it a considerable advantage. 
Moreover, IceCube strings surrounding AMANDA can be used as an active veto against 
cosmic ray muons, making the combined IceCube -i- AMANDA detector more effective for 
low energy studies than AMANDA alone. 

AMANDA-II has been taking data between 2000 and 2004 and has collected 4282 up- 
going neutrino candidates with an estimated background fraction of --^ 5% from wrongly 
reconstructe d dow nward going muons. The analysis for point sources in the Northern hemi- 
sphere sky p204'| yielded no statistically significant point source of neutrinos. Assuming 
as usual a source of muon neutrinos with energy spectrum of E^'^, an upper limit was 
placed averaged over declination in the Northern hemisphere sky at 90% confidence level: 
E~'^d^/dE < 5.5 X 10"^ GeV cm^^ s"^ in the energy range from 1.6 TeV up to 2.5 PeV. 
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Over the same period of time, a search f or neu trino emission from 32 pre-selected 
specific candidate sources has been performed f2041. No statistically significant evidence 
for neutrino emission was found, see Fig. [34] The highest observed significance, with 8 
observed events compared to 4.7 expected background events, is at the location of the GeV 
blazar 3C273. 

In addition to searches for individual sources of neutrinos, AMANDA-II data taken be- 
tween 2000 and 2003 have been used to set a limit on possible diffuse flux of neutrinos. As 
described in ^3.3[ this diffuse flux can be distinguished from the background of atmospheric 
neutrinos due to its harder spectrum. This study relies on the number of triggered OMs 
which serve as an energy estimator. A limit of E~'^d^/dE < 7.4 x 10^^ GeV cm~^ s^^ 
sr^^ was placed (see Fig. IT]) on the diffuse muon neutrino in the energy range from 16 TeV 
to 2.5 PeV at 90% confidence level 1 112, 209|. Additionally, AMANDA-II has searched for 
an all-flavor diffuse flux from the Southern sky, setting a limit of E~'^d^/dE < 2.7 x 10^^ 
GeV cm"2 ^-i the energy range of 2 x 10^ to 10*^ GeV |[TT3l . 



9.3 First results from the IceCube 9 and 22 strings configuration 

The IC-9 dataset has a total livetime of 137.4 days taken between June and November 2006. 
234 neutrino candidates were identified on this data sample with 21 l±76sy st i I4:stat events 
expected from atmospher ic ne utrinos and less than 10% pollution by the background of 
downward going muons 11205 1 . Zenith and azimuth angle distributions of these neutrino 
candidates are shown on Fig. |29] Agreement with simulation is good except for a discrep- 
ancy near the horizontal direction due to a residual contamination of down going muons. 




Figure 30: IC-22 skymap with pre-trial p-values (in colors) and event locations (dots). 
Top: binned method. Bottom: unbinned method. Each method uses a different direction 
reconstruction technique l[206\l . 
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The IC-22 detector took data in 2007-2008, with a hfetime of 276 days ll206l . This 
data has been analyzed to search for extraterrestrial point sources of neutrinos using two 
methods: the binned and the unbinned maximum likelihood method. The binned method 
distinguishes a localized excess of signal from a uniform background using a circular an- 
gular search bin. The search bin radius depends on declination, and the mean value is 2.1°. 
The unbinned maximum likelihood method [207] constructs a likelihood function which 
depends on signal Probability Density Function (PDF) and background PDF, for a given 
source location and total number of data events. Each analysis has followed its own event 
selection criteria, arriving at a final sample of 5114 (2956) events for the unbinned (binned) 
method. From simulation, a sky-averaged median angular resolution of 1.4° is estimated 
for signal neutrinos with E~'^ spectrum. 

The results from the all-sky search for both analyses are shown in Fig. [30l The best 
sky-averaged sensitivity (90% C.L.) is E^'^d^/dE < 1.3 x 10"^ GeV cm'^ s^^ to a 
generic E^^ spectrum of over the energy range from 3 TeV to 3 PeV. No neutrino point 
sources are found from the directions of a pre-selected catalogue nor in a search extended 
to the northern sky. Limits are improved by a factor of two compared to the total statistics 
collected by the AMANDA-II detector and by IC-9 |208 |. A search for neutrinos coming 
from 26 galactic and extragalactic pre-selected objects has also been performed on this 
dataset, with null results. 



9.4 The future 

During the coming years, IceCube will continue to grow starting from present configuration 
of 59 Inlce strings. The capabilities of IceCube will be extended at both lower and higher 
energies in the near future. A compact core of 6 strings using IceCube DOM technology, 
called the DeepCore detector, will be deployed near the center of the main Inlce detector. 
Fig. [27j The inter-string spacing will be of the order of 72 m, allowing for the exploration 
of energies as low as 10-20 GeV. Surrounding IceCube strings will be used as an active veto 
to reduce the atmospheric muon background. The energy range that will be explored is very 
important for dark matter searches that were initiated with AMANDA 112 10112 111 . 

Moreover, the ability to select contained events opens the search for downgoing astro- 
physical neutrino signals at low energies. This will allow looking above the current horizon 
of IceCube, even opening the possibility of looking at the galactic center or at the RX 
J1713.7-3946 source |212|. 

At EeV energies a possible extension of IceCube is also being studied in order to in- 
crease the detection volume and to be sensible to the GZK neutrino using radio and/or 
acoustic detection of the signals generated by neutrino interacting in a huge volume of ice. 
With attenuation lengths of the order of a kilometer for acoustic (kHz frequency range) 
and for radio signals (MHz frequency range), a sparse instrumentation will be sufficient 
for this extension. Two projects are currently being explored: AURA (Askarian Underice 
Radio Array) for the radio signal |213| and SPATS (South Pole Acoustic Test Setup) for 
the acoustic signal. They are currently studying the polar ice and developing the hardware 
necessary to build a hybrid detector enclosing IceCube. 
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10 The underwater neutrino projects in the Mediterranean Sea 
10.1 The ANTARES experiment 

The ANTARES project ll2T4l has been set up in 1996 fTB). Today it involves about 180 
physicists, engineers and sea-science experts from 24 institutes of 7 European countries. 
ANTARES is at present the largest neutrino observatory in the Northern hemisphere, which 
represents a privileged sight of the most interesting areas of the sky like the galactic centre, 
where many neutrino source candidates are expected. 




Figure 31: Schematic view of the ANTARES detector 



From 1996 to 1999 an extensive R&D program has been successfully performed to 
prove the feasibility of the detector concept. Site properties have been studied such as: op- 
tical properties of the surrounding water II216II ; biofouling on optical surfaces 112171 ; optical 
backgrounds due to bioluminescence and to the decay of the radioactive salts present in 
seawater |218|; geological characteristics of its ground. These studies lead to the selection 
of the current site, 40 km off La Seyne-sur-Mer (Fran ce) at a 2475 m depth. 

A mini-instrumented line worked i n 20 05 [2191. The first detector line (also called 
string) was connected in March 2006 112201 . and the second line in September 2006. In 
July 2007, three more lines were connected, and data acquisition with five lines lasted up 
to December 2007. In this month, five more strings plus a dedicated instrumented line for 
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monitoring environmental quantities were connected. Finally, the detector was completed 
with the connection of the last two strings on May, 30th 2008. The strings are made of 
mechanically resistant electro-optical cables anchored at the sea bed at distances of about 70 
m one from each other, and tensioned by buoys at the top. Fig. |3T] shows a schematic view 
of the detector array indicating the principal components of the detector. Each string has 25 
storeys, each of them contains three optical modules (OM) and a local control module for 
the corresponding electronics. The OMs are arranged with the axis of the PMT 45° below 
the horizontal. In the lower hemisphere there is an overlap in angular acceptance between 
modules, permitting an event trigger based on coincidences from this overlap. 

On each string, and on the dedicated instrumented line, there are different sensors and 
instrumentation (LED beacons, hydrophones, compasses/tiltmeters) for timing and position 
calibration. The first storey is about 100 m above the sea floor and the distance between 
adjacent storeys is 14.5 m. The instrumented volume corresponds to about 0.05 km'^. 

The basic unit of the detector is the optical module (OM), consisting of a photomulti- 
plier tube, various sensors and associated electronics, housed in a pressure-resistant glass 
sphere Ii221l . Its main component is a 10" hemispherical photomultiplier model R7081-20 
from Hamamatsu (PMT) glued in the glass sphere with optical gel. A /x-metal cage is used 
to shield the PMT against the Earth magnetic field. Electronics inside the OM are the PMT 
high voltage power supply and a LED system used for internal calibration. 

The total ANTARES sky coverage is 3.57r sr, with an instantaneous overlap of 0.57r sr 
with that of the IceCube experiment. The galactic centre will be observed 67% of the day 
time. 



10.1.1 The ANTARES Data Acquisition System 

The Data acquisition (DAQ) system of ANTARES is extensively described in (222). The 
PMT signal is processed by an ASIC card (the Analogue Ring Sampler, ARS) which mea- 
sures the arrival time and charge of the pulse. On each OM, the counting rates exhibit a 
baseline dominated by optical background due to sea-water '^''K and bioluminescence com- 
ing from bacteria, as well as bursts of a few seconds duration, probably produced by biolu- 
minescent emission of macro-organisms, ^7.3[ Fig. |2T] shows the counting rates recorded 
by two OMs located on different storeys during the 2006-2008 runs. The average counting 
rate increases from the bottom to shallower layers. The baseline is normally between 50 to 
80 kHz. 

Differently from the ^'^K background, bioluminescence suffers from seasonal and an- 
nual variations, see Fig. [2T] There can be large variations of the rate, reaching hundreds 
of kHz in some periods. Since September 2006 to December 2008 the mean counting rate 
is 75% of the time below 100 kHz. A safeguard against bioluminescence burst is applied 
online by means of a high rate veto, most often set to 250 kHz. 

The OMs deliver their data in real time and can be remotely controlled through a Gb 
Ethernet network. Every storey is equipped with a Local Control Module (LCM) which 
contains the electronic boards for the OM signal processing, the instrument readout, the 
acoustic positioning, the power system and the data transmission. Every five storeys the 
Master Local Control Module also contains an Ethernet switch board, which multiplexes 
the DAQ channels from the other storeys. At the bottom of each line, the Bottom String 
Socket is equipped with a String Control Module which contains local readout and DAQ 
electronics, as well as the power system for the whole line. Both the Master Local Control 
Modules and the String Control Modules include a Dense Wavelength Division Multiplex- 
ing (DWDM) system. The DWDM is used in data transmission to merge several IGb/s 
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Ethernet channels on the same pair of optical fibres, using different laser wavelengths. The 
lines are linked to the junction box by electro-optical cables which were connected using a 
unmanned submarine. A standard deep sea telecommunication cable links the junction box 
with a shore station where the data are filtered and recorded. 

All OMs are continuously read out and digitized information (hits) sent to shore. In 
ANTARES, a hit is a digitized PMT signal above the ARS threshold, set around 1/3 of the 
single photoelectron level (Level hits, LO). On-shore, a dedicated computer farm performs 
a global selection of hits looking for interesting physics events (DataFilter). This on-shore 
handling of all raw data is the main challenge of the ANTARES DAQ system, because of 
the high background rates. The data output rate is from 0.3 GB/s to 1 GB/s, depending on 
background and on the number of active strings. Particular conditions define a subset of LO 
for triggering purpose. This subset (called Level 1 hits, or simply LI) corresponds either 
to coincidences of LO on the same OM triplet of a storey within 20ns, or to a single high 
amplitude LO (typically > 3 p.e.). The DataFilter processes all data online and looks for 
a physics event by searching a set of correlated LI hits on the full detector on a ~ 4 /xs 
window. In case an event is found, all LO hits of the full detector during the time window 
are written on disk, otherwise the hits are thrown away. 

The trigger rate is between 1 to 10 Hz, depending on the number of strings in data 
acquisition. Most of the triggered events are due to atmospheric muons, successive ly re - 
constructed by track-finding algorithms. If ANTARES rece ives ex ternal GRB alerts II222I . 
all the activity of the detector is recorded for few minutes, i ]3.2.2l In addition, untriggered 
data runs were collected on a weekly base. This untriggered data subset is used to monitor 
the relative PMT efficiencies, as well as to check the timing within a storey, using the ^'^K 
activity. 




Figure 32: (a) Zenith and (b) azimuth distributions of reconstructed tracks. Black points 
represent data. Lines refer to MC expectations, evaluated with two different simulation. 
The shadowed band (applied on one of the MC simulation) represents the systematic error 
f~ ±30%) due to environmental and geometrical parameters. No quality cuts (QC) are 
applied to reduce the wrongly-reconstructed upward going tracks in the region between 0° 
and 90° of zenith angle. 
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10.1.2 Time and positioning calibration systems 

Differently from the strings in ice, the ANTARES hnes are flexible and move with the sea 
current, with displacements of the order of a few meters at the top for a typical sea current of 
5 cm/s. The reconstruction of the muon trajectory is based on the differences of the arrival 
times of the photons between OMs. The ANTARES detector is expected to achieve an 
angular resolution of < 0.3° for muon events above 10 TeV, through timing measurements 
with precision of the order on the ns [223 ]. This requires the knowledge of the OMs position 
with a precision of ~ 10 cm (Ught travels 22 cm/ns in water). Pointing accuracy thus 
is limited by: i) precision with which the spatial positioning and orientation of the OM 
is known; ii) accuracy with which the arrival time of photons at the OM is measured; 
in) precision with which local timing of individual OM signals can be synchronized with 
respect to each other. 

The positions of the OMs are measured on a real-time, typically once every few min- 
utes, with a system of acoustic ti^ansponders and receivers on the lines and on the sea bed, 
together with tiltmeters and compasses. The shape of each string is reconstructed by per- 
forming a global fit based on these information. Additional information needed for the 
line shape reconstruction are the water current flow and the sound velocity in seawater, 
which are measured using different equipments: an Acoustic Doppler Current Profiler; a 
Conductivity-Temperature-Depth sensors; a Sound Velocimeter. 

Relative time resolution between OMs is limited by the transit time spread of the sig- 
nal in the PMTs (about 1.3 ns) and by the scattering and chromatic dispersion of light in 
seawater (about 1.5 ns for a light propagation of 40 m). The electronics of the ANTARES 
detector is designed to contribute less than 0.5 ns to the overall time resolution. 

Complementary time calibration systems are implemented to measure and monitor the 
relative times between different components of the detector within the ns level. These time 
calibrations are performed by: 

i) the internal clock calibration system. It consists of a 20 MHz clock generator on shore, 
a clock distribution system and a clock signal transceiver board placed in each LCM. The 
system also includes the synchronization with respect to universal time, by assigning the 
GPS timestamp to the data. 

ii) The internal Optical Module LEDs: inside each OM there is a blue LED attached to the 
back of the PMT. These LEDs are used to measure the relative variation of the PMT transit 
time using data from dedicated runs. 

Hi) The Optical Beacons [224], which allows the relative time calibration of different OMs 
by means of independent and well controlled pulsed light sources distributed throughout 
the detector. 

iv) Several thousands of down-going muon tracks are detected per day. The hit time resid- 
uals of the reconstructed muon tracks (see Fig. |25] ) are used to monitor the time offsets of 
the OM, enabling an overall space-time alignment and cahbration cross-checks. 

10.1.3 Measurement of atmospheric muons and atmospheric neutrinos 

Atmospheric muons were an important tool to monitor the status of the detector and to 
check the reliability of the simulation tools and data taking. In ANTARES, two different 
Monte Carlo (MC) simulations are used to simulate atmospheric muons: one based on a 
full Corsika simulation II225II . and another based on a parameterization of the underwater 
muon flux [226 1. 

The full MC simulation [163 1 is based on Corsika v.6.2, with the QGSJET ll227l pack- 
age for the hadronic shower development. Muons are propagated to the detector using the 
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MUSIC II228II code, which includes all relevant muon energy loss processes. 

The second MC data set is generated using parametric formulas 1 156], obtained with 
a full MC tuned in order to reproduce the underground MA CRO flux 1.229. .230 1. energy 
spectrum 112311 12321 and distance between muons in bundle ||233i The characteristics of 
underwater muon events (flux, multiplicity, radial distance from the axis bundle, energy 
spectrum) are described with multi-parameters formulas in the range 1.5 5.0 km w.e. and 
up to 85° for the zenith angle. Starting from this parametrization, an event generator (called 
MUPAGE) was developed [191] in the framework of the KM3NeT project [9] to generate 
underwater atmospheric muon bundles. 

In both simulations, muons that enter inside the surface of a virtual underwater cylin- 
der (the can) are propagated using a GEANT-based program. The can defines the limit 
inside which charged particles in MC codes produce Cherenkov photons 111631 ). Then, the 
background (extracted from real data) is added and the events are feed to a program which 
reproduces the DataFilter trigger logic. After this step, simulated data have the same format 
of the real ones. 
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Figure 33: Distribution of the sinus of the zenith angle of reconstructed events (black 
crosses): 173 days of active time with 10 to 12 lines. Blue line: MC simulation of at- 
mospheric neutrinos. Red line: atmospheric muons. The region between — 1 < sm.9 < 
corresponds to upgoing particles. A contamination form atmospheric muons is present near 
the horizontal direction. From H226\l . 



The main advantage of the full MC simulation (which is very large CPU time consum- 
ing) is the possibility of re-weighting the events according to any possible primary Cosmic 
Rays flux model. The main advantage of the MUPAGE simulation is that a large sample is 
produced with a relatively small amount of CPU time (much less than the time needed to 
simulate the Cherenkov light inside the can), and it is particularly suited for the simulation 
of the background for neutrino events. 

Fig. [32]shows the zenith and azimuth distributions of reconstructed muon tracks. Black 
points represent experimental data. The solid [234| and the dotted |235| lines refer to 
Monte Carlo (MC) expectations obtained using the full MC simulation and two different 
CR composition models. The dashed-dotted line refers to the fast simulation II191II . The 
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shadowed band gives an estimate of the systematic errors, due to the uncertainties on the 
environmental parameters, like water absorption and scattering lengths in the ANTARES 
site, and on the geometrical characteristics of the detector In particular, given the fact that 
OMs are pointing downwards, at an angle of 135" w.r.t. the vertical, knowledge of the 
OMs angular acceptance at these large angles is critical for an accurate determination of the 
muon flux. 

A different analysis is necessary when selecting neutrinos. A set of more severe quality 
cuts must be applied in order to remove downward-going tracks wrongly reconstructed as 
upward. Data presented in Fig. [33] [226] were collected during the 10-12 line configuration 
period, from December 2007 to December 2008. Atmospheric neutrino events are simu- 
lated using the Bartol flux [157|. Only events detected at least by two lines and with at 
least 6 floors are considered. Restricting to the upward-going hemisphere (neutrino candi- 
dates) the number of events are 3.4 per day for data, and 2.9 per day for simulations. The 
shadowed band represents the sum of theoretical and systematic uncertainties. 

Upper flux limit from the direction of selected candidate sources were also evaluated 
using still more stringent criteria for the selection of upward-going muons II236II . The data 
with 5 lines were used and 140 active days. Even with less than half a detector, these limits 
are the best ones for experiments looking at the Southern hemisphere; they are shown as a 
function of the declination of the sources, and are compared with other experiments in Fig. 
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Figure 34: The sensitivity of the ANTARES experiment for point-like sources l[236\l . The 
results equivalent to one year of data taking as a function of declination is presented as 
blue dotted line. Magenta squares are 90% c.l. upper limits obtained by MACRO [237]. 
Red circles are similar limits by AMANDA-II [204], while the red dotted line is the average 
limit. The black dotted line shows the one year expected sensitivity for IceCube. 
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The energy of the crossing muon or of secondary particles generated by neutrino in- 
teractions inside the instrumented volume is estimated from the amount of light deposited 
in the PMTs. Several estimators based on different techniques were developed [238]. MC 
studies show that this resolution is between logio{aE/E) = 0.2 ^ 0.3 for muons with en- 
ergy above 1 TeV. The event energy measurement is a mandatory requirement for the study 
of the diffuse flux of high energy neutrinos. MC simulations indicate that after 3 years of 
data taking ANTARES can set an upper hmit for diffuse fluxes of E"^^ < 3.9 x 10^® GeV 
cm^^ s^^ sr^^ (see Fig. |71l. 

10.2 The NEutrino Mediterranean Observatory 

The NEutrino Mediterranean Observatory (NEMO) is a project P 2391 12401 of the Italian 
National Institute of Nuclear Physics (INFN). The activity has been mainly focused on the 
search and characterization of an optimal site for the detector installation; on the develop- 
ment of key technologies for the km^ underwater telescope; on a feasibility study of the 
km'^ detector, which included the analysis of all construction and installation issues and 
optimization of the detector geometry by means of numerical simulations 12411 . 

The validation of the proposed technologies via an advanced R&D activity, the proto- 
typing of the proposed technical solutions and their relative validation in deep sea is carried 
out with two pilot projects NEMO Phase- 1 and Phase-2. 

Since 1998, the NEMO collaboration conducted more than 20 sea campaigns for the 
search and the characterization of an optimal site where to install an underwater neutrino 
telescope. A deep site with proper features in terms of depth and water optical properties 
has been identified at a depth of 3500 m about 80 km off-shore from Capo Passero (36° 16' 
N 16° 06' E), see S]l2]i242i. 

The main feature of a km^ telescope is its modularity. The proposed NEMO basic 
element is the instrumented NEMO-tower (see Fig. [35l) : it is about 700 m high, and it is 
composed of 16 floors, 40 m spaced; each floor is rotated by 90°, with respect to the upper 
and lower adjacent ones, around the vertical axis of the tower. Each floor is equipped with 
two OMs (one down-looking and one horizontally looking) at both extremities. In addiction 
to the OMs, the tower hosts several environmental instruments plus the hydrophones for the 
acoustic positioning system. The tower structure is anchored at the sea bed and it is kept 
vertical by an appropriate buoyancy on the top. 

The NEMO Phase- 1 project allowed a first validation of the technological solutions 
proposed for the km'^ detector |240|. The apparatus included prototypes of all the critical 
elements: the Junction Box and a reduced version (one fourth) of the tower, called the 
mini-tower. On December 2006, both the Junction Box and the mini-tower were deployed 
and successfully activated at a test site at 2000 m depth near the Catania harbour. The 
underwater detector was connected to the shore station via a 28 km electro-optical cable. 

The NEMO Phase- 1 Junction Box was built following the concept of double contain- 
ment. Pressure resistant steel vessels were hosted inside a large fiberglass container, which 
was filled with silicon oil to compensate the external pressure. This solution has the advan- 
tage of decoupling the two problems of pressure and corrosion resistance. The electronic 
components capable of withstanding high pressures were installed directly in the oil bath. 

The mini-tower was equipped with sixteen 10" Hamamatsu R7081-SEL PMTs, 
mounted on 15 m long floors. The floors were spaced 40 m one from the other, with an 
additional spacing of 150 m from the base. 

In addition to the OMs, the instrumentation included several sensors for calibrations and 
environmental monitoring (see Fig. [36l upper part). In particular two hydrophones were 
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Figure 35: The key-elements of the NEMO proposed version for the km^ detector: the 
Junction Box, for power and data distribution, and the Tower with the bars and OM. 



mounted on the tower base and at the extremities of each floor. These, together with acoustic 
beacons placed at the tower base and on the seabed, were used for the determination of the 
OM positions with a precision of about 10 cm. Moreover, a time calibration system was 
linked to the OMs through external optical fibers. The time offsets were determined with 
the precision of 1 ns. 

One important technical choice of the design of the data acquisition system for NEMO 
Phase- 1 is the scalability to a much bigger apparatus II243I . The electric signal from the 
PMTs was sampled with 2 Flash-ADC (100 MHz each) staggered by 5 ns, for a total 200 
MHz sampling and a low power consumption. Each PMT generating an over threshold 
pulse {hit) was characterized by its time-stamp, total integrated charge and sampled sig- 
nal waveform. The latter allowed an off-line reconstruction of the hit time with precision 
of ~1 ns 1 244 1 . The hits from the four PMTs of each floor were continuously collected 
by the Floor Control Module (FCM) boards, converted into optical signals by an electro- 
optical transceiver and sent to shore through one of the optical-fiber of the 28 km cable 
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Figure 36: Up: sketch of the Phase- 1 mini-tower: the OMs on the floors and other en- 
vironmental instruments: the Acoustic Doppler Current Profiler (ADCP), for profiling the 
sea current velocities; the C*, for measuring the water optical properties; the CTD, for 
salinity, temperature and density of water Bottom: Distribution of the reconstructed events 
as a function of the cosine of the zenith angle, cos{'d^^^), compared with the MUPAGEMC 



by using the Dense Wavelength Division Multiplex protocol. On shore, twin FCM boards 
de-multiplexed the incoming signal and distributed the data to the on-line trigger for a first 
raw selection of data. The trigger was based on coincidences occurred on near OMs within 
20 ns and on lai^ge amplitude single hits. When a trigger seed was found, all hits occurred 
within a time window of ±2 ^s centered on the seed time were recorded. 

A data analysis was done on a small sample of selected events, recorded during and 
24*^^ January 2007, corresponding to a livetime of 11.3 hours. From the analyzed data set, 
2260 atmospheric muon events were reconstructed and their angular distribution measured 
112451 (see Fig. [H bottom). 

The Phase- 1 project provided a fundamental test of the technologies proposed for the 
realization and installation of the detector. Some problems occurred in Phase- 1 after some 
months of functioning. Buoyancy of the tower decreased with the time (due to the construc- 
tion process of the buoy), producing a lowering of the tower position. Another problem was 
related to a malfunction inside the JB that required the recovery for a full diagnosis, which 
pointed out a malfunctioning of the optical penetrator. 

The Phase-2 was planned to validate the new solutions at the depths of the site of Capo 
Passero. In July 2007, a 100 km Alcatel electro-optical cable was laid on the seabed linking 
the 3500 m deep sea site to shore. The cable is a 10 kV DC, along a single electrical 
conductor, allowing a power transport larger than 50 kW. The DC/DC converter, which 
converts the high voltage coming from shore into the 400 V required for the detector, is 
produced by Alcatel and will be deployed by the end of 2009. The data transmission is 
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provided through 20 single mode optical fibers 112461 . The shore station, located inside the 
harbor area of Portopalo di Capo Passero, was completed during 2009 in an ancient winery 
building. 

A complete mechanical tower, a fully equipped NEMO mini-tower and a reduced ver- 
sion of an ANTARES string are planned to be installed on the Capo Passero Site by the 
end of 2009. Some features of the tower and mini-tower were modified according to the 
experience obtained from Phase- 1: (/) floor length is reduced from 15 m to 12 m; (//) higher 
simplification of the floor electric connections; (///) new time-calibration flashers embed- 
ded directly into the OMs; (iv) new electronic components which allow a lower power 
consumption; (v) a new acoustic positioning system with special road-band hydrophones 
able to measure the environmental acoustic background noise at 3500 m depth. 

The cable and the shore station were proposed to the KM3Net consortium as well suited 
to host the Mediterranean km^ detector. 

10.3 NESTOR 

The Neutrino Extended Submarine Telescope with Oceanographic Research 
(NESTOR) collaboration has developed an approach to operating a deepsea station, in the 
Southern Ionian Sea off the coast of Greece at depths exceeding 3500 m, permanently con- 
nected to shore by an insitu bidirectional cable, for multidisciplinary scientific research. 

The basic element of the proposed NESTOR detector is a 32 m diameter hexagonal 
floor (star). A central casing supports aim diameter spherical titanium pressure housing 
which contains the data acquisition electronics, power converters, monitoring, control and 
data transmission equipment. Attached to the central casing there are six arms built from 
titanium tubes to form a lightweight but rigid lattice girder structure. The arms can also be 
folded for transport and deployment. Two OMs are installed at the two end of each arm, 
one facing upwards and the other downwards: OMs are also installed above and below 
the central casing making a total of 14 units per floor. Using the OMs in pairs gives 47r 
coverage, enhancing discrimination between upward and downward going particles. In the 
NESTOR version, the tower will consist of 12 of such floors, spaced vertically 20-30m (see 
Fig. |37l left). 

In January 2002 a prototype was completed and deployed at a depth of 4100 m (project 
LAERTIS). The station transmitted the acquired data to shore from temperature and pres- 
sure sensors, compass, light attenuation meter, water current meter and an ocean bottom 
seismometer. After a period of some months, the station was recovered. 

In March 2003, the NESTOR collaboration successfully deployed a test floor of the 
detector tower, fully equipped with 12 OMs, final electronics and associated environmental 
sensors 112471 (see Fig. [37j right). The detector continuously operated for more than a 
month. For about 1.1% of the total experimental time, bioluminescent activity was observed 
around the detector. This caused about 1% dead time. The prolonged period of running 
under stable operating conditions made it possible to measure the cosmic ray muon flux as 
a function of zenith angle and to derive the deep intensity relation [|248 J . 



11 The KM3NeT Consortium 

KM3NeT is a future deep-sea research infrastructure planned to host a neutrino telescope 
with a volume of at least one cubic kilometre to be constructed in the Mediterranean Sea. In 
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Figure 37: Left: the proposed ftill NESTOR Detection Unit. Right: one floor of the NESTOR 
tower during its deployment for the first test in March 2003. 



February 2006, the Design Study for the infrastructure, funded by the EU FP6 framework, 
started. The KM3NeT research infrastructure has been singled out by ESFRI (the European 
Strategy Forum on Research Infrastructures) to be included in the European Roadmap for 
Research Infrastructures. The primary objective of the Design Study is the development of a 
cost-effective design for a cubic-kilometre sized deep-sea infrastructure housing a neutrino 
telescope with unprecedented neutrino flux sensitivity at TeV energies and providing long- 
term access for deep-sea research. In April 2008 the Conceptual Design Report for the 
KM3NeT infrastructure was made public [61]. 

The Preparatory Phase of the infrastructure, funded by the EU FP7 framework, started 
in March 2008. The primary objective of the KM3NeT Preparatory Phase is to pave the way 
to political and scientific convergence on legal, governance, financial engineering, aspects 
concerning the choice of the site and of the infrastructure and to prepare rapid and efficient 
construction once it gets approved. Reconciliation of national and regional political and 
financial priorities with scientific and technological considerations will be a major issue, 
as has become apparent in the KM3NeT Design Study. The construction of the KM3NeT 
infrastructure is foreseen to start after the three year Preparatory Phase, which has been 
organised in work packages. Each work package has its own coordinator and executive 
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Figure 38: The three sites candidate to host the KMSNet telescope. 



committee. 

Design, construction and operation of the KM3NeT neutrino telescope will be pursued 
by a consortium formed around the institutes currently involved in the ANTARES, NEMO 
and NESTOR pilot projects (see Fig. [38l ). Based on the leading expertise of these research 
groups, the development of the KM3NeT telescope is envisaged to be achieved, after the 
Preparatory Phase, within a period of about four years for construction and deployment. 

The KM3NeT facilities will provide support to scientific, long term and real-time mea- 
surements, also to a wide range of other associated earth and marine sciences, like oceanol- 
ogy, geophysics and marine biology. 

The KM3Net detector is expected to exceed IceCube in sensitivity by a substantial 
factor, exploiting the superior optical properties of seawater as compared to the Antarctic 
ice and an increased overall PMT area. In Table 1 , we summarize the required angular and 
energy resolution needed by the future KM3NeT detector, according to different types of 
astrophysical neutrino sources. 

The KM3NeT neutrino telescope will be composed of a number of vertical structures 
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Figure 39: Three geometry layouts for the KMSNeT detector: left, a squared grid; center, 
clustered Detection Units; right, ring configuration. 



(detection units) which are anchored to the sea bed and usually kept vertical by one or sev- 
eral buoys at their top. Since there is still a variety of viable design options, corresponding 
simulation studies are rather generic, concerning both assumed neutrino fluxes and detector 
properties. Fig. [39] shows three possible layouts: a squared grid of Detection Units (left), 
clustered (middle) and ring (right) configurations. Another configuration could consist of 
arranging the Detection Units in a homogeneous hexagon. 

Each detection unit will carry the photo-sensor and possibly further devices for cal- 
ibration and environmental measurements on mechanical structures which can be like 
ANTARES storeys, NEMO's y^oor^' or NESTOR's stars. Such structure will support the 
necessary sensors, supply interfaces, data lines and electronic components where applica- 
ble. The basic photo-sensor unit remains the Optical Module (OM) [249, 250], which can 
host one or several PMTs, their high-voltage bases and their interfaces to an acquisition 
system of nanosecond-precision. 

Whereas all of the current neutrino telescope projects use OMs composed of a single 
large (typically 10") standard PMT per OM, alternative solutions are also under investiga- 
tion for KM3NeT. In addition to the classical solution described above (see Fig. |40l case a), 
various tests with multi-cathode PMTs (see Fig. |40l case b), multi-PMT OMs (see Fig. |40l 
case c), and large spherical hybrid PMTs (see Fig. |40]case d) 11251 1 are performed together 
with computer simulation for studying the telescope response accordingly. 

The data transport devices and power harness of each Detection Unit is planned to be 
connected via the anchor to a deep-sea cable network. This network can contain one or more 
junction boxes and one or several electro-optical cables to shore. It also provides power and 
slow-control communication to the detector. On shore, a station equipped with appropriate 
computing power is required for collecting the data, applying online filter algorithms and 
transmitting the data to mass storage devices (see Fig. |4T]where the trigger and DAQ system 
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Table 1: Target sources, neutrino energy range, interaction channels and resolution con- 
straints for the KMSNeT telescope. 



is sketched). 

The deployment of the Detection Units on the sea bed and their maintenance along the 
years of the telescope live time require the development of appropriate machines and infras- 
tructures. For Detection Units deployment, the NESTOR Institute has developed a central- 
well, ballasted platform called Delta Berenike. For completing the detector construction 
(junctions between underwater connectors) and maintenance down in the deep-site, under- 
water remotely operated robotic devices will be necessary. 



12 Conclusion 

In this work, we have reviewed most of the phenomenological aspects of neutrino astro- 
physics, the status and future prospectives of neutrino Cherenkov detectors. 

Neutrino telescopes will cover the high-energy part of the neutrino spectrum from astro- 
physical or cosmogenic origin. The history is started [2531 with the detection of neutrinos 
{Ey '^1 — 10 MeV) from the Sun: it corresponds to the direct observation of the nuclear re- 
actions occurring in the core of stars 112541 12551 . Different experimental technologies were 
used to study solar neutrinos: radiochemical 112561125711258 1. water [259] and heavy- water 
112601 Cherenkov detectors and liquid scintillator 1126111 detectors. All experiments agree 
with the hypothesis of neutrino flavor oscillations. In addition to dedicated solar v ex- 
periments, other underground apparatus in the eighties were motivated by the proton decay 
prediction of most grand unified theories (GUTs); they were also MeV-GeV neutrino detec- 
tors [262 ]. Those in acquisition during 1987 were lucky enough to catch neutrinos from the 
supernova 1987 A in the Large Magellanic Cloud. The measurement of neutrinos emitted 
by SN1987A provided a direct observation of the core-collapse of a star 112631 . The size of 
these neutrino experiments (also the biggest ones, as Super-Kamiokande and MACRO) is 
too small to have the hope to discover high energy neutrinos from cosmic sources. With the 
advent of neutrino telescopes, only the very low energy relic neutrinos from cosmological 
origin seems to be far to be recognized. 

Some fundamental questions are still open after a century from the discovery of Cosmic 
Rays: which are the astrophysical engines in our Galaxy able to accelerate CRs up to ~ 10^ 
TeV and which are those in the Universe able to accelerate CRs up to more than 10^ TeV? 
lACT telescopes have found a variety of galactic and extragalactic TeV 7-ray sources whose 
spectrum extends up to some tens on TeV, where the intrinsic sensitivity of the present 
generation of instruments fall down. 
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Figure 40: The studied optical modules: (a) One single 10" PMTin a "benthos" sphere; (b) 
Multi-cathode PMT, with mirror separations to subdivide the PMT acceptance; (c) Multi- 
PMT optical module made by 20 PMTs, 3" each; (d) Spherical geometry X-HPD (8" pro- 
totype). 



At present, no definitive proofs exist that hadrons are also accelerated at the sites where 
TeV 7-ray are observed. If hadrons are present, a strong relationship between the CR energy 
spectrum and that of secondary 7-rays and neutrinos must exist. Compelling evidence 
for CR accelerators would be the detection of cosmic TeV neutrinos. In particular, no 
exponential cutoff is observed in the CR spectrum above ~ 10 TeV, despite what is reported 
by TeV 7-ray experiments for some galactic sources (see Fig. |3]l. The exponential cutoff 
would suppress the expected flux of neutrinos above 10 TeV: in this energy range, the 
detectors have a large neutrino effective area (Fig. [T6l ) and the v ^ jjL kinematic production 
allows a measurement of the neutrino direction within 0.2° in water (Fig. l24l ). With this 
angular resolution, the atmospheric neutrino background can be largely reduced. However, 
galactic CRs can reach 10^ TeV: the exponential cutoff observed in the TeV gamma ray 
spectrum of some sources can be due to an instrumental/lack of statistic effect, to gamma 
absorption effects at sources, or because the sources of CRs above few tens of TeV are 
hidden to gamma rays. In any case, engines that accelerate hadrons up to 10^ TeV exist in 
the Galaxy, and they are presumably accompanied by high energy neutrinos. 

Only a neutrino telescope in the Mediterranean Sea, with effective area > 25 times 
larger than that of the ANTARES detector, can probably solve this puzzle. The next years 
will thus be decisive for neutrino astronomy. While the IceCube experiment at the South 
Pole is in an advanced stage of construction, the technological challenges to build a huge 
neutrino telescope in deep sea have been surmounted by experiments like ANTARES, NE- 
MO and NESTOR, which have motivated the approval of the design study of the kilometer- 
scale version (the KM3NeT consortium). These new generation of neutrino telescope ex- 
periments will achieve effective volumes which will be able to explore the Northern sky (the 
IcuCube experiment in the South Pole) and the Southern sky (the underwater Mediterranean 
experiment) in a way never seen before. 

The advantage of the telescope in the Mediterranean Sea with respect to the IceCube 



72 



T. Chiarusi and M. Spurio 




Figure 41: The Trigger and Data Acquisition System (TRIDAS) for the KMSNeT project. 



detector is a better exposure to the galactic center region (Fig. |2l), and a much better angular 
resolution on the measurement of the neutrino direction . The latter has two consequences: 
the association possibility with known celestial objects; the reduction of the background 
from atmospheric neutrinos, which increases as the square of the search-window opening 
angle. On the other hand, the IceCube experiment has a better chance to observe the diffuse 
flux of neutrinos of extragalactic origin, being the theoretical upper limits accessible within 
few years (Fig. |7]l. Concerning the possibility that the detector in the Mediterranean Sea 
(if national and regional political and financial priorities cannot be easily conciliated with 
scientific and technological considerations) will be distributed in more than one site (the 
multisite option), this will not affect the sensitivity for the CC channel but reduces the 
confinement of the neutrino induced electromagnetic and hadronic showers and the possi- 
bility to detect with high efficiency the Vr channel. 

Identification of neutrino from cosmic accelerators has not been claimed so far; neutrino 
flux can only be evaluated using models and this sets the scale of the detector instrumented 
volume. The value of ~ 1 km-^ seems the minimum in order to observe neutrinos both 
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of galactic and extragalactic origin. The hunt for the first high energy neutrino of cosmic 
origin has started. 



Acknowledgments 

The attention of the high energy physics community is concentrated to the beginning of 
the LHC operations. One of the hidden motivation for this work is to present some of the 
most interesting open questions in astroparticle physics, particularly to young researchers 
and post-graduate students. For this reason, we primarily acknowledge many young col- 
leagues and PhD students, mainly from the ANTARES and NEMO collaboration and from 
the KM3NeT consortium, whose analyses and figures where used and quoted in this work. 
The authors wishes to thank many colleagues for useful discussion about the content of 
this paper from all the "water and ice" neutrino community and in particular our colleagues 
in Bologna (M. Bazzotti, S. Biagi, G. Carminati, S. Cecchini, G. Giacomelli and A. Mar- 
giotta). We thank E. Guzman De La Hoz for reading the final version. Finally, the authors 
are grateful for useful comments and suggestions from the anonymous referees. This work 
was partially supported (T.C.) by a grant of the Preparatory Phase of KM3NeT infrastruc- 
ture, funded by the EU FP7 framework. 

References 

[1] M.A. Markov, Proceedings Int. Conf. on High Energy Physics, p. 183, Univ. of 
Rochester (1960). 

[2] F. Halzen. Astroparticle physics with high energy neutrinos: from AMANDA to Ice- 
Cube. Eur. Phys. J. C 46 (2006) 669-687 

[3] J. Babson et al. (The DUMAND Coll.). Cosmic-ray muons in the deep ocean. Phys. 
Rev. D42 (1990) 3613. 

[4] V. Aynutdinov et al. (Baikal Collaboration). The Baikal neutrino experiment: Status, 
selected physics results, and perspectives. Nucl. Instrum. Meth. A588 (2008) 99-106. 

[5] E. Andres et al. (the AMANDA Collaboration). Observation of high-energy neutrinos 
using Cherenkov detectors embedded deep in Antarctic ice. Nature 410,(2001) 441. 

[6] J. A. Aguilar et al. (ANTARES Collaboration). Astrop. Phys. 26 (2006) 314. 

[7] P. Piattelli (NEMO collaboration). The Neutrino Mediterranean Observatory project. 
NucLPhys. Proc. Suppl. 143 (2005) 359. 

[8] G. Aggouras et al. (NESTOR collaboration). A measurement of the cosmic-ray muon 
flux with a module of the NESTOR neutrino telescope. Astrop. Phys. 23 (2005) 377- 
392. 

[9] http://www.km3net.org 

[10] A. M. Hillas Recent progress and some current questions, astro-ph/0607109. The fig- 
ure is due to T. Gaisser. 



74 



T. Chiarusi and M. Spurio 



[11] P.L. Biermann, T. K. Gaisser , T Stanev. Origin of galactic cosmic rays. Phys. Rev. D 
51 (1995) 3450. 

[12] T. K. Gsdsser.Cosmic Rays and Particle Physics. Cambridge University Press, 1990. 
ISBN 0521326672. 

[13] P. K. F. Grieder. Cosmic Rays at Earth: Researcher's Reference Manual and Data 
Book. Elsevier, 2001. ISBN 0444507108. 

[14] E. Fermi. On the origin of cosmic rays. Phys. Rev. 75(1949) 1169. 

[15] E. Fermi. Galactic magnetic fields and the origin of cosmic radiation. Astrophys. J. 
119 (1954) 1-6. 

[16] V. L. Ginzburg and V. S. Ptuskin. On the origin of cosmic rays: Some problems in 
high-energy astrophysics. Rev. Mod. Phys. 48 (1976) 161. 

[17] A. M. Hillas. Can diffusive shock acceleration in supernova remnants account for 
high-energy galactic cosmic rays? J. Phys. G: Nucl. Part. Phys. 31 (2005) 39. 

[18] E. G. Berezliko and L. T. Ksenofontov. Composition of cosmic rays accelerated in 
supernova remnants. JETP 89 (1999) 391. 

[19] V.L. Ginzburg, Y.M. Khazan and V.S. Ptuskin. Origin of cosmic rays - Galactic 
models with halo. I - Proton nucleon component. Astrophysics and Space Science 
68 (1980) 295-314. 

[20] J. R. Horandel. Models of the knee in the energy spectrum of cosmic rays. Astrop. 
Phys. 21 (2004) 241-265. 

[21] D. Kazanas and A. Nicolaidis. Cosmic Ray Knee: A Herald of New Physics?. In Pro- 
ceedings of the 27th International Cosmic Ray Conference, Hamburg, Germany, 2001 . 
astro-ph/0103147. 

[22] T. Antoni et al. (KASCADE Coll.) Electron, muon, and hadron lateral distributions 
measured in air-showers by the KASCADE experiment. Astrop. Phys. 14 (2001) 245- 
260. 

[23] C. E. Fichtel and J. I. Trombka. Gamma-ray astrophysics: New insight into the Uni- 
verse. Technical report, 1997. NASA-RP-1386. 

[24] N. Gehrels and P. Michelson. GLAST: The next-generation high energy gamma ray 
astronomy mission. Astrop. Phys. 11 (1999) 277- 282. 

[25] A. Abdo et al. (Fermi LAT Coll.) Fermi Large Area Telescope Bright gamma-rays 
source list. Astrophys. J. Suppl.183 (2009) 46-66. arXiv:0902.1340 [astro-ph.HE] 

[26] T. C. Weekes et al. Observation of TeV gamma rays from the Crab nebula using the 
atmospheric Cherenkov imaging technique. Astrophys. J. 342 (1989) 379-395. 

[27] R. Mirzoyan et al. The first telescope of the HEGRA air Cherenkov imaging telescope 
array. Nucl. Instr. and Meth. Phys. Res. A 351 (1994) 513-526. 



High-Energy Astrophysics with Neutrino Telescopes 



75 



[28] J. Holder et al. Status and Performance of the First VERITAS telescope. Proceedings of 
the 29th International Cosmic Ray Conference, Pune, India, 2005. astro-ph/050745 1 . 

[29] R. Enomoto et al. Design study of CANGAROO-III, stereoscopic imaging atmospheric 
Cherenkov telescopes for sub-TeV gamma-ray. Astrop. Phys. 16 (2002) 235-244. 

[30] J. A. Hinton. The status of the HESS project. New Astron. Rev. 48 (2004) 331-337. 

[31] T.C. Weekes et al. VERITAS: the Very Energetic Radiation Imaging Telescope Array 
System Astropart. Phys. 17 (2002) 221-243. 

[32] J. Albert et al. Physics and Astrophysics with a ground-based gamma-ray telescope of 
low energy threshold. Astroparticle Physics 23 (2005) 493-509. 

[33] A. de Angelis et al. Very-high-energy gamma astrophysics. Riv. Nuovo Cimento 31 
(2008) 187. 

[34] F. Aharonian et al. High energy astrophysics with ground-based gamma ray detectors. 
Rep. Prog. Phys. 71 (2008) 096901. 

[35] F. Halzen and D. Hooper, High-energy neutrino astronomy: The Cosmic ray connec- 
tion. Rept. Prog. Phys. 65 (2002) 1025-1078. 

[36] F. A. Aharonian. Very High Energy Cosmic Gamma-Ray Radiation. World Scientific 
PubUshing, 2004. 

[37] V. Berezinsky et al. Astrophysics of Cosmic Rays. North-Holland, Amsterdam (1990). 

[38] C. Stegmann, A. Kappes, J. Hinton and F. Aharonian. Potential neutrino signals in 
a northern hemisphere neutrino telescope from galactic gamma-ray sources. Astro- 
physics and Space Science 309 (2007) 429. 

[39] F. Aharonian et al. Detection of TeV gamma-ray emission from the shell-type super- 
nova remnant RX J0852.0-4622 with HESS. Astron. & Astrophys. L7 (2005) 437. 
astro-ph/0505380. 

[40] M. D. Kistler and J. F. Beacom. Guaranteed and prospective galactic TeV neutrino 
sources. Phys. Rev. D 74 (2006) 063007. astro-ph/0607082. 

[41] O. Reimer and M. Pohl. No hadronic TeV gamma-rays from SNR RX JI7I3.7-3946. 
Astron. & Astrophys. L43 (2002) 390. 

[42] H. Muraishi et al. Evidence for TeV gamma-ray emission from the shell-type SNR RX 
JI7I3.7-3946. Astron. & Astrophys. L57 (2000) 354. 

[43] R. Enomoto et al. The acceleration of cosmic-ray protons in the supernova remnant 
RXJI71 3.7-3946. Nature, 416 (2002)823. Astron. & Astrophys. 449 (2006) 223. 

[44] J. Alvarez-Muniz and F. Halzen. Possible high-energy neutrinos from the cosmic ac- 
celerator RXJ 17 13. 7 -3946. Astrophys. J. L33 (2002) 576. astro-ph/0205408. 

[45] M. L. Costantini and F. Vissani. Expected neutrino signal from supernova remnant 
RXJ 17 13. 7 -3946 and flavor oscillations. Astrop. Phys. 23 (2005) 477-485. 



76 



T. Chiarusi and M. Spurio 



[46] F. Aharonian et al. (HESS Collaboration). A detailed spectral and morphological study 
of the gamma-ray supernova remnant RX J 17 13. 7 -3946 with HESS. As\ion.8c Asti"o- 
phys. 449 (2006) 223-242. 

[47] D. Homs, F. Aharonian, A. Santangelo, A. I. D. Hoffmann, and C. Masterson. Nucle- 
onic gamma-ray production in VelaX. Astron. & Astrophys., 451 (2006) 51. 

[48] D. Guetta and E. Amato. Neutrino flux predictions for galactic plerions. Astrop. Phys. 
19 (2003) 403. asti-o-ph/0209537. 

[49] F. Aharonian et al. (HESS Collaboration). Very high energy gamma rays from the 
direction of Sagittarius A*. Asti^on. & Astrophys., 425:L13-17, 2004. 

[50] F. Aharonian et al. (HESS Collaboration). Very high energy gamma rays from the 
composite SNR GO.9+0.1. Astron. & Astrophys. L25-29 (2005)432. 

[51] F. Aharonian et al. (HESS Collaboration). Discovery of very high energy gamma-rays 
from the galactic centre ridge. Nature, 439 (2006) 695. astro-ph/0603021. 

[52] S. Chaty. Microquasars and Jets. Proceedings of Rencontres de Moriond, Very High 
Energy Phenomena in the Universe, La Thuile, Italy, 2005. astro-ph/0506008. 

[53] S. Migliari, R. Fender, and M. Mendez. Iron emission lines from extended X-ray jets 
in SS 433: Reheating of atomic nuclei. Science, 297 (2002) 1673. 

[54] J. Albert et al. (MAGIC Collaboration). Variable very high energy gamma-ray emis- 
sion from the microquasar LS I-i-61 303. Science, 312 (2006) 1771. astro-ph/0605549. 

[55] F. Aharonian et al. (HESS Collaboration). Discovery of very high energy gamma rays 
associated with an X-ray binary. Science, 309 (2005) 746. 

[56] H. Anderhub et al. (MAGIC Collaboration). Correlated X-ray and Very High Energy 
emission in the gamma-ray binary LS 1 +61 303. arXiv:0910.4381vl 

[57] D. F. Torres and F. Halzen. LS 1+61 303 as a potential neutrino source on the light of 
MAGIC results, asti-o-ph/0607368. 

[58] F. Aharonian, L. Anchordoqui, D. Khangulyan, and T. Montaruli. Microquasar LS 
5039: a TeV gamma-ray emitter and a potential TeV neutrino source. J. Phys. Conf. 
Ser., 39:408. 

[59] S. Aiello et al. (The NEMO Collaboration), Sensitivity of an underwater Cherenkov 
km^ telescope to TeV neutrinos from Galactic microquasars. Astrop. Phys. 28(2007) 
1-9. 

[60] R. Atkins et al. (MILAGRO collaboration). Evidence for TeV Gamma-Ray Emission 
from a Region of the Galactic Plane. Phys. Rev. Lett. 95 (2005) 251 103. 

[61] KM3neT: Conceptual Design Report. http://www.km3net.org/CDR/CDR- 
KM3NeT.pdf 

[62] W. Bednarek, G. F. Burgio, and T. Montaruli. Galactic discrete sources of high energy 
neutrinos. New Astron. Rev., 49 (2005) 1. astro-ph/0404534. 



High-Energy Astrophysics with Neutrino Telescopes 



77 



[63] G. D'Ali Staiti, The ARGO-YBJ experiment in Tibet. Nucl. Instr. Meth. A588 (2008) 
7. See also: http://argo.na.infn.it/ 

[64] MILAGRO project, http://www.lanl.gov/milagro 

[65] Slate: http://www.slate.com/id/2081042/ 

[66] F. Halzen et al. The highest energy cosmic ray. Astrop. Phys., 3 (1995) 151. 

[67] B. Peters. Primary cosmic radiation and extensive air showers. Nuovo Cimento, 22 
(1961) 800. 

[68] P. Sokolsky. Introduction to Ultrahigh Energy Cosmic Ray Physics, Westview Press 
Advanced Book Program, 2004. 

[69] M. Ostrowski. Mechanisms and sites of ultra high energy cosmic ray origin, Astrop. 
Phys. 18 (2002) 229-236. 

[70] K. Greisen. End to the cosmic ray spectrum, Phys. Rev. Lett. 16 (1966) 748. 

[71] G. T. Zatsepin and V. A. Kuzmin. Upper limit of the spectrum of cosmic rays. Sov. 
Phys. JETP Lett. 4 (1966) 78. 

[72] M. Takeda et al.. Extension of the cosmic ray energy spectrum beyond the predicted 
Greisen-Zatsepin-Kuz'min cutoff, Phys. Rev. 81 (1998) 1163. 

[73] T. Abu-Zayyad et al. (HiRes Collaboration). Measurement of the flux of ultrahigh 
energy cosmic rays from monocular observations by the high resolution Fly 's Eye 
experiment. Phys. Rev. Lett. 92 (2004) 151101. 

[74] See the web site of the Pierre Auger observatory: http://www.auger.org/ 

[75] The Pierre Auger Collaboration: Upper limit on the cosmic-ray photon flux above 
10^^ eV using the surface detector of the Pierre Auger Observatory Astrop. Phys. 29 
(2008) 243-256, arXiv:0712.1147. 

[76] J. Abraham et al. Correlation of the highest-energy cosmic rays with nearby extra- 
galactic objects. Science, 318 (2007) 938-943. 

[77] A. Achterberg et al. (IceCube Collaboration and PL. Biermann). On the selection of 
AGN neutrino source candidates for a source stacking analysis with Neutrino Tele- 
scopes. Astrop. Phys. 26 (2006) 282-300. 

[78] R. Engel, D. Seckel, and T. Stanev. Neutrinos from propagation of ultra-high energy 
protons. Phys. Rev. D, 64 (2001) 093010. astro-ph/0101216, and references therein. 

[79] E. Waxman. Extra galactic sources of high energy neutrinos. Phys. Scripta, 
T121(2005) 147. astro-py0502159. 

[80] R. V. E. Lovelace. Dynamo model of double radio sources. Nature, 262 (1976) 649. 

[81] K. Mannheim. High-energy neutrinos from extragalactic jets. Astrop. Phys., 3 (1995) 
295. 



78 



T. Chiarusi and M. Spurio 



[82] F. Halzen and E. Zas. Neutrino fluxes from active galaxies: A model independent 
estimate. Astrophys. J., 488 (1997) 669. 

[83] R. J. Protheroe. High energy neutrinos from blazars. ASP Conf. Sen, 121 (1997) 585. 

[84] J. K. Becker, P. L. Biermann, and W. Rhode. The diffuse neutrino flux from FR-II radio 
galaxies and blazars: A source property based estimate. Astrop. Phys. 23 (2005) 355- 
368. astro-ph/0502089. 

[85] F. Halzen, A. O'Murchadha. Neutrinos from Auger sources. arXiv 0802.0887 [astro- 
ph]. 

[86] R.C. Hartman et al. (EGRET Coll.) The Third EGRET Catalog of High-Energy 
Gamma-Ray Sources. The Astrophysical J. Suppl. Series 123 (1999) 79-202. 
ftp://gamma.gsfc.nasa.gov/pub/THlRD_CATALOG/ 

[87] F. Aharonian et al. (H.E.S.S. collaboration). A low level of extragalactic background 
light as revealed by big gamma-rays from blazars. Nature 440 (2006) 1018-1021. 

[88] W. Benbow et al. (H.E.S.S. collaboration). A spectacular VHE gamma-ray out- 
burst from PKS 2155-304 in 2006. Proc. 30th International Cosmic Ray Conference, 
Merida, 2007. 

[89] T. J. Galama et al. Discovery of the peculiar supernova 1998bw in the error box of 
GRB980425. Nature, 395 (1998) 670. 

[90] K. Z. Stanek et al. Spectroscopic discovery of the supernova 2003dh associated with 
GRB030329. Astrophys. J. L17 (2003) 591. 

[91] J. Hjorth et al. A very energetic supernova associated with the gamma-ray burst of 29 
March 2003. Nature, 423 (2003) 847. 

[92] http://heasarc.gsfc.nasa.gov/docs/sax/sax.html 

[93] E. Waxman and J. N. Bahcall. High energy neutrinos from cosmological gamma-ray 
burst flreballs. Phys. Rev. Lett. 78 (1997) 2292. 

[94] K. Murase, High energy neutrino early afterglows from gamma-ray bursts revisited. 
Phys. Rev. D76 (2007) 123001. 

[95] http://heasarc.gsfc.nasa.gov/docs/swift/swiftsc.html 

[96] D. Guetta et al. Neutrinos from individual gamma-ray bursts in the BATSE catalog. 
Astrop. Phys., 20 (2004) 429. astro-ph/0302524. 

[97] M. Bouwhuis. Detection of neutrinos from gamma ray bursts. PhD the- 
sis July 2005. Amsterdam university- Amsterdam, The Netherlands. 
http://antares.in2p3.fr/Publications/ 

[98] J.A. Aguilar et al. (The ANTARES Coll.). The data acquisition system for the 
ANTARES neutrino telescope. Nucl. Instr. Meth. A570 (2007) 107 . 

[99] M. Bouwhuis (ANTARES Coll.) Search for gamma-ray bursts with the ANTARES 
neutrino telescope, Proceedings of the 31st ICRC, Lodz 2009. arXiv: 0908.0818. 



High-Energy Astrophysics with Neutrino Telescopes 



79 



[100] H. Dannerbauer et al. Follow-up near-infrared spectroscopy of ultraluminous in- 
frared galaxies observed by ISO. Astron. & Astrophys., 441 (2005) 999. 

[101] V. Springel et al. Modeling fredback from stars and black holes in galaxy mergers. 
Astron. Soc, 361 (2005) 776. 

[102] F. Acero et al. (H.E.S.S. Coll.) Detection of Gamma Rays from a Starburst Galaxy. 
Science Express, Sept. 24, 2009 DOl: 10. 1126/science. 1178826. arXiv:0909.4651 

[103] N. Karlsson (VERITAS Coll.). VERITAS Detection of Starburst Galaxy M 82. 2009 
Fermi Symposium. 

[104] A. Loeb and E. Waxman. The cumulative background of high-energy neutrinos from 
starburst galaxies. JCAP, 0605 (2006) 003. astro-ph/0601695. 

[105] P. Sreekumar et al. (EGRET Collaboration). EGRET Observations of the Extra- 
galactic Gamma-Ray Emission. Ap. J. 494 (2000) 523-534. 

[106] E. Waxman and J. Bahcall. High energy neutrinos from astrophysical sources: An 
upper bound. Phys. Rev. D59 (1998) 023002-1; J. Bahcall, E. Waxman. High energy 
astrophysical neutrinos: The upper bound is robust. Phys. Rev. D64(2001) 023002-1 

[107] K. Mannheim, R. J. Protheroe, J. P. Rachen. Cosmic ray bound for models of extra- 
galactic neutrino production, Phys. Rev. D 63 (2000) 023003. 

[108] J. Zornoza. Sensitivity to diffuse fluxes and energy spectrum reconstruction in the 
Antares neutrino telescope. PhD thesis January 2005. Universitat de Valncia- Valncia, 
Spain. http://antares.in2p3 .fr/PubUcations/ 

[109] W. Rhode et al. (Frejus Collaboration). Limits on the flux of very high energy neutri- 
nos with the Frjus detector. Astrop. Phys. 4 (1996) 217 . 

[110] M. Ambrosio et al. (MACRO Coll.) Search for diffuse neutrino flux from astrophysi- 
cal sources with MACRO. Astrop. Phys. 19 (2003) 1. 

[Ill] V. Aynutdinov et al. Search for a Diffuse Flux of High-Energy Extraterrestrial Neu- 
trinos with the NT200 Neutrino Telescope. Astrop. Phys. 25 (2006) 140-150. astro- 
ph/0508675. 

[112] A. Achterberg et al. (IceCube Coll.). Multiyear search for a diffuse flux of muon 
neutrinos with AMANDA II. Physical Review D, 76 (2007) 042008. 

[113] M. Ackermann et al. (AMANDA Collaboration). Search for Ultra-High-Energy Neu- 
trinos withAMANDA-II. Ap. J- 675 (2008) 1014. 

[114] G. Fiorentini, A. Naumov, F.L. Villante. Atmospheric neutrino flux supported by 
recent muon experiments. Phys. Lett. B510 (2001) 173-188. 

[115] G. t'Hooft. Magnetic monopoles in unified gauge theories. Nucl. Phys. B 79 (1974) 
276. 

[116] A.M. Polyakov, Sov. Phys. JETPLett. 20, 194(1974). 



80 



T. Chiarusi and M. Spurio 



[117] E.N. Parker. The Origin of Magnetic Fields. Ap. J. 160(1970) 383. M.S. Turner 
et al. Magnetic monopoles and the survival of galactic magnetic fields. Phys. Rev. 
D26(1982)1296. 

[118] M. Ambrosio et al. (MACRO Coll.). Final results of magnetic monopole searches 
with the MACRO experiment. Eur. Phys. Jour. C25 (2002) 511-522. 

[119] B. van Rens. Detection of MM below the Cherenkov limit Phd. D. thesis, Amsterdam 
U., 2006. See: http://antares.in2p3.fr/Publications/ 

[120] J. Derkaoui et al. Energy losses of magnetic monopoles and dyons in the earth. 
Astropart. Phys. 9 (1998), 173. 

[121] S. D. Wick et al. Signatures for a cosmic flux of magnetic monopoles. Astropart. 
Phys. 18 (2003) 663-687. 

[122] K. Antipin et al. (BAIKAL Coll.). The Baikal Neutrino Telescope: Selected Physics 
Results. arXiv:07 10.3064V 1 [astro-ph] 

[123] P. Niessen. Search for relativistic MM with the AMANDA detector. Phd. D. thesis, 
Humboldt U., Berlin, 2001. 

[124] S. Orito et al. Search for supermassive relics with a 2000 m^ array of plastic track 
detectors. Phys. Rev. Lett. 66 (1991) 1951. 

[125] A. De Rujula, S.L. Glashow. Nuclearitesa novel form of cosmic radiation. Nature 
312(1984)734 

[126] M. Ambrosio et al. (The MACRO coll.). Nuclearite search with the MACRO detector 
at gran Sasso. Eur. Phys. J. C13 (2000) 453-458. 

[ 1 27] G. Pavalas (ANTARES Coll.) Search for Exotic Physics with the ANTARES Detector. 
arXiv:0908.0860 [astro-ph.IM] 

[128] G. Bertone, D. Hooper, J. Silk. Particle dark matter: evidence, candidates and con- 
straints. Phys. Rep. 405 (2005) 279-390. 

[129] H.E. Haber and G.L. Kane. The search for supersymmetry: Probing physics beyond 
the standard model. Phys. Rep. 117 (1985) 75. 

[130] N. Fomengo. Status and perspectives of indirect and direct dark matter searches. 
Adv. Space Res. 41 (2008) 2010-2018. arXiv: astro-ph/0612786vl 

[131] S. Desai et al. (SK Coll.). Search for dark matter WIMPs using upward through- 
going muons in Super-Kamiokande. Phys. Rev. D 70 (2004) 083523. 

[132] M. Ambrosio et al. (MACRO Coll.) Limits on dark matter wimps using upward 
going muons in the MACRO detector. Phys. Rev. D 60 (1999) 082002. 

[133] D. Hubert (IceCube Coll.). Neutralino Dark Matter Searches with Neutrino Tele- 
scopes: AMANDA Results and IceCube Prospects. Nucl. Phys. B (Proc. Suppl.) 173 
(2007) 87-90. 



High-Energy Astrophysics with Neutrino Telescopes 



81 



134] V. M. Aynutdinov et al. (Baikal Coll.). The Baikal Neutrino Experiment: Status and 
beyond. Physics of Atomic Nuclei 67, Number 6 (2004) 1161-1171. 

135] G.M.A Lim (ANTARES Coll.). First results on the search for dark matter in the Sun 
with the ANTARES neutrino telescope. arXiv:0905.2316, 2009. ICRC 2009. 

136] C. Giunti. Neutrino Flavor States and the Quantum Theory of Neutrino Oscillations. 
J. Phys. G: Nucl. Part. Phys. 34 (2007) R93-R109, arXiv:hep-ph/0608070. See also 
the very interesting and exhaustive web page on neutrinos: http://www.nu.to.infn.it/ 

137] W-M Yao et al. (Particle Data Group). Review of Panicle Physics. J. Phys. G 33, 1 
(2006) and 2007 partial update for the 2008 edition. 

138] J. G. Learned and S. Pakvasa. Detecting nu-tau oscillations as PeV energies. Astrop. 
Phys., 3 (1995) 267. 

139] G.L. Fogli et al. Testing violations of special and general relativity through the en- 
ergy dependence of — Uj- oscillations in the Super-Kamiokande atmospheric neu- 
trino experiment. Phys. Rev. D 60 (1999) 053006; M.C. Gonzalez-Garcia, M. Maltoni. 
Atmospheric neutrino oscillations and new physics. Phys. Rev. D 70 (2004) 033010. 

140] D. Morgan et al. Neutrino telescope modelling of Lorentz invariance violation in 
oscillations of atmospheric neutrinos. Astropart. Phys. 29 (2008) 345-354. 

141] O.W. Greenberg. CPT Violation Implies Violation of Lorentz Invariance. Phys. Rev. 
Lett. 89 (2002) 231602. 

142] G. Battistoni et al. Search for a Lorentz invariance violation contribution in atmo- 
spheric neutrino oscillations using MACRO data. Phys. Lett. B615 (2005) 14. 

143] J. L. Kelley (IceCube Coll.) Testing Lorentz Invariance using Atmospheric Neutrinos 
andAMANDA-II. astro-ph/0701333 

144] M. Gasperini. Testing the principle of equivalence with neutrino oscillations. Phys. 
Rev. D 38 (1988) 2635. 

145] V. Van Elewyck (ANTARES Coll.). Searching for high-energy neutrinos in coinci- 
dence with gravitational waves with the ANTARES and VIRGO/LIGO detectors, e- 
Print: arXiv: 0908. 2454 [astro-ph.IM] 

146] G. Amehno-CameUa, Int. J. Mod. Phys. D12 (2003) 1633; S. Choubey and S.F. King, 
Phys. Rev. D67 (2003) 073005; U. Jacob and T. Pkan, Nature Phys. 3 (2007) 87. 

147] F. Acernese (VIRGO Collaboration). Virgo status. Class. Quant. Grav. 25 (2008) 
184001. See also: http://www.virgo.infn.it/. 

148] D. Sigg (LIGO Collaboration). Status of the LIGO detectors. Class. Quant. Grav. 25 
(2008) 114041. See also: http://www.Hgo.caltech.edu/. 

149] P. Sapienza and G. Riccobene. High-energy neutrino astronomy. Rivista del Nuovo 
Cimento v.32, n. 12 (2009). 



150] R. Gandhi et al. Ultrahigh-energy neutrino interactions. Astrop. Phys. 5 (1996) 81. 



82 



T. Chiarusi and M. Spurio 



[151] J. Pumplin et al. New generation of parton distributions with uncertainties from 
global QCD analysis. JHEP, 07 (2002) 012. hep-ph/0201 195. 

[152] See for instance: http://durpdg.dur.ac.uk/hepdata/pdf.html 

[153] M. Ambrosio et al. Neutrino astronomy with the MACRO detector. Astrophys. J., 546 
(2001) 1038-1054. 

[154] K. Abe et al. High energy neutrino astronomy using upward-going muons in Super- 
Kamiokande-I. Astrophys. J., 652 (2006) 198-205. 

[155] M. Bazzotti. Studies of the atmospheric muon flux with the ANTARES detec- 
tor. PhD thesis Mars 2009. Universit degh studi di Bologna, Bologna, Italy. 
http://antares.in2p3.fr/Publications/ 

[156] Y. Becherini, A. Margiotta, M. Sioli and M. Spurio. A parameterisation of single and 
multiple muons in the deep water or ice. Astrop. Phys. 25 (2006) 1-13. 

[157] V. Agrawal, T.K. Gaisser, P. Lipari, T. Stanev. Atmospheric neutrino flux above 1 
GeV. Phys. Rev. D53 (1996) 1314. 

[158] F. Halzen, D. Saltzberg. Tau Neutrino Appearance with a 1000 Megaparsec Base- 
line. Phys. Rev. Lett. 81 (1998) 4305. 

[159] B. D. Hartmann. Reconstruction of neutrino-induced hadronic and electromagnetic 
showers with the ANTARES experiment. PhD tesi, Erlangen U. (2006). FAU-PIl- 
DISS-06-001 . http://antares.in2p3.fr/ 

[160] S.L. Glashow. Resonant scattering ofantineutrinos. Phys. Rev., 118 (1960) 316. 

[161] L. Landau and I. Pomeranchuk, Dokl. Akad. Nauk SSSR 92, 735 (1953); [Enghsh 
translation: L. Landau, The Collected Papers of L.D. Landau (Pergamon, New York, 

1965), p. 589]. 

[162] A.B. Migdal. Bremsstrahlung and pair production in condensed media at high en- 
ergies. Phys. Rev. 103 (1956) (6). 

[163] J. Brunner. Antares simulation tools. VLVnT Workshop on Technical Aspects of a 
Very Large Volume Neutrino Telescope in the Mediterrenean Sea, Amsterdam, The 
Netherlands, 5-8 Oct 2003. http://www.vlvnt.nl/proceedings/ 

[164] D.E. Groom et al. Muon stopping power and Range Tables 10 MeV-100 TeV. Atomic 
Data and Nuclear Data Tables 78, vol. 2 (2001) 183-356. 

[165] P. Lipari and T. Stanev. Propagation of multi-TeV muons. Phys. Rev. D 44 (1991) 
3543. 

[166] P. A. Cherenkov. Visible Radiation Produced by Electrons Moving in a Medium with 
Velocities Exceeding that of Light. Phys. Rev. 52 (1937) 378. 

[167] F. Aharonian et al. (HEGRA Collaboration). The Crab nebula and pulsar between 
500 GeV and 80 TeV. Observations with the HEGRA Stereoscopic Air Cherenkov Tele- 
scopes. Astrophys. J. 614 (2004) 897. 



High-Energy Astrophysics with Neutrino Telescopes 



83 



[168] A. Achterberg, Y.A. Gallant, J.G. Kirk, A.W.Guthmann. Particle acceleration by 
ultra-relativistic shocks: theory and simulations. Mon. Not. Roy. Astron. Soc. 328 
(2001) 393. astro-ph/0107530. 

[169] T.C. Weekes. et al. Observation of TeV gamma rays from the Crab nebula using the 
atmospheric Cerenkov imaging technique. Astrophys. J. 342 (1989) 379. 

[170] F. Ahronian et al. (The HESS Coll.), Observations of the Crab Nebula with H.E.S.S., 
Astron. Astrophys. 457 (2006) 899-915. 

[171] T. K. Gaisser, F. Halzen and T. Stanev. Particle Astrophysics with High-Energy Neu- 
trinos. Phys. Rep. 258 (1995) 173-236; Erratum-ibid. 271 (1996) 355-356, arXiv:hep- 
ph/9410384vl. 

[172] D. J. L. Bailey. Monte Carlo tools and analysis methods of understanding the 
ANTARES experiment and predicting its sensitivity to Dark Matter. PhD thesis. Uni- 
versity of Oxford, United Kingdom, 2002. 

[173] M. Ackermann et al. Optical properties of deep glacial ice at the South Pole. J. 
Geophys. Res. 111(D13) (2006) 1-4. 

[174] Y. R. Whitehead. On the Properties of Ice at the IceCube Neutrino Telescope. Mas- 
ter of Science thesis. University of Canterbury, Christchurch, New Zealand, 2008. 
http://www.icecube.wisc.edu/science/ 

[175] Y. He and P. Price. Remote sensing of dust in deep ice at the South Pole. Journal of 
Geophysical Research, 103(D14) (1998) 17041-17056. 

[176] G. Riccobene et al. Deep seawater inherent optical properties in the southern Ionian 
sea. Astroparticle Physics 27 (2007) 1-9. 

[177] R.C. Smith, K.S. Baker. Optical properties of the clearest natural waters (200800 
nm). Appl. Opt. 20 (1981) 177. 

[178] http://www.wetlabs.com 

[179] CD. Mobley, Light and Water. Academic Press, San Diego, 1994. 

[180] M. Circella, The construction of ANTARES, the first undersea neutrino telescope. 
Nucl. Instrum Meth. A602 (2009) 1-6. 

[181] R.M. Pope, E.S. Fry. Absorption spectrum (380700 nm) of pure water II. Integrating 
cavity measurements. Appl. Opt. 36 (1997) 8710-8723. 

[182] L. Kou, D. Labriel, P. Chylek. Refractive indices of water and ice in the 0.65 to 2.5 
pm spectral range. Appl. Opt. 32 (1993) 3531. 

[183] P. J. Henning. Bioluminescence in action. Academic Press, London, 1978. 

[184] E.A. Widder, l.M. Latz and J.F.Case. Marine bioluminescence spectra measured with 
an optical multichannel detection system. Biol. Bull. 165 (1983) 791. 

[185] E.J.V. Gillibrand et al. Seasonal development of a deep pelagic bioluminescent layer 
in the temperate NE Atlantic Ocean. Marine Ecology Progress, Series 341 (2007) 37. 



84 



T. Chiarusi and M. Spurio 



[186] I.G. Priede et al. Bioluminescence in the deep sea: free-fall lander observations in 
Atlantic Ocean off Cape Verde. Deep Sea Research I, 53 (2006) 1272. 

[187] J. Lundberg et al. Light tracking through ice and water - Scattering and absorption 
in heterogenous media with PHOTONICS. Nucl. Instrum. and Meth. A581 (2007) 
619-631. 

[188] J. Dumm, H. Landsman for the IceCube Collaboration: IceCube - First Results. 
Journal of Physics: Conference Series 60 (2007) 334-336. 

[189] J. Ahrens et al. (the AMANDA Collaboration). Muon track reconstruction and data 
selection techniques in AMANDA, Nucl. Instrum. Meth. A524 (2004) 169-194. 

[190] A. Heijboer. Track reconstruction and point source searches with Antares. 
PhD Thesis (2004) Universiteit van Amsterdam, Amsterdam, The Netherlands. 
http://antares.in2p3.fr/Publications/ 

[191] G. Carminati, M. Bazzotti, A. Margiotta and M. Spurio. Atmospheric MUons from 
PArametric formulas: a fast GEneratorfor neutrino telescopes (MUPAGE). Comput. 
Phys. Commun. 179 (2008) 915. arXiv:0802.0562[physics.ins-det] 

[192] G. Carminati, A. Margiotta, M. Spurio. MUons from Parametric formulas: A fast 
Generator of atmospheric fx-bundles for neutrino telescopes( MUPAGE). Nucl. Instr. 
Methods A602 (2009) 95. 

[193] http://www.phys.hawaii.edu/dmnd/dumand.html 

[194] A. Roberts. The birth of high-energy neutrino astronomy: A personal history of the 
DUMAND project. Rev. Mod. Phys. 64 (1992) 259. 

[195] http://baikalweb.jinr.ru/ 

[196] L. Kuzmichev (Baikal Collaboration). The Baikal neutrino experiment: from NT200 
to NT200+.ICRC 2005. 

[197] S. R. Klein. IceCube: A Cubic Kilometer Radiation Detector. arXiv:0807.0034v2 

[physics.ins-det] 

[198] J.Ahrens et al. (IceCube Collaboration). IceCube preliminary design document. 
http://icecube.wisc.edu 

[199] T. Gaisser et al. (IceCube Collaboration). Performance of the IceTop Array. ICRC 
2007 proc, Merida. arXiv:07 11.0353, p. 15-18. 

[200] C. Portello-Roucelle. High Energy Neutrinos from the Cold: Status and Prospects of 
the Icecube Experiment. arXiv:0805.3546 [astro-ph] 

[201] http://www.icecube.wisc.edu/science/ 

[202] A. Achterberg et al. (AMANDA Collaboration). First year performance of the Ice- 
Cube neutrino telescope. Astrop. Phys. 26 (2006) 155-173. 

[203] F. Halzen, I.E. Jacobsen, E. Zas. Ultra-Transparent Antarctic Ice as a Supernova 
Detector. Phys.Rev. D53 (1996) 7359-7361, arXiv:astro-py9512080vl 



High-Energy Astrophysics with Neutrino Telescopes 



85 



[204] A. Achterberg et al. (IceCube Collaboration). Five years of searches for point 
sources of astrophysical neutrinos with the AMANDA-II neutrino telescope. Phys. Rev. 
D75 (2007)102001. arXiv:astro-ph/0611063. 

[205] A.Achterberg et al. (IceCube collaboration). Detection of atmospheric muon neutri- 
nos with the IceCube 9-string detector. Phys. Rev. D 76 (2007) 027101. 

[206] J.L. Bazo Alba. (IceCube collaboration). Search for neutrino point sources with Ice- 
Cube 22-strings. arXiv:0811.4110vl [astro-ph] 

[207] J. Braun et al. Methods for Point Source Analysis in High Energy Neutrino Tele- 
scopes. Astrop. Phis. 29 (2008) 299. 

[208] C. Finley et al. ( IceCube Coll.). Nine-String IceCube Point Source Analysis. 
ICRC2007 proc, Merida. arXiv:071 1.0353. p. 107-110. 

[209] K. Hoshina, J. Hodges, G. C. Hill (IceCube Coll.). Searches for a Diffuse Flux of 
Extra-Terrestrial Muon Neutrinos with AMANDA-II and IceCube. ICRC2007 proc, 
Merida. arXiv: 07 11.0353. p. 67-70. 

[210] A. Achterberg et al. (AMANDA Coll.). Limits on the Muon Flux from Neutralino 
Annihilations at the Center of the Earth with AMANDA. Astrop. Phys. 26 (2006) 129. 

[211] M. Ackermann et al. (AMANDA Coll.). Limits to the muon flux from neutralino 
annihilations in the Sun with the AMANDA detector. Astrop. Phys. 24 (2006) 459. 
arXiv:astroph/0508518. 

[212] J. Alvarez-Muniz, F. Halzen, Astrophys. Possible high-energy neutrinos from 
the cosmic accelerator RX J1713.7-3946. J. 576 (2002) L33-L36. arXiv:astro- 

ph/0205408v3 

[213] H. Landsmann et al. (IceCube Coll.). Radio Detection of GZK Neutrinos - AURA 
Status and Plans . ICRC2007 proc, Merida,arXiv:07 11.0353. p. 163-166. 

[214] For more information and for a complete list of the ANTARES members see: 
http://antares.in2p3.fr 

[215] P. Amram et al. Background light in potential sites for the ANTARES undersea neu- 
trino telescope. Astropart. Phys. 13 (2000) 127-136. 

[216] J. A. Aguilar et al. Transmission of light in deep sea water at the site of the Antares 
neutrino telescope. Astropart. Phys. 23 (2005) 131-155. 

[217] P. Amram et al. Sedimentation and Fouling of Optical Surfaces at the ANTARES 
Site. Astropart Phys. 19 (2003) 253-267. 

[218] P Amram et al. The ANTARES optical module. Nucl. Instr. Meth. A484 (2002) 369. 

[219] M. Ageron et al. Studies of a full scale mechanical prototype line for the ANTARES 
neutrino telescope and tests of a prototype instrument for deep-sea acoustic measure- 
ments. Nucl. Instrum. Meth. A581 (2007) 695-708. 

[220] M. Ageron et al. Performance of the First ANTARES Detector Line. arXiv:08 12.2095 
[astro-ph]. Submitted to Astrop. Phys. 



86 



T. Chiarusi and M. Spurio 



[221] J.A. Aguilar et al. Study of large hemispherical photomultiplier tubes for the 
ANTARES neutrino telescope. Nucl. Instr. Meth. A555 (2005) 132-141. 

[222] J.A. Aguilar et al. he data acquisition system for the ANTARES neutrino telescope. 
Nucl. Instr. Meth. A570 (2007) 107-116. 

[223] J. Aguilar. Analysis of the optical beacon system and search for point-like sources in 
the ANTARES neutrino telescope. PhD thesis December 2007. CSIC, Valencia, Spain. 
http://antares.in2p3.fr/Publications/ 

[224] M. Ageron et al. The Antares optical beacon system. Nucl. Instr. Meth. A578 (2007) 
498-507. 

[225] A. Margiotta (ANTARES Coll.). Systematic uncertainties in MonteCarlo simulations 
of the atmospheric muon flux in the 5-line ANTARES detector. Nucl. Instrum. Meth. 
A602 (2009) 76. arXiv: 0809. 0268 [astro-ph]. 

[226] M. Spurio (ANTARES Coll.). ANTARES neutrino telescope: status, first results and 
sensitivity for the diffuse neutrino flux. 2nd Heidelberg Workshop on High-Energy 
Gamma-rays and Neutrinos from Extra-Galactic Sources, arXiv:0904.3836 [astro-ph] 

[227] N.N.Kalmykov, S.S.Ostapchenko, A.l.Pavlov. Quark-Gluon-String Model and EAS 
Simulation Problems at Ultra-High Energies. Nucl. Phys. B (Proc. Suppl.) 52B (1997) 
17. 

[228] P. Antonioli et al. A three-dimensional code for muon propagation through the rock: 
MUSIC. Astrop. Phys. 7 (1997) 357-368. 

[229] M. Ambrosio et al. (MACRO Coll.). High-energy cosmic ray physics with the 
MACRO detector at gran sasso: part 1. Analysis methods and experimental results. 
Phys. Rev. D56 (1997) 1407. 

[230] M. Ambrosio et al. (MACRO Coll.). High-energy cosmic ray physics with the 
MACRO detector at Gran Sasso: Part 2. Primary spectra and composition. Phys. 
Rev. D56(1997) 1418. 

[231] M. Ambrosio et al. (MACRO Coll.). Measurement of the energy spectrum of under- 
ground muons at Gran Sasso with a transition radiation detector. Astrop. Phys. 10 
(1999) 11. 

[232] M. Ambrosio et al. (MACRO Coll.). Measurement of the residual energy of muons 
in the Gran Sasso underground laboratories. Astrop. Phys. 19 (2003) 313. 

[233] M. Ambrosio et al. (MACRO Coll.). High statistics measurement of the underground 
muon pair separation at Gran Sasso. Phys. Rev. D60 (1999) 032001. 

[234] E. V. Bugaev et al. Atmospheric muon flux at sea level, underground, and underwater. 
Phys. Rev. D58 (1998) 054001. 

[235] J. Horandel. On the knee in the energy spectrum of cosmic rays. Astrop. Phys. 19 
(2003) 193. 

[236] C. Bigongiari (ANTARES Coll.). First results of the ANTARES neutrino telescope, 
J.Phys. Conf. Ser. 173(2009) 012024. 



High-Energy Astrophysics with Neutrino Telescopes 



87 



[237] M. Ambrosio et al. (MACRO Coll.). Neutrino astronomy with the MACRO detector 
Ap. J. 546 (2001) 1038. 

[238] A. Romeyer, R. Bruijn, and J.D. de Zornoza. Muon energy reconstruction in 
ANTARES and its application to the diffuse neutrino flux. Proceedings of the 28th 
International Conference on Cosmic Rays, vol. 3, 2003, pp. 1329-1332. 

[239] A. Capone et al. Status of the NEMO Project. Nucl. Instrum Meth. A602 (2009) 
47-53. 

[240] E. Migneco et al. Recent achievements of the NEMO project. Nucl. Instr. Meth. A588 
(2008) 111. 

[241] R Sapienza. Status of simulations in NEMO, http://www.vlvnt.nl/proceedings, Am- 
sterdam (2003). 

[242] A. Capone et al. Measurements of light transmission in deep sea with the AC9 
trasmissometer. Nucl. Instrum. Meth. A487 (2002) 423. 

[243] F. Ameli et al. The data acquisition and transport design for NEMO phase I. EEE 
Trans. Nucl. Sci. 55 (2008) 233. 

[244] F. Simeone Data taking system for the NEMO experiment. Nucl. Instrum. Methdods 
A588 (2008) 119 . 

[245] I. Amore et al. First results from the NEMO Phase- 1 experiment. Nucl. Instrum Meth. 
A602 (2009) 68-71. 

[246] M. Sedita (NEMO Coll.). Electro-optical cable and power feeding system for the 
NEMO Phase-2 project. Nucl. Instrum. Meth. A567 (2006) 531. 

[247] E.G. Anassontzis et al. The optical module for the NESTOR neutrino telescope. 
Nucl. Instrum. Meth. A479 (2002) 439. 

[248] G. Aggouras et al. A measurement of the cosmic-ray muon flux with a module of the 
NESTOR neutrino telescope. Astrop. Phys. 23 (2005) 377. 

[249] S. Kuch. Design studies for the KMSNeT Neutrino Telescope. PhD thesis 2007, 
University of Erlangen, Germany. FAU-PIl-DISS-07-001. 

[250] J. Carr et al. Conflguration studies for a cubic-kilometer deep sea neutrino telescope 
- KM3NeT- with NESSY, a fast and flexible approach. Proceedings of 30th Int. Cosmic 
Ray Conf. July 2007. Merida, Mexico. 

[251] R.I. Bagduev et al. Nucl. Instr Meth. A320 (1999) 139. 

[252] U. Katz, Status of the KMSNeT Project, Nucl. Instrum Meth. A602 (2009) 40-46. 

[253] R. Davis, Nobel Lecture: A half-century with solar neutrinos. Rev. Mod. Phys. 75 
(2003) 985-994. 

[254] J.N. Bahcall. Neutrino astrophysics. Cambridge University Press. New York, 1989. 



88 



T. Chiarusi and M. Spurio 



[255] J.N. Bahcall. Solar Models and Solar Neutrinos: Current Status. Phys. ScriptaTlll 
(2005) 46-50. (doi: 10. 1088/003 1-8949/2005/T121/006) (Proc Nobel Symposium 129: 
Neutrino Physics; arXiv:hep-ph/0412068). 

[256] B. T. Cleveland et al. (Homestake Coll.), Measurement of the solar electron neutrino 
flux with the Homestake chlorine detector. Astrophys. J. 496 (1998) 505. 

[257] W. Hampel et al. (GALLEX Coll.), GALLEX solar neutrino observations: Results 
for GALLEX IV. Phys. Lett. B447 (1999) 127; M. Altmann et al. (GNO Coll.), Com- 
plete results for five years of GNO solar neutrino observations. Phys. Lett. B616 
(2005) 174. 

[258] J. N. Abdurashitov et al. (SAGE Coll.), Measurement of the solar neutrino cap- 
ture rate with Gallium metal, Part IIL Phys. Rev. C80 (2009) 015807. Also: 
arXiv:090 1.2200. 

[259] J.P. Cravens et al. (Super-Kamiokande Coll.), Solar neutrino measurements in Super- 
Kamiokande-IL Phys. Rev. D78 (2008) 032002. 

[260] A.B. McDonald, Solar Neutrino Measurements. New J. Phys. 6 (2004) 121. Also: 
arXiv:astro-ph/0406253. B. Aharmim et al. (SNO Coll.). Determination of the and 
Total Solar Neutrino Fluxes with the Sudbury Neutrino Observatory Phase I Data 
Set. Phys. Rev. C 75 (2007)045502. 

[261] G. Bellini et al. (BOREXino Coll.), Measurement of the solar neutrino flux with 
246 live days of Borexino and observation of the MSW vacuum-matter transition. 
arXiv:0808.2868. 

[262] M. Koshiba. Observational neutrino astrophysics. Phys. Rep. 220 (1992) 229-381. 

[263] G. Raffelt. Supernova neutrino oscillations. Phys. Scripta T121 (2005) 102-105. 
(doi:10.1088/0031-8949/2005/T121/014) (Proc Nobel Symposium 129: Neutrino 
Physics; arXiv:hep-ph/0501049). 



